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Average H—H, Repulsive Potential from Collision Cross Section Measurements* 
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(Received December 9, 1942) 


The average repulsive potential between a hydrogen atom and a hydrogen molecule has been 
computed from experimental values of elastic collision cross sections for fast hydrogen atoms 
(with energies of 200 to 800 electron volts) scattered in room temperature molecular hydrogen. 
It was possible to solve the classical theory scattering equations explicitly for the repulsive 
potential by proper analytical representation of the experimental cross section values as a 
function of the relative kinetic energy of the scattering system. The potential function so ob- 
tained predicts total collision cross sections within the experimental error and has the form 
V(r) =A exp (—ar?)+B exp (—bdr?) where r, the distance between the mass points representing 
the hydrogen atom and the hydrogen molecule, extends, in the present case, from 0.28 to 0.70A. 
For 7 in A units, A =0.846X 10-" erg, a=24.9A~2, B=0.211 107" erg and b=2.40A~. 





HE repulsive potential between two helium 
atoms at interaction distances of 0.55 to 
1.05A has been estimated! from measured values 
of total collision cross sections of high velocity 
helium atoms elastically scattered by room tem- 
perature helium. Since experimental values are 
available? for total collision cross sections of high 
velocity hydrogen atoms scattered by room tem- 
perature hydrogen, it should be possible to 
estimate the average repulsive potential between 
a hydrogen atom and a hydrogen molecule at the 
separation distances involved in the scattering, 
0.28 to 0.70A. 
The original curve of H — H: cross sections as a 
function of the hydrogen atom energy indicated 


* Contribution No. 470 from the Research Laboratory 
of Physical Chemistry, Massachusetts Institute of Tech- 
nology. Presented before the Division of Physical and 
Inorganic Chemistry of the American Chemical Society 
at the Buffalo Meeting, September 8, 1942. 

— and H. Pearlman, J. Chem. Phys. 9, 503 

* 1. Amdur and H. Pearlman, J. Chem. Phys. 8, 7 (1940). 


the presence of some inelastic scattering for beam 
particle energies between 500 and 700 volts. 
Since only elastic collision cross sections can be 
used in the present calculation, it is necessary to 
draw the best smooth curve through the original 
experimental points in order to read off the 
appropriate elastic collision cross sections in the 
range 500-700 volts. The values so obtained are 
listed in Table I where the original experimental 
values which show the effects of inelastic scatter- 
ing at 600 and 700 volts appear in parentheses. 


TABLE I. Total elastic collision cross sections for high 
velocity hydrogen atoms scattered by room temperature 
hydrogen. 








Cross section, Sg 
(10-16 cm?) 


1.455 (1.455) 
0.618 (0.618) 
0.410 (0.410) 
0.310 (0.310) 
0.268 (0.302) 
0.243 (0.250) 
0.230 (0.230) 


Percent error 
Voltage, V + 











157 
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The estimated percent errors are those originally 
recorded for the values in parentheses and shall 
be taken to represent the probable reproduci- 
bility of the total elastic collision cross sections, 
since only the 600-volt value shows appreciable 
correction for inelastic scattering. 

It has been shown? that the elastic scattering 
of the high velocity hydrogen atoms may be 
treated by classical methods since the average 
relative angular aperture of the detecting system, 
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Fic. 1. Analytical representation of total elastic collision 
cross section as a function of beam particle energy. 


3.37° (2.25° in apparatus coordinates) far exceeds 
the maximum relative angular aperture, 1.15° to 
1.72°, required for detection of effects due to 
quantum theory diffraction.® ; 
The general equations and method of pro- 


cedure for estimating repulsive potentials of the — 


form V(r) =A exp (—ar”) have been discussed in 
connection with the calculations for helium.! The 
reduced mass yp and the relative angle of scatter- 
ing 8 appearing in the equations are equal, in 
the present case, to two-thirds the mass of the 
beam particle and to three-halves the angular 
deflection of the beam particle with respect to 
the beam axis, respectively. 

Attempts to reproduce the tabulated cross sec- 
tions by variation of the constants, A, a, and n, 
in the general potential V(r)=A exp (—ar’) 
were unsuccessful. Accordingly, the potential 
V(r) =A exp (—ar")+B exp (—dr") was investi- 
gated. For n=2, the relative angle of scattering 


3H. S. W. Massey and C. B. O. Mohr, Proc. Roy. Soc. 
A141, 434 (1933). 
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is given in closed form by 


2ro ¢*® LV (1%) — V(r) Jrdr 2(am)'roA 


uv? Jr 





(r?—19?)3 uv? 


2(br)'roB 
Xexp (—aro”?) +-————— exp (— dro”), (1) 


pv 


‘where 79 is the distance of closest approach of the 


beam particle to the scatterer and 2, the relative 
initial velocity of the colliding system. Equation 
(1) assumes that the beam particle and the 
scatterer are mass points in a spherically sym- 
metrical force field. 

It is possible to express 7o in terms of the col- 
lision cross section Se since, for a given voltage 
and average angular aperture, 


V(ro) 





So= | i fe = be’, 


uv? 
where be is the impact parameter defined by the 
above equation. From attempts to find a satis- 
factory single term potential function, it is esti- 
mated that [1—V(ro)/3uv?] will have a value 
between 0.96 and 0.94 for voltages between 200 
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Fic. 2. Average repulsive potential between a hydrogen 
atom and a hydrogen molecule considered as mass points. 


and 800 volts, so that to a first approximation 
So = 0.95270. 


The accuracy of this approximation can be deter- 
mined after the potential function has been 
found, by computing values of [1— V(ro)/3y0* | 
for the 7o distances of interest. Since 0 is equal 
to 0.0589 radian, Eq. (1) reduces to 


pv? = 34.8a'A Sot exp (—0.335aSe) 
+34.8bIBSo! exp (—0.335bSe). (2) 
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If, then, it is possible to represent the relative 
kinetic energy of the system, 3uv’, as a function 
of the total collision cross section, within the 
estimated error, by an equation of the form 


pv?/So'= A’ exp (—a’So) +B’ exp (—b’Soe), (3) 


the problem of finding a satisfactory repulsive 
potential is solved since the constants A, B, a, 
and 6 in the two-term potential function may be 
directly computed from the numerical values of 
A’, B’, a’, and b’ of Eq. (3). 

It was possible by graphical methods to evalu- 
ate the constants of Eq. (3) and thus to represent 
analytically within the average estimated error 
of +2.8 percent, the elastic cross sections of 
Table I by the equation 


a 146 8.255, 
358 sat exp (—8.25Se) 


10-"° erg 
+11.3 exp (AK (4) 


where é is the charge on the electron. 

Figure 1 shows that the individual Sg values 
scatter well about the smooth curve represented 
by Eq. (4). 





Comparison of Eqs. (2) and (4) yields directly 
values for the constants A, B, a, and 6 so that 
the first approximation to the desired potential 
energy function is 


V(r) =[0.846 exp (—24.6r?) 
+0.211 exp (—2.37r?) ]K10-"“ erg. (5) 


If the Se values of Table I are divided by 
0.957 a series of first approximation 7? values are 
obtained for substitution into Eq. (5). From 
the V(ro) values so obtained it is found that 
[1— V(ro)/3uv? ] varies between 0.969 and 0.958 
so that the second approximation to the general 
expression for Sg becomes 


Soe = 0.963219? 


and the corresponding improved expression for 
V(r) 


V(r) =[0.846 exp (— 24.97?) 

+0.211 exp (—2.40r?) ]X10-' erg. (6) 
Further approximations produce insignificant 
improvement so that Eq. (6), which is repre- 
sented graphically in Fig. 2, is the final repulsive 


potential energy function derivable from the 
experimental collision cross sections. 
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The vibration-rotation energies of the linear X-Y-Z type molecule are derived according 
to a method suggested by Professor Harald H. Nielsen in a recent paper. By the proper sub- 
stitutions the energy expressions are reduced to the ones appropriate for the XY2 type 
molecule, and a comparison with the work of Professor D. M. Dennison is effected. 





I. INTRODUCTION 


N a recent paper! Harald H. Nielsen showed 

that the vibration-rotation energies of a 
molecule composed of N atoms could be expressed 
to second order of approximation in perfectly 
general terms. Following his suggestion that it 
would be simpler to arrive at the energies of a 
particular model by the general method than to 
derive them independently, a solution of the 
linear X—-Y-Z type, of which HCN and COS are 
examples, is given here. 

The proper modifications are made to trans- 
form these energy expressions to the XY2 type 
and these are compared with similar expressions 
given by Dennison? for this type of molecule. 


Il. THE NORMAL COORDINATES 


In Fig. 1 is given a diagram of the linear 
X-Y-Z type molecule. The coordinate system is 
fixed in the molecule with the origin at the center 
of gravity of the system. The particles are 
numbered from left to right. In the positional 
coordinates x;, y:, 2; the kinetic energy of vibra- 
tion is 


3 
2T=> m(é2+92+2,?). 


t=1 


The condition that the center of gravity of the 
system remain at rest in these coordinates is 


3 
Zz. m;a;=0 
i=1 


(where a; takes the values x, y, z) and the con- 
dition that to zero order the internal angular 


1 Harald H. Nielsen, Phys. Rev. 60, 794 (1941). 
2D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 


momentum is zero may be written 
. 
dX m;(a;°68;—B;°da;) =0 
i=l 


(a and @ taking values x, y, 2 but a#8). The 
superscript (°) indicates the equilibrium value of 
the coordinates. 

As intermediate coordinates in the deter- 
mination of the normal coordinates let 


30x /2>0 =x2— (mix1+M3x3)/2me, 
3ay/20 =yo— (miyitmsys)/2m2, 
302/2>- =Ze— (m42;+m323)/2mo2, 


gd =21— 23. 


(1) 


The coordinates x;, y:, 2; can be expressed in 


Paul 
/ 21-3 
-~----@- --e ¥--@------ —Z 


Fic. 1. The coordinate system of the linear 
X-Y-~-Z type molecule. 


---------% 


terms of x, y, 2, and q as follows: 


X1 = — Myo(Z2°—23°)x/m)>._ (21° — 23°), 
x2=0x/>i, 

X3= —Mea(219— 22") x/my,>_ (21°—23°), 

Y1= —Me0 (22° —23°)y/m1y (21° — 23°), 
ye=oy/), (2) 


Y3= —mMe20(219—22°)y/m1> (21°— 23°), 
21= —m22/>.+m3q/e, 

Z2.=02/>., 

23= —m22/>,—migq/e, 
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where 
o=m,+MmM3, > =M,+M2+MmMs3. 
The kinetic energy may now be written in these 
coordinates 
27 = pig? +ue(e?+y") +432? 
where 
M1=Mm3/o,  b2= mso?I,,/mym3>, (2;°—2;°), 
(3) 
w3=moa/)>.. 
In the displacement coordinates the potential 
energy is 
2V=K,iQ01?+ K2Q2?+ K3Q03?+2K 120102 
+2K130103+2K230203+Ka(x?+y’). (4) 


Let Q:=22—21, Q2=23—22, Q3=23—21:. Making 
these substitutions into (4) and utilizing the 
relations (2), the potential energy may be written 
2V=Rig? +R 32+ 2Rs2q+ Ka(x? +") 

in which 
Ri = (Kym;3?/0?+ Kem,?/0?+ Kz 

+2Ki2mym3/o?+2Kism3/o+2Ke23m1/a), 
R3= (Kit+K.—2Kw), (S) 
R= (—Kym3/0o+K2m,/o 

— Ky2(m,—m3)/o— Ki3+Kaz). 

Now let (u:)'g=qi', (u2)*x=qo1, (u2)*y=qoe, 


(u3)*z=q3'. In these coordinates the kinetic and 
potential energies become 


2T = Gr'?+ (Gar +402) +49", ” 

2 V =hi'qu'?-+ he! (gor +922) +hs'qs’?-+2ha'q1'gs’; 
where . 

ky’ =R1/m1, ko! = K4/ po, ks’ =R3/us, 
ka’ =R4/(urms)!. 

If gi’=qisin y+qs cosy and q;’=—gqi cos y 
+q3sin y are substituted into (6) and sin y 
= (3)8{1+[(kr’ —ks’)?/(4ka? + (ki’—ks’)? ]}* the 
cross product in q1’q3’ vanishes and qi, qa1, G22, Ys 


are the normal coordinates. The resulting form 
for the kinetic and potential energies is 


2T=q:°+ (Gort G22) +43’, 
2V= kigi?-++Ro(qo1-+922) +ksq3?, 


LINEAR X-Y-Z TYPE 





MOLECULE 


in which 

ki =k,’ sin? y+k;’ cos? y—2k,’ sin y cos y, 
ko=k?’, (7) 
k3=ky’ cos? y+k;’' sin? y+2k,’ sin y cos y. 


The displacement coordinates can be expressed 
as linear combinations of the normal coordinates 
as follows: 


X= LiseGso, Vi= LM izoGso, i= LN iseGse, 


where 7 is summed over the particles, s is summed 
over the frequencies, and o takes the values 1, or 
1,2 depending on whether or not the frequency 
is a degenerate one. An example of these sum- 
mations is X;=1)11911 +1 i21G21 +1 i22922 +1 31931. The 
only coefficients having values other than zero 
are 


Lie1=M122= — [moms(22—25°)?/mid Tez}, 
° e 

L201 = M222= [mym3(21°— 23°) 2/me2>_ 122 |', 

L321 =M322= — [myme(21°—22°)?/ms>_ } pay 


111 =[(m2/a_2)* cos y+(m3/om)' sin y ], 

1131= —[(m2/o>_)' sin y—(m3/am,)' cos y ], (8) 
No11= —(o/m2)_)' cos ¥, 

N231= (¢/me2),)' sin y, 

2311 = [.(m2/o>_.)! cos y—(mi/oms)! sin y |, 

N331= —[(m2/o>.)' sin y+(m1/oms)' cos y ]. 


With the aid of the Eqs. (8) the coefficients a,, 
and Ages’e" OCCurring in the expressions for the 


. . ° , , 
instantaneous moments of inertia J,,=J,, may 
be determined. The values of these coefficients 
follow: 


41 = —2(Iez) *¥2081, d21=A22=0, 
31 = 2(I¢2) E2011, 
2 
Ajin1=1— £2011, A 2121 =A 2222=0, 


A sisi =1- £0031, (9) 
where 


f2011= { (myms/aIz2)*(21° — 29°) cos 
+-(m2>./oIz2)'22° sin y}, 
foos1 = { — (mym3/oIz2)*(21°—23°) sin y 


+(m2>/oIz2)*z2° cos 7}. 


In the normal coordinates the first- and second- 
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order potential energy expression may be written 


Vi=inigi?+1i3g1°¢s +P 1229p" 
+k1339193" + o28p*@st+Rs339s°, 
Vo=harii1g1t+Rii1391°9gs+F112291"p” 
+ k1133917¢s?+h1ss39193° + k2222p" 
+ ko233p7Gs"+R:s33393' (10) 
in which 
2 2 
p?= (g21 +22). 


With the constants derived the explicit ex- 
pressions for Ho, H; and H2 may now easily be 
written. 


Ill. THE VIBRATION-ROTATION 
ENERGY EXPRESSIONS 


The vibration-rotation energies may be con- 
veniently written as 


(E/he) = (Evin/he) + (Erot/he), 
where (Eyin/hc) = 
(Eo/he) +01(01 +3) +o2(v2+1) +o3(vs+3) 
+211(01 +3)? +x22(v2+1)?+433(03 +3)? 
+ (x12+%21) (01 +3) (0241) 
+ (x13+%31) (01 +3) (v3 +2) 
+ (x23+232)(v2+1)(vst+3)+xl2t!2’, (11) 
where 
x11 = (h/32e4c%w;2) { 3kars 
— 15k 111/84 %C%01?— hi1s/40°C%05? 
+hi13/8n%c*(4e?—ws") }, 
X22 = (h/32m'c?w*) { 3ko200 
—his2/4n°C%w1?—ho09/40°C%ws? 
+ hi22/89%C?(4e2?— ws?) 
+ k223/8m%c?(4ers?— ws?) }, 
X33 = (h/322'c8ws*) { 3k3333 
_ 15k333/89°C%3? _ kiss/44°c%o1? 
+hiss/8m°c?(4us?—w1?)}, 
X12 = (h/322*c8wywe) { Rir2e 
— 6Ri11k122/497C*w? 
—ki13k223/4m°c*w3” 


+49C?(w?/Toz) foo11}, 
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X01 = (h/322'c%w we) { R1122 — Roosku3/42°c?ws" 
— Ri2e/°c?(4eos? —w”) 
+49c?(w1?/Iz2) Esau}, 11) 
X13 = (h/322'c®wyws) { R1133 — OR111k133/49*c7w1? 
= bi13/mc?(4en1? —ws?) I, 
X31 = (h/322'*c%wiws) { R1133 — OR 333k 113/49°C7ws” 
= kis3/°C?(4u3?— 1?) - 
X23= (h/322'c?wows) { Rooss = Ry22k133/407c*w1" 
— koos/°c?(Aeoe? —w3") 
+4 nc? (wo3?/T sz) £2231}, 
X32 = (h/322*c?wows) { Ro233 — ORs33h203/4m2c?ws" 
—ki3sk122/49°c*w," 
+-4¢?(w2/Te2) foos1}, 
Xlgle = (—h/32r'c®we?) { Reo20 
+ ki22/8x°c?(4w2?—w1”) 
+ ho25/80%C?(4u2" —ws*)}, 
and 
(Exor/he) =[J(J +1) -2]B, 
—[J(J+1)-P]}D (12) 
in which 
B,=B,.— (vi +})a1— (¥2+1)a2—(vs+3)as, 
B,=h/8n*cIz2, 
a= B.(h/49°cws1 22) {1— Bi ceses*/ (w1? — we?) 
- Atoos+ 6(Izz)*Es2s1k111/40°C%w1? 
+2 (Ize) *€2211k113/40°C%w3"},, 
as= B.(h/8x*cw2I 22) { e 4(Iz2) *£0011k223/4m7C7w3" 
+4(Iz2)*£o0s1k129/49°C%w1? 
ti (3w2?+w1”) Eoo11/ (we? ~ 1”) 
am (3w2?+w3") E0s1/ (ws? —ws?) ’ (13) 
ag= B.(h/4x°cws] vz) { i— Atsos1ws?/ (ws? — we”) 
on Agee — 6(Iz2)8Es211k333/40°C%ws? 
+2(Iee)*€s2stk133/44°c%an"} , 
D=(h?/1289°c3(Iz2)*) { Eoo1/eos?-+ Eoo1/w1?} ’ 
l=]. 
IV. REDUCTION TO THE XY, TYPE MODEL 


The accuracy of the above expressions may be 
tested by reducing them to the energies for the 
XY2 type model and effecting a comparison with 
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these results? which are well known to be correct. 
The following simplifications may be made: 


m,=mM3, ¢)°= — 23°, Z°=0, sin y= 1, 
cosy=0, Ki=Ke, Kis=Ko;, 
Ri =piki’=0, £5011=0, €oo31= —1, 


P (14) 
a= —2(Iz2)', 021 =422=431=0, 


Ain=1, A131 =A2121=A 2222 =0, 
ki13= koo3 = R333 = 0. 
When these substitutions are made in (11) and 
(13), these expressions reduce to 
211 = (h/32m4c%12) { 31111 — 15k 111/84 2C%17}, 
Xo2 = (h/32m'c%we?) { 3Ro200 
—ki22(8a2? — 301”) /8x°c%e1"(4u02*—w1") |, 
x33 = (h/32m'c*ws") { 3k3333 
— kis3(8ws? — 3w12) /8°e%w12(4es?—w1?) }, 
X12 = (h/322'c®w we) 
X { Ritee— 6Ri11k122/42%c7w1"}, 
Xo1 = (h/322'c®wiwe) 
x { R1122— he120/°C?(4u0s?— ws?) }, 
X13 = (h/322'c®wiws) 
X { R1133— 62 111k 133/42c?w1"}, 
X31 = (h/322'c3 ww) 
x {Rissa — hiss/m°c?(4es?— wr?) }, 
X23 = (h/322'c?wows) 
X {Ree33— Riveki33/49°c7w1* 
+49°c?(w3*/T22)}, (15) 
X32= (h/322'c? wows) 
X {Re233— Riveki33/49°C*w 1? 
+4n%c?(ws*/Iz2)}, 
Xlgle = (—h/322'c8 wo?) 
x { koa00-+hi20/8m°C%(4es?—wo1) }, 
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a= —B,(h/44°cwT 22) 
x {3+6(1;2)*Ri/49°c7w1'}, 
a= —B,(h/8x°cwsl +2) 
X {+4(22) i22/4°c%w1”) 
+ (3w2?-+w3*) /(w2?—ws")}, 
a3= —B,(h/4n°cwsIz2) 
X {2(Tez)'133/49°c%w1? 
+ (3w3?-+- wo?) / (ws? it we”) } ’ 
D= h?/128%c3(I..2) 3w1?, 
B.=h/8n°cI-», 
B,=B,.— (vit })ai— (vo+1)a3— (vst+4)as. 
When it is observed that the relation between 


Dennison’s coordinates o, p, ¢ and those given 
in this discussion, q1, p, 93, is 

o=2ncw1'qi/(hc)?, p=2xcws'p/(hc)}, 
¢ = 2mcws!qs/ (hc)? 


it can easily be shown that the constants (a: - -7) 
in his V; and V2 take the following values in 
terms Of Rese) ANd Reser e’’! 


@=ky31(he)*/82*c8w;}, 
b=ky22(hc)*/82c3w wo, 
C= k133(hc)*/82%c3w $s, 
d=kyyh/167'c8w,?, 


€ = Ro202h/162'c%wo’, 
f=ks333h/162'c8ws”, 
2=Ry22h/1679'c8wwo, 
h=ki133h/1679'c®wws, 
4 = koo33h/162'c? wows . 


The substitution of these constants into Eqs. (15) 
cause them to reduce exactly to the constants x;; 
given by Dennison. 

The author expresses his gratitude to Professor 
Harald H. Nielsen, of the Ohio State University, 
who suggested the problem and with whom 
valuable discussions were held during its solution. 

3An exception may be noted in that Dennison has 
neglected to calculate the term D due to centrifugal 
stretching of the molecules. It should also be noted that 


the term /?B, occurring in Dennison’s x,; has here been 
included in the Erot. 
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The Quenching of Cadmium Resonance Radiation (3261A) by 
Hydrocarbons and Other Gases* 


E. W. R. STEACIE AND D. J. LERoy 
Division of Chemistry, National Research Laboratories, Ottawa, Canada 


(Received February 10, 1943) 


An investigation has been made of the quenching of cadmium resonance radiation (3261A) 
by a number of hydrocarbons and other gases. The following values of the quenching cross 


sections were obtained (in cm? X 10!): 


He 3.54 1-C4Hs 

D2 1.80 2-C4Hs 

NHs 0.052 CoHe 

CH, 24.9 C2He 

C;He 29.1 cyclopropane 


A few rough measurements have also been made with 2288. 





INTRODUCTION 


N the last few years we have investigated a 
series of cadmium and zinc photosensitized 
reactions of hydrocarbons,'~ using the cadmium 
resonance lines 3261A and 2288A, and the zinc 
lines 3076A and 2139A. For a proper interpreta- 
tion of the results obtained it is desirable to have 
information on the quenching of these lines by 
the gases employed. Since such data do not 
exist, it was decided to investigate the quenching 
of cadmium and zinc resonance radiation by a 
number of gases. The present paper is primarily 
concerned with the quenching of the 3261A 
cadmium line (Cd 5°P,;—Cd 51S). 

The only previous work on the quenching of 
\3261 is that of Bates,> Bender,® and Lipson and 
Mitchell.? Bates and Bender made purely quali- 
tative observations of the quenching by hy- 
drogen. Lipson and Mitchell made a quantitative 
investigation with a number of gases and found 
the following quenching cross sections: 


He 0.67 X 10-6 cm? NH; 0.041 X 10-"* cm? 
De 0.19 Ne 0.021 
CH, 0.012. 


Co 0.14 


* Contribution No. 1104 from the National Research 
Laboratories, Ottawa, Canada. 
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Phys. 9, 306 (1941). 

2 E.W.R. Steacie, Annals N. Y. Acad. Sci. 41, 187 (1941). 

3E. W. R. Steacie and D. J. LeRoy, J. Chem. Phys. 10, 
22 (1942). 

4H. Habeeb, D. J. LeRoy and E. W. R. Steacie, J. 
Chem. Phys. 10, 261 (1942). 

5 J. R. Bates, Proc. Nat. Acad. Sci. 14, 849 (1928). 

6 P. Bender, Phys. Rev. 36, 1535 (1930). 

7H. C. Lipson and A. C. G. Mitchell, Phys. Rev. 48, 
625 (1935). 


35.2 C:He¢ 

30.6. C;Hs 

28.4 n-C4Hi0 Ca. 0.02 

22.0 iso-C4H 19 

0.71 
EXPERIMENTAL 


The experimental method was in principle the 
same as that used by Lipson and Mitchell. Two 
quartz cells were employed, one of which con- 
tained cadmium vapor alone, while the other 
contained a foreign gas as well. A beam of \3261 
was sent through both cells, and the resonance 
radiation emitted at right angles to the incident 
beam was photographed. The intensities of the 
resonance radiation from the two cells were then 
compared by photographic photometry. 

The experiments were performed in such a way 
that the light path through the cadmium vapor 
was a minimum, the pressure of cadmium vapor 
was very low, and the pressure of the foreign gas 
was as low as possible. Under these conditions 
the Stern-Volmer equation holds, and the quench- 
ing cross sections can be calculated from the 
relation® 


1 2x#N/ 1 1 4 
rot=— / 2666.64] —+—)| 
QO-p kT \M, Mo 


where o@? is the quenching collision cross section 
in cm?, Q is the quenching, i.e. the ratio of the 
resonance radiation emitted in the presence of a 
foreign gas to that which would be emitted in its 
absence, 1/(Q-p) is the slope of the curve ob- 
tained by plotting 1/Q against the pressure of the 
foreign gas in mm, ¢ is the mean life of the 5°P; 
state of cadmium, which was taken to be 

8 A. C. G. Mitchell and M. Zemansky, Resonance Radia- 


- a Excited Atoms (The Macmillan Co., New York, 
1934). 
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2.5X10-* sec.,8 N is Avogadro’s number, F is 
the gas constant per molecule, and M,; and M, 
are the molecular weights of the colliding par- 
ticles. 


The Light Source 


The source of cadmium resonance radiation 
was of a type which has been previously de- 
scribed.* The emitting part of the lamp was 
made of fused quartz. To the ends of this, 
through graded seals, were attached Pyrex elec- 
trode chambers containing ‘‘coated”’ iron elec- 
trodes. The lamp contained a small piece of 
cadmium and was filled with neon at a pressure 
of 3 mm at room temperature. It was operated 
from a sign transformer at 120 milliamperes. 
The lamp, including the electrode chambers, 
was contained in an electric furnace at 280°C, 
the radiation emerging through a fused quartz 
window. 

The radiation from the lamp consisted mainly 
of the resonance lines at 2288 and 3261A. The 
intensity of these lines was very high, and the 
output of the lamp was unusually steady, the 
maximum variation being less than 5 percent 
over long periods of time. 


The Optical System 


The radiation from the lamp passed through a 
slit, and through a Corex D plate which removed 
\2288 completely, and was then focused into a 
convergent beam by two quartz lenses. The beam 
then passed a shutter and entered the furnace 
containing the resonance cells through a quartz 
window. The resonance cells were situated one 
above the other, and the image of the source 
gave a narrow line which was long enough to 
cover both cells. Opposite each cell was a slit 
which allowed a narrow pencil of radiation to 
enter the cell. 

The fluorescent radiation was taken off at 
right angles to the incident beam through another 
slit. The incident and emergent slits were situated 
as close as possible to a common edge of the reso- 
nance cells in order to provide the shortest 
possible light path through cadmium vapor, and 
thus diminish effects due to imprisonment. 

The beams emerging from the cells passed out 
of the furnace through a quartz window and 
entered the camera where they were focused on 


an Eastman 40 plate. The method of photometry, 
etc. will be described later. 


The Resonance Cells, etc. 


The resonance cells were made of fused quartz, 
with two polished windows about 2 cm in diam- 
eter meeting at a common edge. For the dual 
purpose of maintaining uniformity of tempera- 
ture, and the elimination of scattered light, the 
cells were mounted above one another in a metal 
box with heavy walls, provided with slits for the 
incident and fluorescent beams. The cells were of 
the usual “‘light-trap’’ type. 

Short side tubes of wide bore (14 mm) passed 
out from each cell into holes in a single large 
brass block. Each side tube contained small 
pieces of cadmium. The wide, short tubes per- 
mitted rapid diffusion of the vapor, and the brass 
block kept the vapor pressure of cadmium the 
same in both cells. 

A further side tube from each cell led through 
a graded quartz-to-Pyrex seal to a magnetic 
cut-off to prevent diffusion of cadmium vapor to 
the cold tubing outside the furnace. The tubing 
on the far side of the cut-off passed out of the 
furnace and was connected through a stopcock 
with a vacuum system consisting of pumps, 
McLeod gauges, gas reservoirs, etc. 

The whole system up to the far’side of the 
magnetic cut-offs was contained in an electric 
furnace. The distribution of the heating elements 
in the furnace was such that the brass block was 
colder than the rest of the system, and its tem- 
perature could be independently controlled. The 
windows of the resonance cells were kept at a 
slightly higher temperature than the rest of the 
system to prevent the condensation on them of 
cadmium vapor. Five thermocouples in various 
parts of the system, and independent electrical 
circuits for the various heating elements enabled 
the temperature to be accurately controlled to 
within 1° throughout the furnace, and to within 
0.5° or better for the brass block containing the 
side tubes with the cadmium supply. 

The reason for the somewhat elaborate tem- 
perature control system is that the main source 
of error in the measurements arises from varia- 
tions in the cadmium vapor pressure. The in- 
tensity of the fluorescent radiation is very sensi- 
tive to such variations, and the function of the 
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Fic. 1. Quenching of \3261 by hydrogen and deuterium. 


comparison cell is much more to compensate for 
fluctuations in cadmium vapor pressure than for 
fluctuations in the incident light intensity. By 
maintaining the cadmium supplies for the two 
cells at the same temperature equality of vapor 
pressure is achieved in principle. However, in 
operation one cell contains a foreign gas and the 
other does not. If the fluctuations in the tem- 
perature of the cadmium supply are large, or if 
the tubing leading from it to the cells is long and 
narrow, diffusion of cadmium vapor will be much 
slower in the cell containing the foreign gas and 
this may give rise to large errors. For this reason 
the fluctuations in temperature were reduced as 
much as possible, and the connecting tubing 
was short and wide. 


Experimental Procedure 


The procedure in carrying out an experiment 
was as follows. Both cells were evacuated to at 
least 10-* mm. They were then closed off and 
allowed to stand for an hour or more to ensure 
cadmium vapor pressure equilibrium. The shutter 
between the lamp and the cells was then opened, 
and exposures of various lengths of time, from 
about 10 sec. to 100 sec., were made with the 
camera. These served as calibration marks on the 
plate. Trial exposures made with calibrated 
screens showed that the failure of the reciprocity 
law was appreciable (about 4 times the normal 
experimental error), and a correction for this 
failure was therefore made when translating plate 
blackenings into intensities. 

After the calibration marks had been put on 
the plates, the gas under investigation was added 
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to one cell, its pressure measured with a McLeod 
gauge, and the system was then shut off and 
allowed to stand for about an hour. During this 
time the shutter between the lamp and the cell 
was kept closed to avoid possible photochemical 
reactions of the added gas. The radiation from 
the cells was then photographed again, the pres- 
sure of the added gas was altered, and the 
system was shut off and allowed to stand again. 
In this way a series of 4 or 5 pressures were in- 
vestigated in a single experiment. 

The plate was then developed and the blacken- 
ings determined with a microphotometer. The 
light source was very steady and the fluctuations 
in intensity were negligible over the short period 
of time required for the calibration exposures. 
Over the long period required for a whole experi- 
ment some fluctuations occurred, but these were 
compensated for by comparing the blackening of 
the spot from the cell containing the foreign gas 
with the blackening from the evacuated cell. By 
means of the calibration marks blackenings were 
translated into intensities, and the value of 1/Q 
determined by dividing the intensity of the 
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Fre. 2. Quenching of \3261 by ammonia. 


fluorescence in the absence of a foreign gas by 
that in the presence of a foreign gas. 


Materials 


The quenching efficiencies of the gases investi- 
gated varied widely. In the case of substances 
such as the olefins which quenched very strongly 
the presence of traces of impurities could have no 
appreciable effect on the results. The paraffins, 
however, quench only very slightly, and it was 
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therefore of great importance to free them as far 
as possible from traces of gases which quench 
strongly. 

Hydrogen was taken from a commercial cyl- 
inder of ‘“‘pure’’ gas and passed over platinized 
asbestos at 600°C, and through a liquid-air trap. 
In some experiments it was stored over phos- 
phorus pentoxide and in others it was not. 
A number of preparations gave identical results. 
Further, a sample taken direct from the cylinder 
differed only slightly from the purified samples. 

Deuterium was obtained in a cylinder from the 
Stuart Oxygen Company, and was stated to be 
at least 99.5 percent D2. It was purified by the 
procedure used for hydrogen. 

Ammonia was taken from a cylinder and 
purified by an extended series of bulb-to-bulb 
distillations. Ammonia taken directly from the 
cylinder gave almost identical results. 

Ethylene was obtained in cylinders from the 
Ohio Chemical and Manufacturing Company, 
and was stated to be 99.5 percent ‘C2H,. It was 
purified by distillation in a still of the Pod- 
bielniak type. 

Propylene was also 99.5 percent grade, and was 
fractionated similarly. 

Butene-1 and butene-2 were prepared by the 
methods given by Kistiakowsky, Ruhoff, Smith 
and Vaughan.® The resulting gases were frac- 
tionated in the Podbielniak column. 

Cyclopropane was Ohio 99.5 percent grade. 
Since propylene was a possible impurity and 
quenches much more strongly than cyclopropane, 
the purification procedure was designed mainly 


TABLE I. The quenching of 43261. 








Slope of quenching Quenching cross section 





Gas curve mm cm X1016 
He 12.65 3.54 
D2 4.60 1.80 
NH; 0.072 0.052 
CH, 26.3 24.9 
C;He 26.5 29.1 
1-C,Hs 29.0 33.2 
2-C.Hs 25.0 30.6 
CeHe 20.9 28.4 
C:H: 24.0 22.0 
cyclopropane 0.64 0.71 
CH. (0.025) (0.024) 
C;Hs (0.011) (0.012) 
n-C.Hio (0.052) (0.064) 
iso-CyHio (0.037) (0.046) 








°G. B. Kistiakowsky, J. R. Ruhoff, H. A. Smith, and 
W. E. Vaughan, J. Am. Chem. Soc. 57, 876 (1935). 
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for its removal. Cyclopropane from the cylinder 
was shaken for 24 hours with alkaline permanga- 
nate solution, and was then distilled at a high 
reflux ratio in the column. A further sample was 
taken directly from the cylinder and was merely 
given a bulb-to-bulb distillation. This gave re- 
sults identical with those from the preceding 
sample, and it may therefore be concluded that 
no significant amount of strongly quenching im- 
purity was present. 
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Fic. 3. Quenching of \3261 by the paraffins. 


Benzene of American Chemical Standard grade 
was used without purification other than out- 
gassing. It quenched so strongly that impurities 
could have no appreciable effect on the results. 

Acetylene also quenched very strongly, and 
hence a bulb-to-bulb distillation of the com- 
mercial product was a sufficient purification. 

Ethane proved to have a quenching efficiency 
only about 1/1000 that of ethylene, which is one 
of the most likely impurities. Hence, assuming 
the experimental error to be 5 percent, it is 
necessary to reduce the amount of ethylene in 
the ethane to 1/20,000. Ohio cylinder ethane 
which had previously been given a Podbielniak 
distillation was saturated with bromine vapor 
and passed slowly through a vessel illuminated 
with a 1000-watt photo-flood lamp. The gas was 
then passed through a 40 percent potassium 
hydroxide solution, a dry-ice trap to remove 
water, and the ethane was collected in a liquid- 
air trap. The ethane was then again distilled 
through the column. 

Propane from a cylinder of 99.9 percent grade 
was given a treatment similar to that used for 
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Fic. 4. Quenching of \3261 by the olefins. 


ethane. Since the relatively high boiling point 
made it more difficult to remove water, the gas 
was stored in a reservoir over phosphorus pent- 
oxide. 

n- and iso-butane of 99 percent grade were 
treated as described for propane. 


RESULTS 


Preliminary experiments with hydrogen and a 
number of other gases showed that the quenching 
curve was independent of the intensity of the 
incident radiation: this was reduced by a factor 
of about 10 for these experiments by introducing 
a thin soft glass plate into the light path. The 
results were also shown to be independent of the 
position of the entering slit, which was moved 
from its normal grazing position to a distance of 
about 2 mm from the common edge of the cell. 
This indicates that imprisonment of resonance 
radiation is not affecting the results. This was 
confirmed by experiments in which the cadmium 
vapor pressure was varied over a wide range. 
In these experiments the results were shown to 
be independent of the temperature of the cad- 
mium supply over the range 170° to 240°C, 
corresponding to cadmium vapor pressures from 
4X10- to 2X10-* mm. 

The final results were obtained with the reso- 
nance cells at a temperature of 213°C and the 
cadmium supply at 180°C, corresponding to a 
cadmium vapor pressure of 9X10-> mm. It was 
possible to work with this very low vapor pres- 
sure of cadmium on account of the very high 
intensity of the source. 

The results obtained are shown in the form 
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? 
of plots of 1/Q against pressure in Figs. 1 to 6. 
In some cases the origin has been shifted verti- 
cally to avoid overlapping of the curves. The 
slopes of the lines in these diagrams and the 
quenching cross sections calculated from Eq. (1) 
are given in Table I. 


DISCUSSION 


The only previous work with which the above 
results can be compared is that of Lipson and 
Mitchell on hydrogen, deuterium and ammonia. 
Their results are compared with ours in Table IT. 
The agreement in the case of ammonia is as good 
as could be expected. The heavy line in Fig. 2 
represents our results, the dotted line those of 
Lipson and Mitchell. There is, however, a very 
large discrepancy in the case of hydrogen and 
deuterium. In an attempt to discover the cause 
of this a large number of experiments were 
made with hydrogen under various conditions. 
As mentioned previously, the results were inde- 
pendent of cadmium pressure, the length of the 
light path in the vapor and of the intensity. We 
therefore feel confident that our results are not 
affected by errors due to imprisonment of radi- 
ation. 








‘| At 
| bp - 























° 0.15 0.20 


S 010 
PRESSURE - MM 


Fic. 5. Quenching of 43261 by benzene and acetylene. 


In some experiments a spectrograph was sub- 
stituted for the camera and long exposures were 
made. These showed no detectable radiation 
other than 3261, and thus prove that errors 
were not introduced by scattered visible light, 
or by effects: concerned with other transitions 
involving the 5°P, level, i.e. by \\3404, 3466, 
and 3613. 

A number of runs were made with different 
samples of hydrogen. Tank hydrogen gave results 
only slightly different from those with the puri- 
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fied gas. Further, the addition of water vapor or 
oxygen in quantities much larger than the 
possible amounts present as an impurity caused 
a change in the quenching far too small to ac- 
count for the discrepancy. It appears impossible 
to account for the difference in the results on the 
basis of an impurity. The greatest ‘“‘quenching 
efficiency,” i.e. slope of the 1/Q—p curve, which 
we have found is 29.0 for butene-1. If we assumed 
that the true value for hydrogen is 2.45 as found 
by Lipson and Mitchell, and that our value of 
12.65 is due to the presence of butene in the 
hydrogen, it would be necessary to assume that 
40 percent of butene was present as an impurity. 
To account for the result on the basis of any 
reasonable amount of impurity, say even the 
extremely unlikely amount of 2 percent, would 
necessitate a substance with an efficiency 20 
times greater than butene-1, i.e. with a quenching 
cross section of about 70010-'* cm*. To say 
the least this is extremely unlikely, especially 
since such an impurity would have to possess a 
low enough molecular weight to be volatile at 
liquid-air temperatures. 

The only way in which our technique differed 
appreciably from that of Lipson and Mitchell 
was in the source of \3261. They used a cadmium- 
hydrogen discharge, while we used a cadmium- 
neon source. This difference might possibly fur- 
nish an explanation of the discrepancy. The 
cadmium-hydrogen discharge contains, in addi- 
tion to 43261, cadmium hydride bands of con- 
siderable intensity. Since the quenching process 
presumably involves the formation of CdH,!° it 
is possible to have resonance emission of the CdH 
bands from the quenching cell. If this were the 
case in Lipson and Mitchell’s investigation, some 
of the excited atoms quenched by forming CdH 
would subsequently emit CdH bands, and the 
result would be an apparent quenching efficiency 
much less than the true value. The values ob- 
tained for hydrogen, and for deuterium since the 
CdH and CdD bands would largely overlap, 


TABLE II. Quenching cross section, cm? X 10"®, 











Gas Lipson and Mitchell Present paper 
He 0.67 3.54 
2 0.19 1.80 
NH; 0.041 0.052 








10L. O. Olsen, J. Chem. Phys. 6, 307 (1938). 
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would be low, while those for substances which 
do not form CdH in the quenching process, such 
as ammonia, would not be affected. However, it 
seemed unlikely that such a large discrepancy 
could arise from this cause, and after corre- 
spondence with Professor Mitchell it was de- 
cided to reinvestigate the quenching with a lamp 
of the type used by Lipson and Mitchell. 

A lamp of the general type described by 
Ellett!! was therefore constructed. This consisted 
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Fic. 6. Quenching of 43261 by cyclopropane. 


of a quartz tube connected through graded 
seals to Pyrex electrode chambers containing 
large cylindrical aluminium electrodes. Hydrogen 
flowed into the lamp through a saturator con- 
taining a liquid cadmium-tin alloy, heated by an 
electric furnace, passed through the lamp, and 
was then pumped away by a diffusion pump. 
The lamp was operated at 5000 volts and about 
0.5 amp. This source gave \3261 relatively 
weakly, as compared with a cadmium-neon 
lamp, but gave the CdH bands, especially those 
at 44310 and 4500 with good intensity. As 
before, the radiation from the lamp passed 
through a Corex D filter. 

The results obtained were not as reproducible 
as those with the previous source, since fluctua- 
tions in lamp intensity were considerably greater 
than before. However, they showed conclusively 
that no appreciable difference was caused by the 
change in the source. In other words the quench- 
ing values obtained were in complete agreement 
with our previous values with the cadmium-neon 
source. 


u A. Ellett, J. Opt. Soc. Am. 10, 427 (1924). 











170 S. Ws &. 

As a further check, experiments were made 
with the cadmium-hydrogen lamp in which a 
chlorine filter was placed in the fluorescent beam 
emerging from the quenching cells. The filter 
was such that A3261 would be completely ab- 
sorbed, but AA4310 and 4500 would be fairly 
largely transmitted. Any fluorescent emission of 
these bands could thus be easily detected. Long 
exposures were made, but only negligible amounts 
of CdH fluorescence were observed. It may thus 
be concluded that no large error could have been 
introduced into Lipson and Mitchell’s experi- 
ments on account of fluorescence by cadmium 
hydride. 

The discrepancy between the two investiga- 
tions in the case of hydrogen and deuterium 
thus remains unexplained. We have, however, 
considerable confidence in our results since (a) all 
pertinent conditions were varied over wide 
ranges, (b) the cadmium vapor pressures em- 
ployed were so low that imprisonment effects 
cannot have been of importance, (c) a much 
lower quenching cross section for hydrogen and 
cadmium as compared with hydrogen and mer- 
cury would be rather surprising since the energy 
balance for 


Cd(5°P;) +H2=CdH+H 


is almost exact, and (d) there is some theoretical 
support for the higher value.” 


Accuracy of Results 


The percentage accuracy of the results for sub- 
stances which quench strongly is of the order of 
5 percent, and is much greater than that for 
substances which quench only slightly. Further, 
in the case of substances which quench very 
slightly the requirements as far as purity is 
concerned become impossibly exacting. As a re- 
sult the values given for the paraffins in Table I 
cannot be considered to be more than an upper 
limit to the quenching cross sections. Certainly 
the relative order of the quenching for the 
paraffins cannot be relied on. Actually it is prob- 
ably easier to remove propylene from propane 
than any other olefin from the corresponding 
paraffin, since ethylene is less reactive than the 
other olefins and the higher boiling butenes are 
more difficult to handle. The low value obtained 


2K. J. Laidler, J. Chem. Phys. 10, 34, 43 (1942). 
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TABLE III. Comparison of results with cadmium, 
mercury, and sodium. 








Quenching cross section, cm? X10'6 





Gas Cd(5°P1) Hg(6*P1) Na(32P)4 

He 3.54 8.60° 7.4 

CH, 0.0122 0.085% 0.11 

CoHe 0.024 0.59% 0.17 
3H» 0.012 2.32% 0.2 
butane 0.46, 0.64 5.88° 0.3 
CoH, 24.9 48¢ 44.0 
C;He 29.1 — 52.0 
C,Hs-1 35.2 —- 58.0 
C,Hs-2 30.6 — 58.0 
CeHe 28.4 59.9% 75.0 








@ See reference 7 in text. 

6 See reference 8 in text. 

¢ E. W. R. Steacie, Can. J.. Research B18, 44 (1940). 

4R. G. W. Norrish and W. MacF. Smith, Proc. Roy. Soc. A176, 295 
(1940). ‘ 
for propane probably means merely that it was 


the purest of the paraffins used. 


Comparison with Mercury and Sodium 


The values of the quenching cross sections for 
5’P; cadmium atoms obtained here are compared 
in Table III, wherever the data exist, with corre- 
sponding values for Hg(6°P1) and Na(?P) atoms. 

In a general way it will be seen that a paral- 
lelism exists among the results for the three 
substances. Thus hydrogen has a rather large 
quenching cross section in all three cases, the 
olefins quench very strongly, and the paraffins 
are relatively inefficient. In the case of the 
paraffins, however, there is a definite difference 
between the results for mercury and those for 
the other two metals. Thus with sodium and 
cadmium all paraffins quench very poorly. With 
mercury the lower paraffins quench equally 
poorly, but there is a steady rise in cross section 
with increasing molecular weight. Actually, no 
special precautions were taken to free the par- 
affins from olefin impurities in the investigation 
with mercury. Since propane and butane would 
be much harder to free from olefins by a simple 
bulb-to-bulb distillation, the results for these 
gases may well be too high. However, in view of 
the fact that Bates found a cross section for 
heptane much too large to be accounted for on 
the basis of an impurity, it seems certain that 
there is a real and large increase in the quenching 
cross sections with increasing molecular weight. 

The mechanism of the quenching of sodium, 
mercury, and cadmium has. recently been dis- 


1% J. R. Bates, J. Am. Chem. Soc. 52, 3825 (1930) ; 54, 


569 (1932). 
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cussed theoretically by Laidler” from the stand- 
point of potential energy surfaces. It appears 
certain in the case of hydrogen that quenching 
of sodium is of a physical nature, while that of 
cadmium involves the formation of CdH. In the 
case of mercury there is spectroscopic evidence 
against the formation of HgH in the quenching 
process.!° However, Laidler suggests that in 
reality the quenching act involves the formation 
of HgH, but that this possesses sufficient excess 
energy to dissociate and does so immediately. 

In the case of the paraffins a similar mechanism 
appears to hold. It is certain from photochemical 
evidence? that the quenching of Cd(5*P1) by 
paraffins involves the process 


Cd(5°P;)+RH=CdH+R, 


since the evidence is overwhelmingly in favor of 
a C—H bond split in the primary process, and 
the excitation energy of the 5*P; state is insuffi- 
cient to permit the reaction 


Cd(5°P1) +RH =Cd(5'So)) +R+H. 


It seems likely that quenching of Hg(6*P;) pro- 
ceeds by the same mechanism. From this point 
of view it is difficult to see why the cross sections 
of the paraffins increase with increasing molecular 
weight for mercury but not for cadmium. 

The very large quenching cross sections of all 
the olefins must mean that the quenching phe- 
nomenon is associated with the double bond. 
This leads to an interesting photochemical re- 
sult. It has been found that the cadmium photo- 
sensitized reactions of ethane and propane 
proceed readily, the primary step involving the 
formation of CdH and an alkyl radical. The 
strength of the C—H bond in the paraffins is 
certainly higher than the 87 kcal. excitation 
energy of the cadmium atom, so that the above 
reaction is possible only because of the 15 kcal. 
heat of formation of CdH. In the case of ethylene 
the quenching apparently involves the double 
bond, and practically no photosensitized reaction 
occurs. With propylene it might be expected that 
reaction would occur by the mechanism 
Cd(5*P,) +CH;CH = CH: 

=CdH+CH.CH =CH, 
etc. 


4D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 
10, 683 (1942). 
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analogous to the reaction with propane. Actually, 
virtually no photosensitized reaction occurred. 
Apparently quenching occurs again at the double 
bond by a mere transfer of energy. The 87 kcal. 
thus received is then incapable of causing reac- 
tion since even if it is transferred internally to 
a C—H bond it is insufficient to rupture it. 
The photosensitization experiments thus provide 
a very clean-cut example of the difference be- 
tween the two quenching mechanisms. 


Quenching of \2288 


For comparison with the above results and for 
photochemical reasons it would be of consider- 
able interest to have data on the quenching of 
the lower cadmium resonance line \2288. This is, 
however, very difficult to investigate. The 5'P, 
state has a life of only 2X10-® sec., i.e. about 
1200 times less than that of the 5*P; state. 
Hence, other things being equal, the actual 
quenching will be 1200 times smaller for a given 
value of og’. Appreciable quenching will thus 
only be obtained at high foreign gas pressures 
where pressure broadening will introduce large 
errors, or at relatively high cadmium vapor 
pressures where imprisonment will come into 
play. 

In order to avoid imprisonment effects, it is 
necessary to work at low cadmium vapor pres- 
sures since the absorption coefficient of the line 
is very large. A few experiments were made with 
hydrogen and with ethylene at cadmium vapor 
pressures corresponding to temperatures from 
245°C to 170°C. The apparent quenching fell off 
with the cadmium vapor pressure throughout the 
whole range. All that may be concluded from 
the results is that the quenching cross sections 
for ethylene and hydrogen are not greater than 
those with \3261, and are probably somewhat 
less. 

In order to obtain relative values which might 
be of some use, at least for photochemical pur- 
poses, a few measurements were made under 
comparable conditions for C2.H,, C3Hs, and He 
at 245°C. If it is assumed that imprisonment 
effects will be identical in all three cases, the 
quenching cross sections for the three gases are 
roughly in the ratio 


CH, ° H. : C;H;=1 tee F 
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The transference numbers for aqueous solutions of potassium bromide at 25°C have been 
determined by the moving boundary method. Cells with cadmium bromide and lithium bromide 
as cation indicators and with potassium iodate as anion indicator were employed. The measure- 
ments are more difficult than with potassium and sodium chlorides; with sheared junctions, a 
thorough investigation of the effect of indicator concentration is necessary for each solution 
studied, and the same is true with both sheared and autogenic boundaries of the effect of current 
density. The results differ slightly from those previously reported by Longsworth, but are 
consistent with the known value of ¢,°. Potassium bromide also differs from the alkali chlorides 
and acetates in that the Longsworth function ¢,” is not linear in the concentration, thus pre- 
cluding an independent extrapolation of the transference data. 





HIS research was originally undertaken to 

provide for potassium bromide transfer- 
ence numbers covering a range of temperature 
analogous to those reported from this laboratory 
for potassium chloride! and sodium chloride.’ 
In this way it was hoped to obtain from the 
recently reported conductances’ for this salt an 
additional check on the limiting ionic mobilities.‘ 
The work has, however, been much interrupted, 
and moreover it has become evident that unlike 
the alkali chlorides and acetates,> potassium 
bromide does not permit an unambiguous extra- 
polation of the transference data per se. We 
therefore report briefly our results for 25° since 
these differ somewhat from the existing trans- 
ference numbers for this salt reported in the 
literature. 

The experimental technique!? and the method 
of preparation of the potassium bromide solu- 
tions* have already been adequately described. 
Autogenic cation boundaries with cadmium 
bromide as indicator (Cell II), sheared anion 
boundaries with potassium iodate as indicator 
(III), and falling cation boundaries with lithium 
bromide as indicator (IV) were employed. It 
was found that while autogenic boundaries were 


1R. W. Allgood, D. J. LeRoy, and A. R. Gordon, J. 
Chem. Phys. 8, 418 (1940). 

2R. W. Allgood and A. R. Gordon, J. Chem. Phys. 10, 
124 (1942). 

3H. E. Gunning and A. R. Gordon, J. Chem. Phys. 11, 
18 (1943). 

4H. E. Gunning and A. R. Gordon, J. Chem. Phys. 10, 
126 (1942). 

5D. J. LeRoy and A. R. Gordon, J. Chem. Phys. 6, 398 
(1938), 7, 314 (1939). 


satisfactory up to 0.07 N, at 0.10 N it was im- 
possible to obtain reproducible results; thus, 
while at the lower concentrations the values of 
t, for a run showed no progression as the bound- 
ary moved up the tube, and were within ex- 
perimental error independent of the current 
used,® at tenth normal this was no longer the 
case. This behavior was apparently associated 
with the appearance of a heavy syrupy liquid 
layer above the cadmium anode, and no current 
density which still gave a visible boundary was 
sufficiently low to prevent its appearance. It 
was for this reason that falling LiBr boundaries 
were used for this solution. As is well known, 
such falling boundaries require careful adjust- 
ment of indicator concentration within narrow 
limits, but if this condition obtains, the resulting 
transference numbers are highly reproducible.’ 
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6 For example, the first and third entries at 0.07 N in 
Table I were run at 1.6 ma, while the second and fourth 
were run at 1.2 and 0.4 ma, respectively. 

7 The lithium bromide concentration used to obtain the 
entries in Table I was 0.069 normal; the first two entries 
were run at 2.3 ma and the third at 1.6. At lower indicator 
concentrations, e.g., at 0.065 N, the measured t, is low 
and tends to increase as the boundary moves down the 


172 





| 
1 








TRANSFERENCE OF KBr 173 


The measurements on anion boundaries with 
potassium iodate as indicator also lead to a 
rather surprising result. With potassium chloride 
solutions, MacInnes and Longsworth® found 
(and measurements in this laboratory! confirmed 
them) that the measured ¢_ was independent of 
indicator concentration within wide limits. 
In contrast, Fig. 1 shows the effect of indicator 
concentration for KBr/KIO; boundaries at 
0.05 N; in all the experiments, the current was 
the same (1.6 ma) and it is apparent that only at 


TABLE I. 








ty Sol. + Vol. 
Cell Cc obs. Corr. by ty 


II 0.014980 0.4835 +0.0004 0.4839 





0.4839 
II 015057 4834 -+0.0004 4838 
II 019892 4835 +0.0003 4838 

0.4837 
II .020169 4833  +0.0003 4836 
Ill .029026 5160 +0.0001 5161 
Ill 029933 5164 +0.0001 5165 

0.4837 


II 029981 4834 
II .030081 4835 
III 049864 5161 
Ill .049970 5163 
III 049970 5162 


+0.0003 4837 
+0.0003 4838 
—0.0001 .5160 
—0.0001 5162 
— 0.0001 5161 


II .050168 4834 +0.0003 4837 0.4839 
II 050120 4837 +0.0003 .4840 
II -050120 4836 +0.0003 4939 
II 050436 4834 +0.0003 4837 
II 070293 4842 +0.0004 .4846 
II 069961 4842 +0.0004 .4846 
II 070136 4838 +0.0004 4842 — 
II .070249 4840 +0.0004 4844 


IV .099807 4849 
IV .099893 4849 
IV 099692 4849 


+0.0003 4852 
+0.0003 4852 0.4852 
+0.0003 4852 








tube; at higher indicator concentrations, there is no pro- 
gression in the apparent ¢, but its value is again low. For 
both higher and lower indicator concentrations, the trans- 
ference number depends on the current used. See also 
reference 13. 

8 D. J. MacInnes and L. G. Longsworth, Chem. Rev. 11, 
171 (1932). 
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Cx103=0.0405 N does t_ equal 1—¢, (indicated 
by the dotted line in the figure). Moreover, it is 
only for this indicator concentration that the 
observed #¢_ is independent of the current—a run 
at 0.95 ma giving t-=0.5162. The results for 
0.03 N were similar to those indicated in the fig- 
ure, while for 0.07 and 0.10 NV we were unable at 
any indicator concentration to obtain anion 
transference numbers independent of the current. 

It would therefore seem that with potassium 
bromide absence of progression in the measured 
transference number as the boundary moves is 
not of itself conclusive evidence as to the reli- 
ability of the result. We therefore suggest, as a 
routine precaution in future investigations of 
other electrolytes, that the effect of changing 
the current and (with sheared junctions) of 
changing the indicator concentration should be 
studied for each solution, and that where pos- 
sible both ¢, and ¢t_ should be measured. 

The results are given in Table I as a func- 
tion of the concentration in moles per liter; 
the solvent and volume corrections® are listed 
in the fourth column. The latter® is given by 
— CAV/1000 where 


Cell II AV= —3Vca+}Vcapr, —t,: 0 KBr 
Cell III AV=+3Vca— 3 Vcapr, +t VkBr 
Cell IV AV= Vag— VaAgBr+t,-V KBr 


since in both cells II and III the cadmium anode 
forms the closed side of the cell, while in Cell IV 
the silver-silver bromide cathode is the closed 


®The volume correction for Cell IIIa of refreence 1 
(p. 420) is incorrectly stated; AV is given by Vag— Vage 
+t,:Vxc1. This has the effect of making the volume 
corrections for the three strongest solutions at 25° which 
employed this cell +0.0001, so that the resulting ¢_ in 
the sixth column of Table I are 0.5099, 0.5101, and 0.5098 
instead of 0.5098, 0.5100, and 0.5097 as printed. 
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TABLE II. TABLE III. 
ic 0.0005 0.001 0.002 0.005 0.010 "8g 0 0.0005 0.001 0.002 0.005 0.010 
ts 0.4846 0.4845 0.4844 0.4843 0.4841 (At)xpr (73.52) 72.50 72.08 71.52 70.46 69.35 
C 0.02 0.03 0.05 0.07 0.10 (At)xe1 73.52 72.48 72.06 71.48 70.40 69.27 
ty 0.4838 0.4837 0.4838 0.4843 0.4852 








side. To a close enough approximation, the 
volume corrections are given by +0.004,C, 
—0.004,C, and +0.002,C for cells II, III, and 
IV, respectively. The results are shown as a 
function of the square root of the concentration 
on the lower curve of Fig. 2, which also gives the 
previously reported data of Longsworth.’® It 
is at once apparent that the difference between 
his results and ours is somewhat greater than 
the apparent precision of the experiments. 
This is disturbing in view of the excellent agree- 
ment between the results obtained in this labora- 
tory for potassium and sodium chlorides and his 
25° data for the two salts, and correspondence 
with Dr. Longsworth, to whom we wish to 
express our thanks, has failed to explain the 
discrepancy. 

There is however some independent evidence. 
The limiting conductance! for potassium bro- 
mide, obtained by two entirely independent 
methods,*” is 151.67; from this and the known 
limiting mobility of potassium ion‘ viz., 73.52, 
t,° is 0.4847, to which our transference data have 
been extrapolated. A glance at the figure will 
show the difficulty of extrapolating Longsworth’s 
numbers to this value whereas ours for the lower 
concentrations actually lie on the theoretical 
limiting slope,!° indicated by the dotted line in 
the figure; while we regard this as fortuitous, 
it nevertheless permits a reasonable extrapola- 
tion. Secondly, our result for 0.1 N is in good 
agreement with the earlier value reported by 
MacInnes and Smith" for this concentration 


10 L. G. Longsworth, J. Am. Chem. Soc. 57, 1185 (1935). 

11 There is an obvious misprint in Table II of reference 3; 
ped first three lines of the table should be labeled 3, o 
and Ao. 

12 Grinnell Jones and C. F. Bickford, J. Am. Chem. Soc. 
56, 602 (1934). 

13 —D. A. MaclInnes and E. R. Smith, J. Am. Chem. Soc. 
45, 2246 (1923). 








viz., 0.485. Finally, it might be added that the 
transference numbers reported here, when used 
in conjunction with e.m.f. measurements on cells 
with transference, lead to values of the activity 
and osmotic coefficients for potassium bromide 
solutions which by isopiestic comparison are in 
excellent agreement with the now well-estab- 
lished thermodynamic data for potassium and 
sodium chloride solutions." 

The upper curve in Fig. 2 shows the Longs- 
worth function? ¢,"%, computed from our data, 
as a function of the concentration (indicated on 
the upper scale of the figure). This method of 
plotting has proven extremely useful with potas- 
sium and sodium chlorides and the corresponding 
acetates, but here the plot is definitely curved 
and thus does not give a satisfactory extrapola- 
tion to infinite dilution. 

Table II gives ¢, for round values of C, cor- 
responding to the continuous curves of Fig. 2, 
while Table III gives (A*)Kkpr, obtained from 
Table II and the conductances.’ The table also 
gives for comparison (A*)Kc1 obtained from 
Tables V and VI of reference 4. It is apparent 
that with increasing concentration there are 
deviations from the rule of independent ionic 
mobilities comparable with those found for 
chloride ion with potassium and sodium as 
cations (see Table VI of reference 4). There is 
however one difference; the anion chloride had 
the greater mobility with sodium, the slower 
cation, while here the cation potassium has the 
greater mobility with bromide, the faster anion. 
Whether this is of general significance or not, 
can only be settled by further experiment. 


4 A. S. Brown and D. A. MacInnes, J. Am. Chem. Soc. 
57, 1356 (1935); T. Shedlovsky and D. A. MaclInnes, 
J. Am. Chem. Soc. 59, 503 (1937); W. J. Hornibrook, 
G. J. Janz, and A. R. Gordon, J. Am. Chem. Soc. 64, 513 
(1942); G. J. Janz and A. R. Gordon, J. Am. Chem. Soc. 
65, 218 (1943). 
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The theory of dielectric polarization of polar liquids is used to calculate the dielectric con- 
stants of water and some alcohols. Quasi-rigid coordination models, based upon x-ray structural 
investigations, are employed in the calculation of a function g which measures the orientational 
correlation between neighbors in the liquid. It is pointed out that the correlation function g is 
closely related to the thermodynamic “abnormality” of polar liquids. 





INTRODUCTION 


N an earlier article a general theory of the 
dielectric polarization of polar liquids was 
formulated.! The role of hindered molecular rota- 
tion was clarified and it was shown how corre- 
lation in the orientation of neighboring molecules 
might augment or diminish the contributions of 
the individual molecular dipoles to the dielectric 
constant of a liquid. 

In the present article, we shall discuss some 
applications of the theory to water, the aliphatic 
alcohols, and certain other hydrogen bonded 
liquids. On the basis of simplified but reasonable 
models of the local structure maintained by each 
molecule in its environment by hydrogen bonding 
with its neighbors, we are able to calculate the 
dielectric constants of these liquids with some 
degree of success. Due to strong orientational 
correlation, the effect of the individual molecular 
dipoles is greatly enhanced. This circumstance 
makes it possible for water and the alcohols to 
have dielectric constants considerably greater 
than those of certain other polar liquids the 
molecules of which possess dipole moments of 
comparable magnitude. 

According to our theory of dielectric polariza- 
tion, the dielectric constant ¢ of a liquid is re- 
lated to the polarization P» by the equation 


(e—1)/3=[3e/(2e+1) ]Po/v (1) 
Po= (4nN/3)La+(y-B/3kT)], 
where v is the molal volume, a the optical 


polarizability, N Avogadro’s number, and kT 
the product of Boltzmann’s constant and the 


1J. G. Kirkwood, J. Chem. Phys. 7, 911 (1939); L. 
Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 


absolute temperature. In the scalar product 
up, wu is the molecular dipole moment and @ is 
the sum of the molecular dipole moment and the 
average electric moment induced by the molecule 
in its environment by hindering the rotation of 
its neighbors relative to itself. When Po/v is 
large compared with unity, Eq. (1) may be 
simplified to 


e—1=(9/2)(Po/v) (2) 


and the dielectric constant becomes proportional 
to the density of the liquid. 


WATER 


If the orienting effect of a water molecule is 
assumed to extend only to its first shell of 
neighbors and if the region of exclusion of the 
first coordination shell is roughly spherical, the 
product u-@ may be expressed in the form 


yu=gn’, 

g=1+2(cos y)w, (3) 
where z is the average number of neighbors in the 
first coordination shell. The coordination number 
z may be estimated from the area under the first 
peak of the radial distribution function of the 
liquid. These approximations will be employed 
in the calculation of the dielectric constant of 
water. 

We shall assume water to be approximately 
tetrahedrally coordinated in accordance with the 
Bernal-Fowler model.? A directed hydrogen bond 
about which there is free rotation but no bending 
is assumed to exist between neighboring mole- 
cules. The bond is assumed to be directed along 


2 J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
(1933). , 
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Fic. 1. Hydrogen bonding in water. Free rotation is 
assumed about the hydrogen bond between adjacent water 
molecules. The short dashes indicate regions of negative 
one and are in a plane perpendicular to the H—O—H 
plane. 


the O—H bond of the one molecule toward the 
oxygen atom of the second molecule passing 
through the concentration of negative charge due 
to one of the non-bonding 2 electrons of the 
latter, as exhibited in Fig. 1. With this quasi- 
rigid model of the local structure, we calculate 


7] 
(cos y) a = cos? > g=1+2 cos? (4) 


H 
where @ is the 0c bond angle in the water 
H 


molecule, which we take to be 105°, the best 
spectroscopic value. The dipole moment yz of a 
water molecule in the liquid differs from its 
moment po in the vapor due to induction of an 
additional moment by the dipole field of its 
neighbors. If @ is the polarizability of the water 
molecule, we have approximately 


rnffi-(omee'fe)) 0 


where a is the distance between neighbors. 

We are now able to estimate the dielectric 
constant of water with the aid of Eqs. (2), (4), 
and (5). The coordination number z and the inter- 
molecular distance a are available from the x-ray 
scattering data of Morgan and Warren.’ In the 
temperature range, 0°C to 83°C, the average 
coordination numbers are non-integral, lying 
between 4 and 5. Although non-integral coordi- 
nation numbers exceeding 4 are not entirely 
consistent with our simplified model of the 
neighbor structure, we shall nevertheless employ 
them in the calculation. This is equivalent to 
approximating (cos y)4 by cos? 6/2 for small 


8 J. Morgan and B. E. Warren, J. Chem. Phys. 6, 666 
(1938). 


G. OSTER AND J. 








G. KIRKWOOD 


deviations from the tetrahedrally coordinated 
structure. 

In Table I the theoretical values of the di- 
electric constant of liquid water are compared 
with the experimental values of Wyman.‘ The 
agreement between the theoretical and experi- 
mental values of the dielectric constant is re- 
markably close at 25°C. At 83°C, the values 
differ by about 13 percent. In view of the obvious 
limitations of the model of local structure upon 
which the theoretical values are based, the agree- 
ment may be considered satisfactory. 

We remark that ‘the experimental values of g 
decrease with increasing temperature whereas 
the theoretical values increase. This trend is in 
accord with failures of the simple model which 
might be expected. A partial destruction of the 
tetrahedrally coordinated structure would natu- 
rally tend to decrease the correlation in orienta- 
tion of neighbors and thus diminish (cos 7)w. 
Bending of the hydrogen bond would have the 
same effect and would be expected to be more 
pronounced the-higher the temperature. 


ALIPHATIC ALCOHOLS 


The structure of the liquid alcohols has been 
extensively investigated by the x-ray scattering 
technique.® These investigations -have led to the 
conclusion that the average coordination number 
of a molecule is 2. It has therefore been suggested 


TABLE I. The dielectric constant of liquid water. 











t(°C) v z a(A)_ g(obs.) g(calc.) e(obs.) e(calc.) 
O 18.02 44 2.90 2.75 2.63 88.0 84.2 
7 1607 45 292 266 241 YO5 782 
62 18.34 4.9 3.00 2.56 2.82 66.1 1235 
83 18.58 4.9 3.05 2.49 2.82 59.9 67.5 








6 =105°; a =1.5A3 (constant over temperature range) ; wo = 1.88 X107!8 
e.s.u. 


that the molecules are arranged in chains by 
means of hydrogen bonds. 

Cole® has discussed the application of the 
present theory to the liquid alcohols, which we 
wish to elaborate here. We shall assume that each 


4J. Wyman, Jr., Phys. Rev. 35,623 (1930). 

5W. H. Zachariasen, J: Chem. Phys. 3, 158 (1935); 
G. G. Harvey, J. Chem. Phys. 7, 878 (1939); Pierce and 
MacMellan, J. Am. Chem. Soc. 60, 779 (1938). See also 
L. Pauling, Nature of the Chemical Bond (Cornell Uni- 
versity Press, Ithaca, 1940), second edition, p. 305. 
6 R. H. Cole, J. Chem. Phys. 9, 251 (1941). 
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alcohol molecule is a member of an infinite chain 
and is correlated in orientation only with mem- 
bers of its own chain. We shall further suppose 
that the average configuration of a specified 
chain, due to curling up by internal rotation, 
is approximately spherical. Directed hydrogen 
bonds of the type existing in water are assumed 
to link each alcohol molecule with its two nearest 
neighbors in the chain, as exhibited in Fig. 2. 
Free rotation around the hydrogen bonds is 
allowed, but bending is excluded. 

For chainwise coordination, the factor g be- 
comes 


g=142 5 (cos Ya)m (6) 


n=1 


where 7, is the angle between the dipole moments 
of a molecule and a chain neighbor separated 
from it by bonds. On the basis of our assump- 
tions concerning the hydrogen bond, we may 
calculate’ 


6 n 
(COS Yn) av =1( cos? -) , 


(7) 
_ (ua+ur) (at ur Cos 8) 


Cun?-+ur?+2urun cos 0) 





where wy and ye are the components of the dipole 


moment of the alcohol molecule ROH along the - 


TABLE II. Dielectric constants of the normal 
alcohols at 20°C. 











v n g(obs.) g(calc.) e(obs.) €(calc.) 
Methyl alcohol 40.4 1.33 2.94 2.57 32.8 29.2 
Ethyl alcohol 58.3 1.36 3.04 2.57 24.6 21.3 


n-Propyl alcohol 75.5 1.38 3.07 2.57 19.5 17.3 
n-Butyl alcohol 91.5 1.40 3.21 2.57 18.0 14.6 
n-Amylalcohol 108.0 1.41 3.43 2.57 15.8 12.8 








0 =105°; pw =1.70 10-18; pp =0.7 10-18; zy =1.6 “10-18, 

O—H and O—R bonds, respectively, and @ is 
R 

the KC, bond angle, taken to be 105°. For the 


infinite chain we obtain by summation of the 
series resulting from the substitution of Eq. (7) 
into Eq. (6) 

g=1+2f cot? (6/2). (8) 


asa and Kirkwood, J. Am. Chem. Soc. 63, 390 


HINDERED ROTATION 





IN DIELECTRICS 


In order to estimate f, we use the bond dipole 
moments which have been assigned to methyl 
alcohol in the vapor phase.® In order to estimate 
the ratio of liquid dipole moment y» to the vapor 
dipole moment yo, we employ Onsager’s® sphere 
formula 

2e+1 n?+2 


| ae HO, (9) 
2e+n? 3 








where 1 is the optical refractive index. The use 


Fic. 2. Hydrogen bonding in alcohol. 


of Eq. (9) in the liquid alcohols is somewhat 
questionable, since the model on which it is 
based scarcely applies. However, it serves for a 
rough estimate of the induction effect from 
neighbors on the dipole moment of a given mole- 
cule, and a detailed calculation of the dipole 
field from chain neighbors is more complicated 
than the limitations of our idealized coordination 
model justifies. 

In Table II the theoretical values of the di- 
electric constants of several of the normal 
aliphatic alcohols are compared with the ob- 
served dielectric constants.'° The agreement with 
experiment is not as good as in the case of water. 
It is, however, remarkably good in view of the 
limitations of the assumed model of local struc- 
ture. The fact that the experimental value of g 
exceeds the theoretical value suggests that there 
is a residual correlation in the orientations of 
molecules in different chains. It will be observed 
that the experimental correlation parameter in- 
creases with increasing chain length in the case 
of the normal alcohols. 


8 Eucken and Meyer, Physik. Zeits. 30, 397 (1929). 

9L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 

10 Dielectric constants and dipole moments values mainly 
from P. Debye and Sack, Tables Annuelles de Constantes et 
Donnees Numeriques, Volume XI, Booklet 2 (1937); 
Volume XII, Booklet 32 (1939). Refractive indices and 
densities from International Critical Tables. 
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TABLE III. Correlation parameter of ethyl alcohol 
as a function of temperature. 











4(°C) g(observed) 

—90 3.33 

— 30 3.24 
20 3.03 
70 2.84 








TABLE IV. Influence of organic radical on the correlation 
parameter of an alcohol ROH at 20°C. 











g(observed) 
Secondary butyl alcohol 2.83 
Tertiary butyl alcohol 2.38 
Isobutyl alcohol 3.36 
Benzyl alcohol 2.08 








The correlation parameter g differs very little 
for the several alcohols. It is interesting to ex- 
amine the effect of temperature and the effect of 
the organic radical on this parameter. In Table 
III the experimental value of g for ethyl alcohol 
is presented as a function of temperature. As one 
would expect, g decreases with increasing tem- 
perature. In Table IV, the effect of the organic 
radical on the correlation parameter of the 
alcohols is exhibited. Here there does not appear 
to be any striking regularity in the manner in 
which the organic radical influences orientational 
correlation. It is, of course, reasonable that g 
should be low for tertiary butyl alcohol, since 
steric interference with hydrogen bond formation 
is likely. 


OTHER POLAR LIQUIDS 


The correlation parameter g is a useful index 
to the degree of hindered rotation in a liquid. 
Polar liquids possessing values of g differing 
appreciably from unity form a class which is 
frequently, although perhaps unfortunately, des- 
ignated as the class of associated liquids. Polar 
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TABLE V. Correlation parameters of several polar liquids. 








e(ob- 
served) g(observed) 





1(°C) po *1018 
Hydrogen cyanide 20 2.80 116 4.1 
Hydrogen fluoride 0 (1.8)* 83.6 3.1 
Hydrogen peroxide 0 2.13 91 2.8 
Ammonia 15 1.48 17.8 1.3 
Ethyl ether 20 1.35 4.4 1.7 
Acetone 20 2.85 21.5 BI 
Nitrobenzene 20 3.90 36.1 1.1 
Ethyl bromide 20 1.80 9.4 1.1 
Pyridine 20 2.20 12.5 0.9 
Benzonitrile 25 3.90 25.2 0.8 








* Estimate by L. Pauling, Nature of the Chemical Bond, p. 292. 


liquids with a value of g near unity fall into the 
class of normal liquids. As a basis for classifica- 
tion of liquids according to their thermodynamic 
behavior, it appears that the correlation param- 
eter g is more significant than the molecular 
dipole moment. 

Orientational correlation is determined not by 
the molecular dipole moment alone, but also by 
the location of the permanent charge distribution 
in the interior of a polar molecule. When the 
charge distribution is exposed, as in the molecules 
of hydrogen bonded liquids, neighbors can as- 
sume relative configurations in which there is 
very strong interaction. When the charge dis- 
tribution is shielded, as in nitrobenzene, such 
configurations are excluded by intermolecular 
repulsive forces, and orientational correlation is 
not great. 

In Table V, the correlation parameter is 
presented for a series of polar liquids. No at- 
tempt will be made at the present time to in- 
terpret g on the basis of specific structural 
models, although this offers interesting possi- 
bilities in several instances. For a more com- 
prehensive survey of the data for other liquids, 
the reader is referred to an article by J. N. 
Wilson." 


1 J. N. Wilson, Chem. Rev. 25, 377 (1939). 
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Diffraction of X-Rays by Sodium Stearate at Room Temperature* 


James W. McBain, ALEXANDER DE BRETTEVILLE, JR., AND SYDNEY Ross 
Department of Chemistry, Stanford University, Stanford University, California 


(Received June 22, 1942) 


The alpha-form of sodium stearate, originally reported by Thiessen and Stauff, and of which 
a sample was made available to us, has been re-examined and the data found to be in general 
conformity with those of the previous investigation. Re-examination of their monoclinic beta- 
form of sodium stearate, obtained by heating the rhombic form above 54°C, yields a different 
set of data. A new, more stable, monoclinic form of sodium stearate has been examined and 


designated as the gamma-form. 





VER since 1923 when the x-ray study of soap 
systems was begun at Bristol University,' 
it has been recognized that this approach gives 
invaluable information. Such work has become 
more urgent since in the last few years it has 
been found that every soap exists in a large 
number of forms, only some of which are typically 
crystalline or liquid crystalline. Even isotropic 
solutions in water show spacings corresponding 
to lamellar micelles.2 Nearly two dozen papers 
dealing with the x-ray examination of soap 
systems have now appeared, coming from many 
laboratories. 
As a preliminary to an extensive investigation 
of soap systems by means of x-ray diffraction, 
being carried out at Stanford University, it was 


necessary to review the material already pub- ° 


lished. Of this the paper by Thiessen and Stauff* 
on crystalline forms of sodium stearate and 
sodium palmitate at room temperature has gone 
furthest towards presenting a complete x-ray 
analysis of these substances. This valuable and 
important paper is marred by a number of 
numerical errors, so that -significant revision is 
here suggested. The present paper contains 


* This manuscript was first completed for publication on 
November 6, 1940; cf. reference 6. 

1S. H. Piper and E. N. Grindley, Proc. Phys. Soc. 
London 35, 269 (1923). 

2K. Hess and J. Gundermann, Ber. 70, 1800 (1937); 
Hess, Kiessig, and Philippoff, Naturwiss. 26, 184 (1938); 
Hess, Philippoff, and Kiessig, Kolloid Zeits. 88, 40 (1939) ; 
J. Stauff, Naturwiss. 27, 213 (1939); H. Kiessig and W. 
Philippoff, ibid. 594 (1939). 

’P. A. Thiessen and J. Stauff, Zeits. f. physik. Chemie 
A176, 397 (1936). 

* They themselves have published a subsequent note, 
Zeits. f. physik. Chemie A177, 398 (1936), correcting a very 
few of the numerical errors and admitting that their 
published analysis of the monoclinic form is no longer in 
complete agreement with all the measurements, including 
new and improved but unpublished data. 


179 


further data on the alpha- and two monoclinic 
forms of sodium stearate, beta- and gamma-. By 
courtesy of Dr. J. Stauff of Berlin the authors 
were able to obtain, through Dr. Vold, a sample 
of the alpha-sodium stearate actually used in the 
first investigation, so that there is no doubt but 
that the present paper deals with some of the 
very same specimen that was examined by 
Thiessen and Stauff. 


MATERIALS 


Two samples of alpha-sodium stearate were 
available, one (A) originating from Kahlbaum 
purified by Dr. J. Stauff in Berlin, and the other 
(B) prepared by Dr. R. D. Vold in this labora- 
tory from Eastman stearic acid No. 402. The 
alpha-crystals are obtained by slow crystalliza- 
tion from 1 percent solution in 95 percent ethyl 
alcohol and subsequent drying in air. The mono- 
clinic beta-form was obtained from the alpha- 
form in both cases by heating the latter above 
54°C. On cooling, no trace of the original alpha- 
form was to be found. 

The new y-monoclinic form was prepared by 
merely melting and cooling the soap. Stearate 
prepared from stearic acid Kahlbaum® and from 
Eastman stearic acid No. 402 was employed. 


METHOD 


The radiation was supplied by a copper-target, 
water-cooled x-ray tube made by Mr. J. F. Lee 
and purchased through a grant from the Cyrus 
M. Warren Fund of the American Academy of 
Arts and Sciences. The tube was operated at 
28 kv, with a current of 28 ma. The radiation 
was filtered through nickel foil to make it es- 


5 McBain, Vold, and Frick, J. Phys. Chem. 44, 1014 
(1940). 
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TABLE I. Comparison of x-ray diffraction lines for alpha-sodium stearate. 























































Results of Thiessen and Stauff, from Results of present investigation, from 
rotation photographs powder photographs 
sin? 0 sin? 0 sin? 0 d/n sin? 6 sin? 6 

hkl calc.* obs. Int. hkl calc. A (A) (B) Int. 

002 0.00089 0.00085 m 002 0.00087 26.1 0.00093 0.00089 m. 

003 0.00200 0.00199 s 003 0.00197 17.4 0.00187 0.00201 m. 

004 0.00355 0.00360 m. 004 0.00349 13.0 0.00344 0.00347 w. 

005 0.00555 0.00553 S. 005 0.00546 10.4 0.00567 0.00547 Ww. 

007 0.0109 0.0108 m. 

008 0.0142 0.0140 Ww 

009 0.0180 0.0184 w. 

110 0.0161 0.0159 V.w. 

150 0.1826 0.0198 V.w. 

020 0.0277 0.0278 v.s. 020 0.0279 4.60 0.0280 0.0279 V.s. 

a 0.0297 0.0299 w. 023 0.0299 4.55 0.0305 0.0305 V.W. 

200 0.0366 0.0365 v.s. 200 0.0368 4.01 

120 0.0369 0.0370 m. | 120 ost 3.99 0.0371 0.0369 v-8. 

203 0.0386 0.0386 m. 

123 0.0389 0.0386 V.W. 

205 0.0422 0.0422 w. 

125 0.0425 0.0420 V.w. 
210 0.0438 3.67 0.0440 0.0438 V.S. 

210 0.0436 0.0474 v.S. 207 0.0475 3.53 0.0468 V.w. 

207 0.0475 0.0481 V.W. on 7 

213 0.0456 0.0493 m. 215 .0492 3.4 . 

215 0.0491 0.0529 w. | 216 0.0816) 3.39 0.0508 0.0509 m. 

217 0.0545 0.0580 V.w. 

220 0.0644 0.0642 s. 220 0.0647 3.02 0.0643 0.0644 s. 

223 0.0664 0.0699 w. 

130 0.0716 0.0718 w. 130 0.0720 2.87 0.0723 0.0721 w. 

225 0.0699 0.0736 V.w. 

310 0.0894 0.0829 m. 229 aaa 2.67 0.0824 V.w. 
310 0.089 2.56 
= Ops08} ye 0.0912 0.0900 w. 
315 0.0952 2.49 0.0941 w. 

230 0.0991 0.0985 s 230 0.0995 2.44 0.0995 0.0997 Ss. 

233 0.1011 0.1010 w 

235 0.1046 0.1040 w. 

320 0.110 0.107 m. 

323 0.112 0.109 V.w 323 0.113 2.28 

237 0.1099 0.1102 V.w. 237 0.110 2.32 0.111 0.110 w. 

040 0.111 0.111 m. 040 0.112 2.30 

325 0.116 0.117 v.v.w. | 325 0.116 2.26 0.116 

239 0.1171 0.1175 v.w. | 239 0.117 2.25 ; ” 

043 0.1130 0.1213 w. 

140 0.1201 0.1215 V.w. 

410 0.1535 0.135 m. 145 0.126 2.16 0.125 0.125 V.w. 

240 0.1476 0.1480 m. 240 0.148 2.00 0.149 0.148 V.w. 
410 0.154 1.96 0.155 0.154 m. 

420 0.174 0.176 m. 420 0.175 1.84 0.177 0.175 w. 

430 0.209 0.202 w. 430 0.210 1.68 0.206 0.206 V.V.W. 
447 0.269 1.48 0.271 0.268 w. 
543 0.343 Loe 0.345 0.341 V.w. 
620 0.359 1.28 0.362 0.369 V.w. 
630 0.394 1.22 0.397 V.W. 


















sentially monochromatic. A circular camera of 
57.2-mm radius was used to obtain the short 
spacings and a flat plateholder at distances up 
to 250 mm to resolve the long spacings. The data 
were scrutinized to ensure that no tungsten or 
Cu K&B lines were included in the tables. 


5a A. de Bretteville, Jr., Rev. Sci. Inst. 13, 481 (1942). 


* Many of these values are corrected to their own equation as published by them, although their equation does not completely accord with the 
data for the unit cell from which it was deduced. It will be noted that the recalculated data in many cases no longer fit the observed data. 


ALPHA-FORM OF SODIUM STEARATE 


Both samples of the alpha-sodium stearate as 
observed by us in 1940° gave x-ray diffraction 
lines that agreed in position with each other and 


6 A. de Bretteville, Jr., Ph.D. thesis, Stanford University, 
1941; obtainable in bibliofilm from the U. S. Department 
of Agriculture. 
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_TABLE II. Comparison of x-ray diffraction lines for beta-monoclinic sodium stearate. 











Results of Thiessen and Stauff, from Results of present investigation, from 
rotation photographs powder photographs 
sin? 6 sin? @ sin? 6 d/n sin? 6 sin? 0 
hkl calc.* obs.** Int. hkl calc. A (A) (B) Int. 
002 0.00111 0.00112 m 002 0.00109 23.3 0.00103 0.00119 Ss. 
003 0.00251 0.00246 s 003 0.00245 15.5 0.00253 0.00238 V.S. 
004 0.00446 0.00383 m 004 0.00435 11.7 0.00462 0.00436 w. 
005 0.00697 0.00696 s 005 0.00679 9.32 0.00690 0.00705 m. 
006 0.0100 0.0103 w. 
007 0.0137 0.0130 m. | 007 0.0133 6.66 0.0131 0.0136 w. 
110 0.0333 0.0335 v.s. | 110 0.0321f 4.29 0.0331 0.0334 V.S. 
200 0.0496 0.0498 V.S 200 0.0393 3.88 0.0393 0.0392 V.S. 
203 0.0418 0.0520 w - 
204 0.0553 3.27 0.0540 V.W. 
210 0.0616 3.09 0.0600 0.0588 w. 
214 0.0775 2.76 0.0761 0.0774 w. 
020 rye 0.0851 s. 020 0.0891 257 0.0894 V.V.W. 
205 -0393 
ad 0.0990) 0.0975 vw. 
310 0.132 0.131 w. 310 0.111 2.31 0.112 0.112 w. 
400 0.198 . 0.192 s. _ 
405 0.193 1.75 0.189 0.189 V.W. 
220 0.133 0.212 m. 
330 0.289 1.43 0.273 0.289 V.w. 
333 0.304 1.39 0.309 V.w. 
335 0.317 1.36 0.317 0.313 V.W. 
040 0.334 0.346 m. 
023 0.0860 
024 oaen) 0.888 w. 











* These values are corrected by us according to their revised equation, 1/d? =0.0209h2 +-0.0352k? +-0.00047/? —0.0029h/, to agree with their unit 


cell dimensions. 
** As tabulated by Thiessen and Stauff. 


t In other photographs this value has been observed as 0.0324, and as low as 0.031. 


with the positions reported by Thiessen and 
Stauff.? Because of certain arithmetical errors in 
the published results it is deemed advisable to 
print a table of indices in which the sin® @ 
values calculated from their equation are cor- 
rected. The first three columns of Table I re- 
print the data of Thiessen and Stauff with these 
corrections. The last five columns report the 
data obtained by the present authors and a 
suggested table of indices, some of which differ 
from those assigned by Thiessen and Stauff. 
Although these changes are based solely on 
powder diagrams, they provide in general a better 
agreement between calculated and observed 
values, until a more advanced technique has 
been applied. The equation used to calculate the 
sin? @ values for the powder photographs is 
based on the mean of three experimental values, 
including that of Thiessen and Stauff, for the 
200 plane, the 020 plane and the 00/ planes. 
The equation is 


sin? @=0.009194h? +-0.006974k? +-0.0002184/?. 


The column headed (A) gives the experimentally 
determined values of sin? @ for the sample ob- 


tained from Dr. Stauff, and that headed (B) 
gives the values for the sample obtained from 
Dr. Vold. We have not suggested indices for lines 


" not obtained from these two powder photographs. 


The dimensions of the unit cell, as reported by 
Thiessen and Stauff, are a=8.04A, b=9.12A, 
and c=51.77A. Their equation yields a=8.04A, 
b=9.24A, and c=51.6A. The calculated density 
for exactly eight molecules in the unit cell is 
1.057. The density measured by them by im- 
mersion in xylene was 1.0361. The mean values 
used in obtaining our equation are a=8.03 
+0.03A, b=9.22+0.02A, and c=52.1+0.3A. 
The calculated density is then 1.052. 


THE BETA-MONOCLINIC FORM OF SODIUM 
STEARATE 

In -a supplementary note‘ to their original 
paper, Thiessen and Stauff revised the mono- 
clinic equation for sodium stearate, but still 
based it on their previously published dimensions 
of the unit cell. When either equation is used to 
test the suggested indices for the observed 
spacings they are found to be inconsistent with 
their tables, and, accordingly, a new investiga- 








182 McBAIN, DE BRETTEVILLE, AND ROSS 





Fic. 1. Long spacings of the y-form of monoclinic sodium 
stearate; flat plate 125 mm from sample. 


tion of their beta-monoclinic form was required. 
This was undertaken on the two samples that 
were available, one (A) obtained from the 
alpha-form prepared by Thiessen and Stauff 
and the other (B) from the alpha-form prepared 
by Dr. Vold. In both cases the beta-monoclinic 
form was prepared from the alpha by heating the 
latter somewhat above 54°C. 

The most significant revision of the data for 
the beta-monoclinic form concerns the directly 
observed position of one of the two very strong- 
est lines. This line, indexed by Thiessen and 
Stauff as 200, is recorded as sin® 6=0.0498, 
whereas we fail to observe any diffraction effect 
at that position. Instead we obtain a very strong 
line at sin? 2=0.0393, not reported by them. 
Those two lines are matched by virtue of their 
intensities, and also by measurement of their 
actual photograph (Fig. 12), and we conclude* 
that the real position of the 200 diffraction 
effect is at sin? @=0.0393. 

We have made use of Thiessen and Stauff’s 
oscillation photograph, Fig. 12, in assigning 
indices. On inspection Fig. 12 exhibits imperfect 
orientation in the sample. Not all spots seen 
on the equator have hO/ indices; because they 
can be traced by a continuous line, of constant 
sin 6, to more intense spots on the first Jayer. 


* The interesting suggestion has been made to us that 
possibly their Fig. 12 was obtained by Fe Ka radiation 
at a crystal to plate distance of 37.8 mm, and that, there- 
fore, their beta-form was not that here found. However, 
if so, the strongest line in Fig. 12 finds no place in their 
table, and the results would not agree with the values given 
by Thiessen and Stauff for the dimensions of the unit cell; 
nor would the results agree with our data. If taken as we 
= with Cu Ka radiation their Fig. 12 agrees with our 

ata. 





Fic. 2. Powder photograph of the y-form of monoclinic 
sodium stearate; flat plate 45 mm from sample. 


Conversely, some of the spots appearing on the 
first layer are similarly traced to more intense 
ones on the equator. 

Table II presents the data of Thiessen and 
Stauff for the beta-form of monoclinic sodium 
stearate, together with the result of the present 
investigation. 

Our equation for the beta-monoclinic form is 


sin? 6=0.009824h? +0.02227k? 
+0.0002717/ —0.001454A1. 


The monoclinic angle 6 can be derived from the 
ratio of the measured value of ¢ sin B to the 
theoretical length 52.4A for the double molecule 

TABLE III. Powder diffraction data for the new gamma- 


form of monoclinic sodium stearate, prepared by melting 
and allowing to cool. 








sin? @ sin? 0 





(Kahl- d/n (East- d/n 
baum) A Int. man) A Int. 
0.000333 42.2 v.S. (n=1) 
0.00119 22.3 m. 0.00125 21.8 (n=2) s. 
0.00274 14.7 s. 0.00266 14.9 (n=3) V.S. 
0.00436 11.7 (n=4) w. 
0.00765 8.80 m. 0.00728 9.02(n=5) m. 
0.0120 7.03(n=6)  v.w. 
0.0142 6.46 w. (n=7) 
0.0266 4.72 v.w. 0.0262 4.74 w. 
0.0288 4.54 Ss. 0.0286 4.54 Ss. 
0.0334 4.21 v.w. 0.0333 4.21 w. 
0.0371 4.00 v.v.s. 0.0366 4.02 é 
0.0392 3.88 s. 
0.0524 3.36 v.w. 0.0515 3.39 V.w. 
0.0594 3.16 v.w. 0.0577 3.20 V.w. 
0.0699 2.91 m. 0.0697 2.91 m. 
0.0732 2.84 m. 
0.102 2.41 v.v.w. 
0.112 208 Wa 292 2.30 w. 
0.118 2.24: vw. @121 2.21 w. 
0.187 1.78 v.w. 0.188 1.77 V.w. 
0.204 1.70 v.w. 0.202 1.71 V.w. 
0.275 1.47 v.w. 0.270 1.48 V.V.W. 
0.383 1.24 v.v.w. 
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of sodium stearate,’ or the experimentally ob- 
served value for the alpha-form, 52.1A, giving B 
equal to 63°. Thiessen and Stauff found 62.4°, by 
the method of Brill and Meyer. 

Hence we have the following values: a sin 8 
=7.77A; a=8.7A; b=5.1A; ¢ sin B=46.6A; 
B=63°; e=b/(asin 8B) =0.66; density =0.99 for 
four molecules in the unit cell. ¢ has an average 
value of 0.64 for a large number of paraffin chain 
compounds. According to molecular models, 
this is the approximate proportion of the cross 
section of a hydrocarbon chain; the zig-zag 
presumably lies in the a-c plane of the crystal. 
The spacing 6, taken here as 5.16+0.11A, appears 
slightly variable as we have found with other still 
unpublished photographs; likewise the layer lines 
in Thiessen and Stauff’s photograph are irregular. 


A NEW, MORE STABLE, MONOCLINIC GAMMA- 
FORM OF SODIUM STEARATE? 

Kahlbaum sodium stearate, when melted and 

allowed to cool (‘‘cast’’) was found to be dis- 


7L. Pauling, Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1940). 
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tinctly different from the beta-monoclinic form. 
The long spacing is 44.6A (Fig. 1), as compared 
with 46.6A for beta-monoclinic and 52.1A for 
alpha-sodium stearate. Likewise the side spacings 
are distinctive, as appears from Table III and 
Fig. 2. 

The new gamma-form resembles in long 
spacings the ‘‘b” form of stearic acid found by 
Francis, Piper, and Malkin,® ¢ sin B=44.0A, 
also observed by Muller as 43.8A, whereas the 
beta-monoclinic form of sodium stearate is 
like their “‘a’’ form of stearic acid 46.6A. The 
polymorph corresponding to their “‘c’’ form, 
39.5-40A, remains to be found, although we 
may have observed it with a sample of thallous 
palmitate given us by Dr. A. S. C. Lawrence; 
and there may be yet other forms.!° 








8 Cf. the subsequent note by A. de Bretteville, Jr., and 
J. W. McBain, Science 96, 470 (1942). 

® Francis, Piper, and Malkin, Proc. Roy. Soc. A128, 214 
(1930). 

10 Cf. the tabulation by L. H. Germer and K.-H. Stork, 
J. Chem. Phys. 6, 280 (1938), Table V. 
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Aerial Fogging of Photographic Emulsions* 
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(Received January 13, 1943) 


Alkaline hydroquinone solutions, oxidized in the presence of a suitable photographic ma- 
terial, exert a fogging action upon that material. The fogging action is due to some transitory 
intermediate in the oxidation process. Various possibilities are considered, and the conclusion 
is reached that the peroxide radical, HO, is the fogging agent. The aerial fogging action of a 
number of essentially different developing agents is explained on the common basis of this 


intermediate. 


INTRODUCTION 


XIDATION processes which do not involve 
direct addition of oxygen probably take 
place in steps of single electron transfer.1~* 
Thus, the first step in the reaction between 
* Communication No. 901 from the Kodak Research 
Laboratories. 
1P. A. Schaffer, J. Phys. Chem. 40, 1021 (1936). 


2 J. Weiss, Naturwiss. 23, 64 (1935). 
’L. Michaelis, An. Rev. Biochem. 7, 1 (1938). 


hydroquinone and oxygen should be 
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The important participation of the semi- 
quinone in the autoxidation of durohydro- 
quinone has been demonstrated by kinetic 
methods.‘ The semiquinone reacts more rapidly 
with oxygen than the hydroquinone. This 
reaction likewise should yield the peroxide 
radical, O2- (or HO2). However, the radical 
has not been detected kinetically in either the 
oxidation of hydroquinone or durohydroquinone. 


Fic. 1. ‘‘Collar effect.” 





It is well known that the oxidation of alkaline 
hydroquinone solutions in the presence of a 
suitable photographic emulsion can produce 
fog.5-8 This phenomenon has been attributed 
to chemiluminescence? and to the action of 
hydrogen peroxide.* Neither suggestion has been 
confirmed. 

The purpose of the present paper is to show 
that neither chemiluminescence nor hydrogen 


4T. H. James and A. Weissberger, J. Am. Chem. Soc. 
60, 98 (1938). 

5 J. I. Crabtree, Am. Ann. Phot. for 1919, p. 20. 

6M. L. Dundon and J. I. Crabtree, Brit. J. Phot. 71, 
701, 719 (1924). 

7E. Fuchs, Phot. Ind. 22, 27, 56 (1924). 

8C, E. K. Mees, Theory of the Photographic Process 
‘yc, Company, New York, New York, 1942), 
p. 458. 
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peroxide can account for aerial fog. The fog must 
be due to the action of a transient intermediate 
formed during the oxidation of the developing 
agent. This transient substance is probably the 
peroxide radical, HOs, since only this radical or a 
compound of similar character possesses the 
necessary properties. 


EXPERIMENTAL 


A modification of Fuchs’ experimental tech- 
nique’ was employed. The film was immersed for 
40 seconds in the solution under examination. 
The film was then removed and suspended in air 
for the desired time (usually 9 minutes). It was 
then bathed for 15 seconds in a one percent 
acetic acid stop-bath, washed in running water 
for 5 minutes, and developed for 4 minutes in a 
standard elon-hydroquinone developer. This 
procedure at least partially separated fog forma- 
tion from development effects. Experiment 
showed that the initial time of immersion was 
not critical, and could vary from 20 seconds to 
1 minute without producing noticeably divergent 
results. The aerial fog before final ‘“‘development”’ 
was usually invisible or of very low density. 
Normal fog produced during the experiment was 
corrected by means of a blank in which aerial 
fog was suppressed by the addition of pheno- 
safranin to the solution. The validity of this 
correction was checked by experiments carried 
out in the absence of oxygen. The normal fog in 
most cases was quite low, seldom exceeding 0.05. 

In additional experiments, the film was al- 
lowed to stand quietly for some specified time 
partly immersed in the test solution. It was then 
removed and processed as before. Aerial fog 
showed up in the form of a dense ‘‘collar’’ at the 
air-liquid interface. An example is shown in Fig. 1. 

The standard solution used for the production 
of aerial fog contained 0.05M hydroquinone, 
0.10M sodium sulfite, 0.10M sodium carbonate, 
and 0.005M potassium bromide. When_ the 
effect of variation in concentration of one com- 
ponent was studied, all other concentrations 
and pH were maintained essentially constant. 
All solutions and air were at 25°C. 


The Chemiluminescence Theory 


Fuchs had established that aerial fog was 
connected in some way with the oxygen oxida- 
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tion of the developing agent. He had shown fur- 
ther that none of the then known oxidation 
products could account for the effect, and that 
partially oxidized solutions did not possess 
enhanced fogging activity. Lenard and Wolf?® 
had observed a chemiluminescence during the 
aerial oxidation of pyrogallol, and Trautz!® 
had observed chemiluminescence during oxida- 
tion by hydrogen peroxide of certain other 
developing agents in alkaline solution. Fuchs 
accordingly attributed aerial fog to the photo- 
graphic action of chemiluminescence. He was 
totally unable, however, to obtain aerial fog 
when the plate was protected by an intermediate 
quartz sheet from direct contact with the 
oxidizing developer. 

In the present investigation, aerial fog densi- 
ties up to 0.5 were obtained by 5 minutes’ 
aeration in the absence of copper, and densities 
greater than 3.0 were obtained in the presence of 
copper. No trace of chemiluminescence was 
observed in any case, although the eyes of the 
observer were well accommodated to faint 
light. Film suspended at 5 mm above the surface 
of the oxidizing solutions for several hours 
showed no fog. Moreover, in experiments on the 
production of fog ‘“‘collars’’ (cf. Fig. 1), the film 
was immersed at various angles of interface. 
The line of demarcation of density consistently 
coincided sharply with the line of wetting of the 
film. No fog was obtained beyond the portion 
of the film wetted by the oxidizing solution. 

In general, the more sensitive an emulsion is 
to light action, the more prone it is to aerial fog. 
This fact was used as an argument in favor of the 
chemiluminescence theory. There are, however, 
exceptions to the rule. Emulsion ‘“‘B,”’ described 
in a previous paper" on fog formation, differs 
essentially from the motion picture positive 
material in that it was ripened in the presence 
of a sulfur sensitizer to the point of appreciable 
ripening fog formation. A standard motion 
picture positive material, “D,’’ was obtained 
which gave almost identical threshold sensitivity. 
Figure 2 gives comparative sensitometric curves 


* Ph. Lenard and M. Wolf, Wied. Ann. 34, 918 (1888) ; 
Eder’s Phot. Mitteil., p. 344 (1887). 

10M. Trautz, Zeits. f. physik. Chemie 53, 1 (1905); 
oo and P. Schorigin, Zeits. f. wiss. Phot. 3, 121 
1 T, H. James, J. Frank. Inst. 234, 371 (1942). 
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for 23 minutes’ development in the standard 
developer. The basic fog of 0.08 obtained with 
emulsion ‘‘B’’ has been subtracted from the 
total densities. The fifteenth step was the last 
visible with each emulsion. Control tests showed 
no significant difference in the reciprocity failure 
of the two emulsions. However, the standard 
hydroquinone solution yielded aerial fog of 
0.20 on emulsion “D,’’ and none on “B.” A 
copper-containing solution yielded a collar-fog 
density of 2.25 on “D,” and 0.07 on “B.” 
Evidently, the correlation between aerial fog 
and light sensitivity fails badly with this pair 
of emulsions. 

Further evidence advanced as favoring the 
chemiluminescence theory lay in the fact that 
certain desensitizers decreased or eliminated 
aerial fog. However, these same desensitizers 
decrease or eliminate the purely chemical 
fog produced by very dilute (ca. 10-§M) solu- 
tions of sodium stannite. 

An apparently genuine case of chemilumines- 
cence fog was reported by Benton.” This was 
produced on a hypersensitive panchromatic 
material being developed in a pyrogallol solution. 
It is interesting to note the differences between 
this and typical aerial fog. The chemilumines- 
cence was clearly visible to the operator. The 
fog was formed on film totally immersed in the 
solution. Finally, the fog was formed in a de- 
veloper which completely fails to give aerial 
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Fic. 2. Sensitometric curves for normal and special 


emulsions. Dev. time, 2} minutes at 25°C. Basic Fog: 
“B” =0.08; ““‘D” =0.01. Last visible step: the 15th. 


12,W. A. Benton, Brit. J. Phot. 86, 223 (1939). 
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fog under normal test conditions (in the absence 
of copper). 


Hydrogen Peroxide 


Elimination of chemiluminescence as_ the 
cause of aerial fog throws the burden upon some 
chemical entity produced during or as a product 
of the oxidation process. None of the stable or 
semi-stable products formed during the oxida- 
tion of hydroquinone-sulfite solutions produces 
significant fog. Hydrogen peroxide, which has 
a very transitory existence in such solutions," is 
known to attack normal emulsions, but no 
direct comparison has been made_ between 
hydrogen peroxide fog and aerial fog. 

Experiments were carried out in which fog 
produced by various concentrations of hydrogen 
peroxide was directly comparable with aerial fog 
produced in the same length of time and at the 
same pH. A second standard motion picture 
positive emulsion, ‘‘C,’’ was employed." Bromide 
and carbonate, but no sulfite, were added to the 
peroxide. Figure 3 gives the results for one set 
of experiments. In order to obtain hydrogen 
peroxide fog equal in intensity to the aerial fog, 
a hydrogen peroxide concentration equivalent 
to almost half of the entire hydroquinone content 
of the standard test solution is required. Since 
hydrogen peroxide cannot be detected after the 
addition of excess sodium sulfite to a hydrogen 
peroxide solution, and hence the two react 
together very rapidly, it becomes evident that 
hydrogen peroxide cannot be responsible for the 
observed aerial fog. Indeed, the addition of 
small amounts of copper salt to the hydroquinone 
solution will increase the aerial fog markedly 
beyond that given by the simple curve in Fig. 3. 


The Free Radicals 


We have established that neither chemilumi- 
nescence nor molecular hydrogen peroxide can 
account for aerial fog. Two possible intermediates 
of radical character must be considered: (1) the 
semiquinone, and (2) the peroxide radical, HOs. 

The semiquinone radical is probably capable 
of producing fog.“ However, this radical cannot 
completely account for the phenomenon of 

1% T. H. James, J. M. Snell, and A. Weissberger, J. Am. 


Chem. Soc. 60, 2084 (1938). 
4 T. H. James, J. Phys. Chem. 45, 223 (1941). 
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Fic. 3. Comparison of aerial fog and hydrogen peroxide fog. 


aerial fog. In a durohydroquinone developer, 
for example, aerial fog is decreased by the addi- 
tion of sodium sulfite. The sulfite does not 
inhibit the oxygen oxidation of durohydro- 
quinone, and there is no evidence that it reacts 
with the semiquinone at a kinetically significant 
rate. Moreover, aerial fog is not restricted to 
developing agents which can form semiquinones. 
For example, aerial fog of roughly the same 
amount was obtained both with the standard 
hydroquinone test solution and with a ferro- 
oxalate developer. Considerable aerial fog was 
obtained also with an ammoniacal cuprous so- 
lution. Neither of the latter solutions forms a 
semiquinone. 

The formation of aerial fog can best be ac- 
counted for by a reaction of the transitory 
peroxide radical. Like hydrogen peroxide itself, 
this radical possesses both oxidizing and re- 
ducing properties, and could readily form free 
silver by reaction with silver bromide. 

Three factors would largely control the rate of 
formation of aerial fog: (1) the rate of formation 
of the peroxide radical; (2) the relative rates at 
which the radical reacts with silver bromide and 
with oxidizable bodies in the immediate vicinity ; 
and (3) the formation in the developer of oxida- 
tion products which more or less tend to offset 
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the fogging action of the radical. All factors may 
vary with the composition of the fogging solu- 
tion and the nature of the developing agent 
whose oxidation gives rise to the fog. 

Pyrogallol furnishes an extreme example of the 
importance of the third factor. Solutions con- 
taining this agent are notably free from aerial 
fogging action. The pyrogallol even protects 
solutions of other developing agents. Thus, the 
addition of 0.001.M pyrogallol to a 0.05M hydro- 
quinone solution which originally gave an aerial 
fog density of 0.53 reduced the fog to 0.06; 
addition of 0.005M pyrogallol eliminated the 
aerial fog completely. The effect is probably 
due to an oxidation product of pyrogallol, 
not to the agent itself. This is shown by the 
fact that antifogging action is obtained if the 
pyrogallol is first oxidized in solution® and then 
added to the hydroquinone solution. It is not 
surprising, therefore, that no aerial fog is ob- 
tained in the oxidation of pyrogallol itself. 

In the absence of sodium sulfite, hydroquinone 
forms little aerial fog. This is probably due to the 
antifogging action of the quinone itself. The 
destruction of catalytic centers by quinone at 
pH 8.3 was observed in previous work on de- 
velopment by hydroquinone.'® A small but 
definite decrease in aerial fog was obtained by 
bathing the fogged film prior to development in a 
freshly prepared solution of quinone, potassium 
bromide and sodium carbonate. (The quinone 
decomposes rapidly in this solution.) Almost 
complete destruction of aerial fog was effected 
in a neutral quinone-bromide solution. 

The first small additions of sodium sulfite to 
the hydroquinone solution. increase the aerial 
fog, due to removal of quinone. Upon further 
additions of sulfite, the aerial fog formation 
passes through a maximum, then rapidly de- 
creases. This second stage is due to the action of 
the sulfite in decreasing the rate of oxidation of 
the hydroquinone, and hence the rate of forma- 
tion of HO. The sulfite may also react with the 
peroxide radical itself, as indicated by the effect 
upon aerial fog in durohydroquinone solutions. 

Finally, the marked increase in aerial fog 
produced by the addition of small amounts of 
copper salts to the oxidizing solution must be 


® T. H. James, J. Phys. Chem. 44, 42 (1940). 
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considered. Aerial fog has been detected when 
copper salts were added to many solutions 
(even those of pyrogallol) which normally do not 
give measurable aerial fog. Copper accelerates 
the oxidation of hydroquinone and other de- 
velopers, and it might be supposed that this is 
the explanation of the enhanced fogging action. 
However, the increase in fog is out of proportion 
to the increase in oxidation rate, so far as can be 
judged from experiments on oxidation rates at 
lower pH. Moreover, the standard fogging 
hydroquinone test solution showed a decrease 
in fogging action when the oxidation rate was 
increased by use of a pure atmosphere of oxygen. 
This effect resulted from accumulation of hydro- 
quinone oxidation products, since it was not 
present at higher sulfite concentrations. 

An explanation of the effect of copper is sug- 
gested by the following experiment. Strips of 
film were partially immersed in an ammoniacal 
solution of cuprous oxalate. The surface of the 
solution had free access to air. As soon as slight 
fog formation became evident on the film portion 
immersed in the solution, the film was removed, 
rinsed in 1 percent acetic acid, and washed for 
10 minutes. It was then ‘‘developed”’ in a stand- 
ard elon-hydroquinone developing solution. A 
heavy collar, similar to that shown in Fig. 1, 
was obtained. In two typical experiments, the 
normal fog densities were 0.24 and 0.33; the 
collar densities were 0.88 and 1.10, respectively. 
Thus, considerable aerial fog arises from oxida- 
tion of cuprous salt. 

In a developing solution of normal composi- 
tion, cupric salts are readily reduced to the 
cuprous state. The cuprous salt, which is prob- 
ably adsorbed to the silver bromide, can readily 


‘ give rise to aerial fog by reaction with oxygen, 


forming HO,. The cuprous salt is reoxidized to 
cupric in the process, but will again be reduced 
to cuprous by reacting with the developing 
agent once more. Since these reactions are quite 
rapid, a small amount of copper salt can and does 
give rise to heavy aerial fog. 
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The quantum yield of fluorescence of diacetyl, excited 
by radiation of \4047 and 44358, has been measured and 
found to be 0.145+0.03. It is approximately independent 
of pressure. On the other hand, the yield with excitation 
of 3650, as first observed by Henriques and Noyes, 
increases apparently from zero at zero pressure to nearly 
the value obtained for the longer wave-length excitations 
at a pressure of 5 cm Hg. This marked difference in pres- 
sure effect for different energies of excitation is explained 
by assuming that predissociation occurs at a level between 


those reached upon \3650 and \4047 excitations. By re- 
moving vibrational energy collisions can save for fluores- 
cence a molecule excited by 43650, but they do not play 
an important role in the case of molecules with excitation 
energy below the predissociation limit. The observed 
phenomena of diacetyl fluorescence can be explained in 
terms of a model which involves a metastable state just 
below the state reached on absorption. This model is 
compared with McMurry’s predictions of states for mole- 
cules of the type of diacetyl. 





INTRODUCTION 


XPERIMENTS on the fluorescence of 

diacetyl show that the photochemistry of 
this molecule involves a number of elementary 
processes.'~* Energy absorbed at A\3650, 4047, 
4358 partly reappears as long-lived green 
fluorescence with maximum intensity at 5120A 
and partly leads to other processes, including 
decomposition. It is important for the descrip- 
tion of what happens to the excited molecule to 
know the quantum yield of the fluorescence. 
Inferences regarding fluorescence and other 
processes can also be drawn from the effect of 
temperature and pressure on the fluorescent 
yield. 

A tentative determination of the yield was 
made in this laboratory’ with the result that at 
pressures of 4-5 cm, it was about 0.03 with 
excitation at \A3650 and 4358, but much less at 
4047. Henriques and Noyes,’ using only \3650 
excitation, found that the relative yield in- 
creased with increasing pressure. This unusual 
behavior had not been observed with \4358 
excitation.! 

To clear up uncertainties in the absolute yield, 
two new determinations of the yield have been 
made and the earlier values have been found to 
be greatly in error. The error was due to the use 


1 Almy, Fuller and Kinzer, J. Chem. Phys. 8, 37 (1940). 
940) M. Almy and S. Anderson, J. Chem. Phys. 8, 805 
1 

3F. C. Henriques and A. W. Noyes, Jr., J. Am. Chem. 
Soc. 62, 1038 (1940). 

‘Fuller, Phillips, ‘and Almy, J. Chem. Phys. 7, 973 
(1939). 


of filters in the calibration of a_ photo-cell 
against a thermopile. The thermopile readings 
were apparently affected either by infra-red 
transmitted by the filter or by heating and 
reradiation of the filters as a consequence of 
absorption of most of the incident energy. The 
sensitivity of the photo-cell used to measure the 
incident energy was thus underestimated, the 
light absorbed overestimated and the calculated 
yields came out too low and not even relatively 
correct for the different wave-lengths. A thermo- 
pile is not, unless special precautions are taken, 
a reliable instrument for measuring the ap- 
parently monochromatic energy transmitted by 
a filter in the visible or ultraviolet spectrum. 

Two new determinations of the yield have 
been made. In the first, the apparatus was un- 
changed except that a monochromator was used 
instead of filters in the calibration mentioned 
above. At diacetyl pressures of 4 to 5 cm the 
quantum yield with \\4057 and 4358 excitation 
was found to agree to within 5 percent with the 
value obtained in the second determination, 
described below, under similar conditions. 

In the second determination, described in 
detail in this paper, a new K-O photo-cell re- 
placed the Cs-O photo-cell for measuring the 
fluorescence intensity and a new thermocouple 
as well as the pile previously used were employed 
in the calibrations. The yield at \A4047 and 4358 
was found to be 0.145. At 3650 the yield at 
5 cm pressure was 0.84 of the yield under similar 
conditions at \\4047 and 4358. At A3650 the 
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Fic. 1. Diagram of optical apparatus, showing source, H-4; improvised monochromator ; fluorescence tube; slit system 
and photo-cell (.S3, S;, and F) for measuring fluorescence intensity ; photo-cell successively at P and Q to measure incident 
and transmitted exciting light. Inset shows enlargement of fluorescence tube and photo-cell F with slits. 


pressure effect found by Henriques and Noyes 
was verified. It does not occur at the longer 
wave-lengths. 

The previously proposed model for the energy 
states and processes of the diacetyl molecule? 
will account for the variation with wave-length 
of the effect of pressure on yield by the assump- 
tion of a simple natural process in the state 
reached immediately on absorption. 


EXPERIMENTAL ARRANGEMENT AND METHOD 
Preparation of Diacetyl 


The diacetyl was supplied by two firms, the 
Lucidol Corporation and Eastman Kodak. 
Each batch used was placed in a blackened bulb 
on the system, frozen down with liquid air and 
pumped for hours. After it had melted approxi- 
mately the first quarter was pumped off. The 
last quarter was also not used. 


Optical System 


The optical system is shown in Fig. 1. The 
fluorescence tube was connected to a suitable 
vacuum system from which diacetyl could be 
admitted at pressures measured accurately by a 
multiplying manometer. The source was a 
General Electric H-4 mercury lamp. The wave- 
lengths used were isolated by means of an 
improvised quartz monochromator which in- 
cluded a large quartz prism from a Steinheil 


spectrograph. The intensity and purity of visible 
spectrum were nearly as good as with filters and 
the uncertainties mentioned in the introduction 
were avoided. 

The beam of fluorescent light striking the 
photo-cell F was limited by two similar apertures, 
S; and S,, 1 cm wide by 2 cm high. The solid 
angle of fluorescence seen by the photo-cell is 
practically constant for all points between L 
and M in the exciting beam. For points outside 
LM the slit S; causes the effective solid angle to 
decrease linearly with the distance from L, or M, 
to zero at K, or N. The beam of exciting light is 
about 1 cm in diameter in the region between 
K and N. 


Instruments and Their Calibration 


The energy incident upon the tube and the 
energy transmitted were measured in turn by 
placing a GM Visitron F 2A photovoltaic cell, 
connected directly to a shunted wall type 
galvanometer, at P and Q (Fig. 1). The cell 
intercepted the whole beam in each case. The 
fluorescent energy passing through S, was 
measured by a GM 1044 TDV K-O vacuum 
photoemissive cell. It was placed in a vacuum 
chamber together with a Western Electric 
D96475 electrometer tube. The latter was used 
in a Barth balanced circuit, with an input grid 
resistance of 11,000 megohms. 
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The sensitivities of the photo-cells in their 
associated circuits, at the appropriate wave- 
lengths, were compared with that of a Moll 
vacuum thermocouple, which had a lens one 
cm in diameter for focusing the incident light 
on the thermojunction. For the photoemissive 
cell calibration a ribbon filament lamp and a 
Kipp and Zonen glass two-prism monochromator, 
followed by a corrected lens placed at its focal 
length from the exit slit, were used to produce a 
nearly uniform monochromatic beam. The beam 
was much larger than the slit S,. The ratio of 
thermocouple and photo-cell deflections (d/dr) 
was determined at a series of wave-lengths from 
4800A to 6100A. With this arrangement plus a 
rotating sector disk the linearity of response of 
the photo-cell circuit and the ratios of the 
galvanometer Ayrton shunt were also checked. 

The ratio of thermocouple and photovoltaic 
cell deflections (d,,/dy) was obtained by placing 
them alternately in a nearly uniform but slightly 
diverging beam produced by removing lens L. 
A diaphragm with a circular opening of known 
area A was placed on the photovoltaic cell for 
the calibration. The ratio of deflections was ob- 
tained at AA3650, 4047, 4358. This calibration 
was repeated with the use of a Coblentz thermo- 
pile in air. The relative sensitivity of the photo- 
cell at 14047 and \4358 was the same with Cob- 
lentz and Moll instruments and in agreement 
with information in a paper by Parker.’ At 
3650, however, the Moll couple indicated an 
increase in sensitivity of the photo-cell, due, no 
doubt, to the absorption of lens and glass en- 
velope in front of the thermojunction. The Cob- 
lentz pile gave results in agreement with an 
extrapolation of Parker’s curves. Hence the 
Coblentz pile was used as the standard in the 
3650A-4358A range and the Moll couple, whose 
greater sensitivity was very helpful in the 
calibration of the photo-cell behind S.4, was 
assumed to measure energy, independent of 
wave-length, in the overlapping range 4047A- 
6100A. 


Method of Measuring Quantum Yield 


Since Beer’s law was observed to hold, the 
energy absorbed between two planes, distant 


5A. E. Parker, Illum. Eng. 35, 833 (1940). 
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x1 and x2 from the window of the tube, is given by 
Fabs — Eince~9(e- 2"! _ e~ «P22) : 


where e~* is the fraction of incident energy 
transmitted by the window, measured for each 
wave-length. Under the conditions of the experi- 
ment, this equation can be replaced by 





Xi tX2 
Favs = Eince*(%2—%1) exp | - ae ; | 


with 
KL+MN 
X_— xX, =——-+- LM = 1.566 cm 
and 

Xi +X 





= 5.55 cm. 


The total fluorescent energy emitted between 
x1 and x2 is the energy passing through S, multi- 
plied by 42/w, where w is the solid angle sub- 
tended at the fluorescence by S, and by e” to 
correct for absorption in the exit window. This 
assumes that the long-lived fluorescence of 
diacetyl has a uniform angular distribution. 

Now if the total incident radiation produces a 
deflection D, on the photovoltaic cell (aperture of 
area A removed) the energy absorbed per sec. 
in the diacetyl is proportional to 


din XitXe2 
p,(—*) -A-e~ap(x2—x1) exp |---| 
N 


v 


If the fluorescence passing through S, produces a 
deflection Dr in the photoemissive cell circuit 
the total energy radiated per sec. is similarly 
proportional to 


dir 
Dr: (2 -area of S4-e+?-42/w. 
dr 5120A 


The above rates of absorption and emission 
could be expressed in watts if the thermopile 
deflection were also obtained in the known 
radiation field (watts per cm*) of a standard 
lamp. This is not necessary, however, in order 
to calculate the absolute quantum yield which 
is given by 


energy radiated Afiu 








energy absorbed lie. 
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Fic. 2. (A) Reciprocal yield of fluorescence when excited 
by \A4358 and 4047, as a function of pressure, for samples 
with different preliminary treatments (see text). (B) Yield 
with 43650 excitation. Dashed curve is Eq. (8a). (C) Re- 
ciprocal yields for 43650 excitation, contrasted with those 
for \A4047 and 4358, transcribed from (A). 


RESULTS 


Readings of both incident energy and energy 
getting through the fluorescence tube, as well 
as readings of diacetyl pressure, were taken 
each time the tube was filled; this gave a large 
number of data from which the variation of 


absorption with pressure could be accurately 
determined. It was found that Beer’s law was 
obeyed. Values of the absorption coefficient 
obtained from these data are: 


r a at 27°C 
(H-4 lamp) per cm Hg per cm path 
4358 0.0215 +0.0005 
4047 0.0193 +0.0005 
3650 0.0048 +0.0005 


The intensity of fluorescence was measured at 
\A4358, 4047, 3650 in a pressure range of 0.2 cm 
to 4.5 cm Hg, by the use of three different 
procedures: 

I. The tube was refilled for each pressure. 

II. The tube was filled and emptied by steps 
and readings were taken at each step. 

III. The tube was filled to the maximum 
pressure of the run, radiated strongly with 
unfiltered light from a glass H-2 mercury lamp 
for two or three hours, and then readings of 
intensity of fluorescence taken as the pressure 
was reduced in steps. The sample also was 
given 10-min. radiations from the H-2 lamp 
between steps. This procedure was used because 
it was found that the variation of yield with 
pressure depended upon the past radiation of the 
sample. 

In Fig. 2A are plotted in detail the reciprocal 
of the yield against the pressure for two runs for 
44358 and two for \4047 by procedure I (tube 
refilled at each point) and three runs by pro- 
cedure III (strong preradiation). Within experi- 
mental error, (1) identical results are obtained 
for excitation by \4047 and by 4358, (2) the 
points fall on a line for each run and all lines 
have the same intercept at p=0, although the 
slopes of the lines vary with the procedure fol- 
lowed: Several additional runs, mostly by pro- 
cedures I and II all gave results similar to those 
plotted: straight lines with different slopes but 
with intercepts within 5 percent of an average 
value of 9.0. 

Upon excitation by \3650, on the other hand, 
the yield increases sharply with increasing 
pressure but becomes nearly constant and nearly 
as large as the yield at \\4047 and 4358 for pres- 
sures above 3 cm (Fig. 2B). Procedure III was 
used since the ‘‘impurity’’ pressure effect, dis- 
cussed in the next paragraph, was then nearly 
eliminated with the longer wave-length excita- 
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tion. The variation with pressure at \3650 is in 
substantial agreement with the result obtained 
by Henriques and Noyes. Their Fig. 2 is tran- 
scribed to our Fig. 2B adjusted to the new scale 
of ordinates. The striking contrast between the 
behavior with \3650 excitation and with the 
longer wave-lengths is shown in Fig. 2C. This 
difference will be discussed in the following 
section. 

The straight lines in Fig. 2A indicate that a 
Stern-Volmer mechanism is followed. The varia- 
tion in slope rules out self-quenching but can be 
understood if we assume that a quenching sub- 
stance is present in amounts proportional to the 
pressure of diacetyl and is destroyed by strong 
radiation. Under these conditions, the yield is 
given by 


ky 
y=, 
k,+katkop 


where k,, ka, and k,p are the probabilities, for an 
excited molecule, of fluorescence, monomolecular 
competing processes, and quenching, respec- 
tively. The reciprocal yield varies linearly with 
the pressure and the intercept at p=0 is the 
reciprocal yield in the absence of impurities. 
The slope depends upon the fraction of impurity 
present and is reduced as the impurity is de- 
stroyed by radiation. 

Since O2 is known to be a strong quencher of 
diacetyl fluorescence! (1 part in 10,000 would 
account for the observed effect) its presence was 
suspected. Drastic outgassing procedures which 
in earlier experiments eliminated entirely the 
dependence of yield on the pressure, at \4358 
excitation, had no effect here. 

At all events, the quantum yield for \\4047 
and 4358 excitation was obtained from the 
reciprocal of the average value of the intercepts 
of plots like those in Fig. 2A. This average is 9.0, 
its reciprocal 0.111, which is the yield if all of 
the fluorescence were at 5120A. We estimated 
that the two other bands of fluorescence, with 
maxima at 5580A and 6100A, are each about one- 
fourth as intense as the main band. Taking into 
account the smaller size of the quanta at these 
wave-lengths and the variation in sensitivity of 
the photo-cell, we obtained a correction factor 
of 1.3. The quantum yield for \\4047 and 4358 
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Fic. 3. Schematic energy level diagrams of models of 
diacetyl molecule. Model A is discussed in detail in text. 
a and m are electronic states with many closely spaced 
vibrational levels. Transition to ¢ corresponds to mona- 
tomic processes competing with fluorescence, including 
predissociation. B is an alternative model which probably 
could also reproduce observed properties of fluorescence. 


thus comes out 0.145 with an uncertainty of 
about 0.03. This estimate of error contains a 
large contribution from the uncertainty in the 
correction factor. It does not include an error 
which may arise from the assumption that the 
response of the Moll thermocouple is independent 
of wave-length in the range from 4047A to 6100A. 

For \3650 excitation the yield extrapolated to 
5.0 cm pressure, obtained by procedures I and 
III, was in each case 0.84 of the yield obtained 
with 4047 excitation under similar conditions. 
Earlier results with a Cs-O cell indicated a ratio 
of about 0.9, at pressures of 4 to 5 cm. 

The results can be summarized as follows. 
When excited by A4358 or 4047 diacetyl 
fluoresces with a quantum yield of 0.145+0.30. 
In the pure substance this value is nearly 
independent of pressure in the range 0-5 cm. 
With excitation of \3650 the quantum yield 
rises with pressure, apparently from zero at 
zero pressure. With increasing pressure it ap- 
pears to be approaching a yield nearly as great as 
obtained with the longer wave-length excitation. 

Measurements of the variation in yield with 
temperature, of samples which had received a 
strong preradiation, gave the result that the 
yield decreased almost linearly with rising 
temperature. The yield at 80°C was one-half 
the yield at 30°C. 
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DISCUSSION OF RESULTS 


The model previously proposed? to account 
for the properties of diacetyl fluorescence can 
readily be adapted to account for the new experi- 
mental results. These are included in the follow- 
ing list: 


1. Quantum yield with 4358A- and 4047A- 
excitation equal to 0.145 at 27°C and inde- 
pendent of pressure. 

2. Quantum yield with 3650A-excitation rises 
apparently from zero at zero pressure to 84 per- 
cent of the yield for longer wave-lengths at 
5-cm pressure. . 

3. VY at 25°C/Y at 65°C=1.65 (4358A ex- 
citation). 

4. Decay of fluorescence is simply exponential, 
with a mean lifetime (7) of 1.40107 sec.,® 
independent of pressure, exciting intensity and 
exciting wave-length.” 

5. 7 at 25°C/r at 65°C =2.0+0.4. 

6. Lifetime of state reached on absorption is 
10-5 sec., estimated from integrated absorption 
coefficient.? 


The model is shown schematically in Fig. 3A. 
An alternative appears in Fig. 3B. Assuming 
model A, consider what happens when a molecule 
absorbs a quantum of \4358 or \4047 and is 
raised to vibrational state b or v of electronic 
state a. Redistribution of vibrational energy and 
collisions bring the molecule to a different and 
lower vibrational state 0, from which it can, 
(1) fluoresce at the rate k,, (2) predissociate or 
dissociate into chemical products (c) by attain- 
ing the necessary critical configuration, at the 
rate k,.,7 or (3) go into the long-lived state m 
at the rate ko». If it goes into m it can return to a 
at a slow rate k». and then has again a chance to 
radiate. The yield is determined by the competi- 
tion between k,, and k, since it is assumed for 
simplicity that no molecules are lost from m 
except by return to o. As will be shown below 
the lifetime depends on several factors, but is 
made much longer than the reciprocal of k, +o. 
by a large ratio of Rom to Rmo. 


®R. D. Rawcliffe, Rev. Sci. Inst. 13, 413 (1942). 

7In place of, or in addition to, this process deactivation 
to the ground state by “internal conversion” could occur 
(cf. ref. 8). It is known, however, that diacetyl is decom- 
posed into gases non-condensible in liquid air by radiation 
of 14358, as well as by the shorter wave-lengths. 





When a quantum of \3650 is absorbed, the 
molecule finds itself in state u, above c, which 
we shall assume is a level at which the molecule 
can predissociate. In this case, the molecule can 
(1) be stabilized toward state o by collisions and 
follow the course outlined for 4358A excitation 
or (2) predissociate at a high rate directly into 
c and be lost as a diacetyl molecule. The low 
fluorescence yield at low pressure is ascribed to 
the loss of excited molecules to c before sufficient 
collisions have occurred to save them for a 
chance at fluorescence. 

It may be remarked that the best data for 
44047 (Fig. 2A, procedure III) show the same 
pressure effect as 3650 but at much lower 
pressures, where the error of measurement is 
nearly as large as the effect. If the departure of 
the 4047 plot from a straight line is real it 
means that collisions also reduce the chance of 
transition:to c with this excitation, but that the 
chance of this transition is about ten times less 
than with \3650 excitation. 

These processes can be described more ex- 
actly. All molecules excited by \\4047 and 4358 
arrive by collisions at state 0. Let nm, be the 
population of electronic state A, vibrational 
state o. Let m» be the population of m. Then 


dn,/dt= No— (Roct+Rom+ky)NotRmotm, (1) 
dNm/ dt = Romno— RmoNm, (2) 


where NV, is the number of molecules excited 
per sec. At equilibrium each rate is zero. It 
follows upon adding Eqs. (1) and (2), each set 
equal to zero, that the quantum yield (Y), which 
is given by k,n,/N, is 


Y=k,/(ky+Roc). (3) 


Since, from experiment, Y=0.145 and k,=10°, 
Eq. (3) leads to kp. =6X10*, at room tempera- 
ture. If we assume that k,.=Z,. exp [—E../RT], 
the observed temperature variation of Y yields 
the values, Z,.=8.7 10° sec.-' and E,,=4400 
cal./mole. 

The transient solution of Eqs. (1) and (2) for 
No is, 


No= C1 exp (Ait) +c2 exp (Aol), (4) 
where 
Rmo(Roc +k») Rmo(Roc tk») 
ee ee == a ee 
B B 
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B=Roctky+Romtkmo. In obtaining \; and dg, it 
is assumed that {6?—4kmo(Roc+k,)}? can be 
adequately replaced by the first term in its 
expansion. It appears in the following paragraphs 
that the rates are such that the approximation is 
a very good one. 

The fluorescent intensity (=k,n.) is observed 
to decay exponentially at a rate of about 10 
sec.—!. Xs, however, is of the order of magnitude 
of 10°. Hence c.“ci; otherwise an initial rapid 
drop, before the observed slow decay, would 
occur. The fluorescence mean life must, there- 
fore, in this model, be, 


T= —1/d1=B/Rmo(Roctk,). (S) 


Substituting for 8, setting r= 1.410, k,. +k, 
=7X10°, Eq. (5) becomes 


1 Rom 
1.4 10-*=—-+ 


eaaeeniatemm, (Sa) 
Rmo 7X10°Rmo 


in which a term or order 10-6 has been dropped. 
Since both terms on the right are positive, 
Rmo>7 X10? and Rom/Rmo< 10°. 

By applying the boundary conditions at -=0, 
namely, 2o=¢Ci+¢2 and Mm=NKRom/Rmo, one finds 
that 


C2 Roctky 
— =——_., (6) 
C1 Romt+ mo 


Since from experiment this ratio is small, say 
less than 0.1, it follows that komn>7X10® But 
Rom/Rmo<10® and, therefore, kmo>7 X10, ten 
times the lower limit deduced from Eq. (5a). 
Hence 1/km. can be dropped from Eq. (5a) and 
it follows that, approximately, Rom/Rmo= 10°. 

Thus lower limits and a ratio are fixed for 
Rom and kmo. Since they measure the rates of 
radiationless transitions in and out of the me- 
tastable state they probably depend upon the 
amount of vibrational energy in the initial state. 

For molecules excited by 3650 the analysis 
given holds except that N, in Eq. (1) is the 
number reaching state o by collision, against the 
competition of immediate transition to c. If the 
number reaching u per sec. is NV, the fraction of 
these reaching 0 is 


No/Nu=ksp/(RuctRsp), (7) 


where k,p is the chance for a stabilizing collision 
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at pressure p. Substituting for N, in Eq. (1) 
and solving for the yield as before one finds 
kno k, 


kp kp 
Ny RetRocRucthsb Ructhep 








(8) 


Y, is the limiting yield at high pressure and 
should be equal to the yield at \\4047 and 4358, 
which was 0.11, before correction for distribution 
in wave-length. From measurements by Ander- 
son? on the diffusion of excited diacetyl mole- 
cules, kinetic theory gives a cross section of 
11.310-* cm? or a collision diameter of 5.9A. 
Assuming the same diameter is required for a 
collision in which enough vibrational energy is 
transferred to stabilize a molecule in state u, 
the number of collisions per excited molecule 
per sec. comes out to be 1.5X10%, with p-in 
cm Hg; that is to say, ks=1.5X108. 

We can estimate k,,. from the experimental 
curve. The fluorescence yield, after a small 
correction is made for the type of pressure 
quenching observed with \\4047 and 4358, is 
0.055 (half of Y..) at p=0.95 cm. Hence, k,,.= 1.4 
X 108. Equation (8) becomes 


1.5p 
Y=0.110-———_—-. 


1.441.5p 


(8a) 


It is plotted as a dotted line in Fig. 2B. It fits 
the observed fluorescence curve well if the latter 
is assumed to be depressed by the “impurity”’ 
pressure effect observed under similar conditions 
with the longer wave-lengths. In fact, Y/Y, 
calculated from (8a) is, at p=5.0 cm, equal to 
0.84, which is the ratio found experimentally 
between the yields at 43650 and at \A\4047 and 
4358, under similar conditions. The experimental 
curve of Henriques and Noyes can be fitted with 
the same type of equation if ku. is set equal to 
2.510%. These authors find that if the total 
pressure is raised by adding acetone the curve 
passes through a maximum, due presumably to 
the superposition of a weak quenching by ace- 
tone on the main pressure effect. 

The lifetime for predissociation is 1/k,. or 
5X10-* sec. Although this is much longer than 
the vibrational periods of the molecule are likely 
to be, it is not an unreasonably long time for a 
polyatomic molecule to get into the critical 

















configuration which may be necessary for 
predissociation. 

The mechanism shown in Fig. 3B could doubt- 
less also account for the observed properties of 
the fluorescence. In this case the fluorescence is 
the transition m—X, with decomposition or 
other monatomic process m—c competing. The 
level at which m dissociates into c may or may 
not be the same as the level into which a dis- 
sociates. 

In either scheme, A or B, the essential special 
features are (1) a long-lived state m at about the 
same energy level as state a, (2) a predissociation 
limit below vibrational level u to account for the 
differences observed with A3650 as against 
\A4047 and 4358 excitations, (3) predissociation 
or other monomolecular decomposition or deac- 
tivation processes, competing with the fluores- 
cence, whatever the exciting wave-length, at a 
rate which depends on the temperature. 

Henriques and Noyes have proposed a mecha- 
nism which will also account for the effect of 
pressure on yield with A3650 excitation. A 
state B® (our a) is reached on absorption which 
may dissociate or go reversibly into B' (our m). 
B' may be transformed by collisions into B? 
which may fluoresce. It is similar to our model 
if B* is identified with a low vibrational level o 
of electronic state a. The essential difference 
is in the role of collisions. In their model, colli- 
sions enhance the transfer to the final fluorescing 
state (B'—B?); in ours, collisions prevent im- 
mediate predissociation upon excitation by 
\3650. It is certainly easier with our model to 
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give a simple and natural explanation of the 
qualitative difference in the effect of pressure 
with 43650 excitation as compared with excita- 
tions by longer wave-lengths. 

The state m may correspond to a tautomeric 
form of excited diacetyl, of the kind suggested 
by Franck and Livingston® to explain certain 
properties of dyes. It may, however, be a state 
which is metastable against electronic transition 
in the usual sense. The electronic states of 
conjugated ketones have been discussed by 
McMurry® and it appears that in molecules 
containing two conjugated 


C=0 
al 


groups, such as diacetyl, the first excited states 
come in pairs (called EZ’ and E” states) which are 
nearly degenerate. Transition to only one of a 
pair is allowed by the electronic selection rules. 
McMurry states, however, that both transitions 
probably do occur and that both owe most of 
their intensities to perturbation due to electronic- 
vibrational coupling. But if we identify states a 
and m in our model with the two states of a 
pair (E;’ and E,’’) and suppose that on account 
of the electronic selection rules the transition 
a—X is more than one hundred times as prob- 
able as m—X, the model proposed fits qualita- 
tively with the theoretical picture. 


8 J. Franck and R. Livingston, J. Chem. Phys. 9, 184 
(1941). 
9H. L. McMurry, J. Chem. Phys. 9, 241 (1941). 
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The conductivity of sintered pellets of lead chromate varies with the temperature in accord- 
ance with the Boltzmann-type law ox exp (—1.48 ev/kT) in the temperature range 700°C to 
350°C and with the law ¢= exp (—0.60 ev/kT) at lower temperatures. At constant temperature 
the conductivity increases with decreasing oxygen pressure, ¢ « P9~!/3-1, which indicates that lead 
chromate probably dissociates in accordance with the equations: 2PbCrO,—~Pb(CrO2z)2+ Pb® 
+202 and Pb°—Pb*+e. The energies of these reactions are 5.28 and 1.20 electronvolts, re- 
spectively. These same reactions would explain optical and thermal darkening as well as the 
conductivity, and they indicate along with transport measurements and Seebeck effect measure- 
ments that lead chromate is an excess electronic semiconductor. 





I. INTRODUCTION 


HROME yellow pigments, which consist 

essentially of lead chromate, have the 
unfortunate property of darkening on exposure 
to light. It was expected that lead chromate is one 
of that interesting group of photo-sensitive 
materials, which includes the silver halides and 
zinc sulfide, that decompose photolytically. 
With this possibility in view the work of the 
present paper on the electrical properties of 
pure lead chromate was undertaken. It was 
hoped that from this work would come an under- 
standing of the fundamental chemical reactions 
that produce darkening. In addition to the 
present paper, studies of the optical and elec- 
trical properties of variously treated lead chro- 
mate pigments are being carried out at this 
laboratory, and the results will be published in 
other papers. 

Lead chromate belongs to that group of 
solids known as semiconductors. The solids 
classified as semiconductors are those which 
have ‘‘a small electronic conductivity that is 
negligible at very low temperatures and in- 
creases with increasing temperature.”! The 
electrical conductivity of a semiconductor is 
very much a property of the method of measure- 
ment. In determining the conductivity of a 
semiconductor, therefore, certain problems must 
be considered: the contact with the semicon- 


* A dissertation in physics presented to the Faculty of 
the Graduate School in partial fulfillment of the require- 


ments for the degree of doctor of philosophy. 
1F, Seitz, The Modern Theory of Solids (McGraw-Hill 
Book Company, Inc., New York, 1940). 
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ductor; the effect of impurities and of the 
mechanical and thermal history on the conduc- 
tivity; the dependence of the conductivity on 
voltage, temperature and oxygen pressure; and 
the determination of the type of conduction. 
These problems are discussed in the order given - 
in the section on results. 


II. PREPARATION OF SAMPLES 


Pure lead chromate pigments were prepared 
for this work by the Krebs Pigment Department 
of E. I. du Pont de Nemours and Company. 
They started with C.P. sodium chromate and 
C.P. lead nitrate and obtained by carefully 
controlled fractional crystallization very pure 
salts from which the lead chromate was obtained 
by precipitation. A spectroscopic analysis made 
at the du Pont experimental station showed that 
the samples contained possibly ‘‘several parts in 
10,000 of Si, Na, and Ca; traces perhaps picked 
up in analysis of Mg and Ba; and no other 
metallic impurities.”” In all of the work carried 
out here the structure of the material used was 
restricted to the monoclinic phase which is 
stable from room temperature to 707°C. 

The pigment was pressed into pellets about 
one-half inch in diameter and one-quarter inch 
thick with the use of a press made from an eight- 
ton hydraulic truck jack. Although a die made 
of hardened tool steel and equipped with two 
hardened plungers was found satisfactory for 
making the pellets used in. rough measurements, 
considerable difficulty was encountered because 
the pigment stuck to the die. In order to reduce 
this and to eliminate the strains produced when 
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removing the pellet from the die a special 
divided die was constructed of extremely hard 
material. About 2 grams of the pigment were 
used to make each pellet, and the optimum 
forming pressure employed was 9000 Ib./in.?. 
The density of the pellets thus formed was about 
50 percent of the single crystal value. It may be 
added that natural single crystals from Siberia 
and Tasmania have a density of 5.92 and 6.00 
grams/cm', respectively, whereas artificial crys- 
tals have a density between 6.1 and 6.3 g/cm*.? 

On exposure to the atmosphere the lead 
chromate pellets were found to absorb up to 0.1 
percent of their weight in moisture. This amount 
of moisture increased the conductivity a hun- 
dred-fold.. If, however, a current was allowed to 
pass through the sample for several hours, the 
conductivity dropped below the value for a 
carefully dried specimen. Subsequent examina- 
tion of the sample showed dark areas at the 
contacts. The conduction increased when these 
areas were removed. For standardization, all 
of the pellets used in the final measurements 
were dried at 105°C for approximately 24 hours. 

Contact was made to the sample either with 
platinum disks which were held against the 
pellet with a strong spring, or with gold films 
that were vacuum-evaporated onto the parallel 
surfaces of the pellet. Platinum disks were usually 
pressed against the gold layer. 

In order to measure the contact resistance 
some of the samples were formed with probes. 
Two small flat spirals of 1-mil Nichrome wire 
were imbedded between three 1-gram sections 
which were formed into a single pellet. The pellet 
was pressed in stages, each section being pressed 
to 2000 Ib./in.? before the spiral and next section 
were placed on it. The entire unit was completed 
at a pressure of 9000 Ib./in.?. 

For lead chromate heated 18 hours at 400°C 
in air the density becomes 5.8 grams/cm* or 93 
percent of single crystal density. Pellets thus 
sintered possess the advantage that many of 
their properties are nearly those of a single crys- 
tal. For some measurements the temperature 
and oxygen pressure during sintering were 
varied and the sample was quenched to room 
temperature. 


(1939) Grossner and F. Mussgnug, Zeits. f. Krist. 75, 410 
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Ill. APPARATUS AND CIRCUITS 


Most of the measurements were made with the 
current-voltage circuit shown in Fig. 1(a). 
To prevent puncturing of the sample as a result 
of Joule heating and the consequent runaway 
conduction, the ballast tube 6D6 was used as a 
series resistor with its grid bias adjusted so as 
to hold the current below 5 milliamperes. 

In order to measure the contact resistance, the 
potentiometer circuit of Fig. 1(b) was added to 
that of Fig. 1(a). The tube 6F5 is operated as an 
impedance changer, and was selected because 
of its large grid resistance and large transcon- 
ductance. The unknown potential difference 
appears across the 1 megohm cathode resistor, 
and is measured with a vacuum tube voltmeter 
(VV). 

Two furnaces were used to determine the 
temperature and oxygen dependence. The high 
temperature furnace (Fig. 2) is a modification 
of that described in the thesis of M. Earle.* 
The Alundum tubes (A), Glowbar heating ele- 
ments (G), and quartz tube (Q) for the vacuum 
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Fic. 1. (a) Upper, circuit for measuring current-voltage 
characteristics. (b) Lower, circuit for measuring contact 
potential drop. Li and Le are two lucite commutator 
switches coupled so as to keep the grid of the tube 6F5 
positive when reversing the current through the sample. 
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Note: w=ohm, 2=megohm. 


3M. D. Earle, Phys. Rev. 61, 56 (1942). 
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Fic. 2. High temperature furnace. 


or oxygen chamber are the same as in his design. 
Only the reading thermocouple and sample 
holder are changed. These are inserted through 
the top end of the quartz tube. A 32-inch diameter 
thermocouple tube 18 inches long (T) is waxed 
into a tightly fitting hole in the 1-inch thick 
brass cap. Two lavite blocks (L) are sealed to the 
thermocouple tube with Alundum cement. A 
third lavite block is sealed to the lower end of a 
12 inch tube of the same diameter. This second 
tube can be moved up and down parallel to the 
first, during which motion the lavite blocks 
serve as guides. They also serve as supports for 
the platinum electrodes. The weight of the 12- 
inch tube is sufficient to give the mechanical 
contact. The leads from the sample are taken 
through the thermocouple tubes. The resistance 
of both lavite and thermocouple tubing (Alun- 
dum) was found to be very high compared with 
the resistance of the sample at any given 
temperature. 

The second furnace is useful in the temperature 
range from room temperature to 220°C. Ni- 
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chrome wire replaces the Glowbar heating ele- 
ment and a Pyrex tube replaces the quartz tube 
of the foregoing furnace. The electrodes are brass, 
supported by a stainless steel tube and insulated 
by Pyrex rods. 


IV. RESULTS 
A. The Electrical Contacts 


As already mentioned, a moist sample through 
which current is flowing develops on its surface a 
layer possessing high resistance. The probe 
method indicates that vacuum-evaporated gold 
electrodes show no contact resistance when 
placed on a thoroughly dried sample, whereas 
vacuum-evaporated aluminum electrodes show 
a contact resistance. For example, at 105°C the 
resistance varied from a value equal to three 
times the volume resistance of the sample at 5 
volts to about 5 percent of the volume resistance 
at 50 volts. On the other hand, at 105°C a 
contact made mechanically with a platinum 
electrode has a contact resistance which is 
independent of the voltage and is 5 times that of 
the volume resistance of the sample. Mechanical 
electrodes of lead foil show the same high con- 
tact resistance as platinum. The method is not 
too accurate since the exact position of the probes 
is difficult to ascertain. Also in these porous 
unsintered samples the evaporated gold possibly 
penetrated the surface of the sample. See Table I. 

For sintered pellets at 105°C platinum and 
Nicrome mechanical contacts have about twice 
the resistance of evaporated gold or lead foil 
contacts. Aquadag contacts, on the other hand, 
seem to have half this resistance, perhaps be- 
cause Aquadag introduces moisture into the 


TABLE I. Voltage drop across sections of unsintered 
sample at 105°C. 











1-2 2-3 3-4 1-4 
Pt-Pt 5.5 0.7 3.8 10.0 
mechanical 33 3.8 23.1 59.2 

3.7 0.6 3.0 7.0 
Au-Au a2 3.8 2.9 9.8 
evaporated 16.1 18.2 13.8 47.6 
Pt-Au 
mech. evap. 44.5 2.8 iJ 49.1 
Al-Al 2.45 0.40 2.40 5.25 
evaporated 17.2 14.6 17.1 49.1 
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sample. The probe method, which for the set-up 
used becomes accurate only at higher tempera- 
tures (150°C) when the resistance of the sample 
is sufficiently low, shows that evapoated gold 
electrodes have no contact resistance. Platinum 
mechanical electrodes, however, have a high re- 
sistance until about 500°C when they become as 
good as gold electrodes. 


B. The Effect of Impurities on the Conductivity 


Moisture increases the conductivity, the elec- 
trolysis of the water contributing the added 
conduction. 

A series of samples were measured at several 
temperatures below 200°C. The material from 
which the samples were made contained metal 
oxides added after precipitation of the lead 
chromate. The conductivity of the treated 
material was found to decrease by as much as a 
factor of 10,000 and depended directly on the 
amount and type of oxide added. This behavior 
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Fic. 3. Conductivity as a function of temperature. The run 
either up or down took 80 minutes in each case. 
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Fic. 4. Resistance-temperature curve for 
oxygen equilibrium. 


is readily explained by assuming that the oxides 
form layers or coatings of high resistance on the 
particles of lead chromate. It was found that the 
light-fastness of these materials in general 
increased with decreasing conductivity, showing 
that there is probably a connection between the 
two properties. 


C. The Effect of the Thermal History on the 
Conductivity 


An unsintered sample has a conductivity of 
1.2X10-5 ohm cm at 105°C. The conductivity 
of a sintered pellet, on the other hand, also at 
105°C, varies with the sintering temperature and 
oxygen pressure. In fact values that range from 
10-® to 10-* ohm cm have been obtained. 
The value of 5X10-° ohm-! cm found by 
Kapp‘ for a natural single crystal lies within this 
range. 


D. The Conductivity as a Function of 
Voltage 


The dependence of the conductivity of an 
unsintered sample on the voltage is related to the 
contacts, the temperature, and the impurities. 
Gold electrodes extend the range of Ohm’s law 


4G. Kapp, Ann. d. Physik 22, 257 (1935). 
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to lower voltages and the temperature and 
impurities have a marked effect on the high 
voltage end of the range. On the other hand, a 
carefully prepared sintered pellet shows none of 
the foregoing behavior, and is closely ohmic over 
a large range of voltage. 


E. The Conductivity as a Function 
of Temperature 


The dependence of the conductivity « on the 
absolute temperature T is given for certain 
temperature ranges by the Boltzmann-type 
equation 


o=Ae- kT 


where k is the Boltzmann constant, ¢ a constant 
called the activation energy, and A is relatively 
independent of T. The name ‘‘activation energy” 
will be restricted to k times the slope of the 
log o versus 1/T curve. The significance of the 
activation energy will be discussed in the next 
section. 

The most significant data concerning tempera- 
ture dependence have been taken on samples 
sintered at 700°C. In this case the activation 
energy is 1.48 ev, as shown in Figs. 3 and 4, in 
the temperature range from 700°C to about 
350°C when the rate of temperature change is 
sufficiently slow so that oxygen equilibrium is 
maintained. At lower temperatures the <activa- 
tion energy fluctuates widely about the value 
0.60 ev. During very quick runs in which the 
specimen is taken from 700°C to 375°C in about 
4 minutes, the activation energy is 0.65 ev 
below 20 cm of Hg oxygen pressure and is 0.96 
ev in some samples at higher pressures. 

As already mentioned, the conductivity at 
low temperatures depends very much on the 
temperature and oxygen pressure at sintering. 
These factors affect the log o versus 1/T curves 
on the initial heating of a quenched specimen 
from room temperature in the manner shown in 
Fig. 3. It may be seen that there is a dip, or at 
least a plateau, in the curve at about 300°C. 
The precise character of this dip depends upon 
both the temperature and oxygen pressure at 
sintering. The low temperature conductivity, 
despite variations, is in all cases lower after 
heating than before. 
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F. The Conductivity as a Function of 
Oxygen Pressure 


There are a few semiconductors in which at 
higher temperatures there is observed an in- 
trinsic electronic conductivity.’ It frequently 
happens, however, that at least part of the 
conductivity depends on the surrounding gas 
atmosphere.® In such cases the conductivity is 
related to the presence of a stoichiometric excess 
of one of the constituent atoms. This concentra- 
tion of stoichiometric-excess atoms in the semi- 
conductor is proportional to some power of the 
partial pressure of one of the constituents of the 
semiconductor in the surrounding atmosphere. 
In excess semiconductors in which there is 
conduction by electrons this is a negative power, 
whereas in defect semiconductors in which holes 
conduct it is a positive power. The particular law 
which applies in a given case gives a hint con- 
cerning the chemical reaction which produces the 
conducting units in the semiconductor. In zinc 
oxide, for example, the conductivity varies 
inversely as the fourth root of the oxygen pres- 
sure. This indicates that zinc atoms are formed 
and oxygen molecules released. In the case of 
zinc oxide, however, it is not possible to tell 
immediately whether nearly all or only a few of 
the zinc atoms lose their electrons to the conduc- 
tion band, since both possibilities lead to the 
one-fourth law.’ 

The case of lead chromate is more complicated 
than the usual case of the metallic oxide, 
because a chromate radical replaces the simpler 
oxygen ion. However, the steps used to compute 
the oxygen dependence associated with an 
assumed reaction are the same as in the simpler 
cases. As we shall see, there is reason for assuming 
that a small number of lead chromate molecules 
dissociate according to the reaction: 


2PbCrO,—Pb(CrO2)2+ Pb°+ 202 
or 


Pb*+++2Cr0O,--—2CrO.—+ Pb®+ 20>. 
There are formed in this way mp» neutral lead 


5H. Diinwald, H. H. von Baumbach, and C. Wagner, 
Zeits. f. physik. Chemie B22, 226 (1933). 

6K. Baedacker, Ann. d. Physik 22, 749 (1907); 29, 566 
(1909) ; Baedacker, Physik. Zeits. 9, 431 (1908); 13, 1080 
(1912); C. Wagner, Physik. Zeits. 36, 721 (1935); B. 
Gudden, Ergeb. d. exact. Naturwiss. 13, 223 (1934). 

7™P. H. Miller, Phys. Rev. 60, 890 (1941). 
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containing vessel ; see F. Seitz, reference 1, p. 462. 





CONDUCTIVITY OF LEAD CHROMATE 201 





' 4 10 40 100 
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Fic. 5. Resistance-pressure curves. 


atoms which probably occupy interstitial posi- 
tions. 2» may be calculated as follows. 

We make the assumption that the crystal is in 
thermodynamic equilibrium at the temperature 
T. This requires that the free energy shall be a 
minimum. The free energy A is made up of 
the following terms: (1) the energy required to 
produce the above reaction; (2) the free energy 
of the vacancies in the positive ion lattice; (3) 
the free energy associated with the chromite, 
CrO.-, ions; and (4) the free energy of the oxygen 
molecules. Thus 


A=npent+kT ny In (no/N)+2kTn, In (20,/N) 
+kT(No+2nz) In No+2n,/C, 


where ¢; is the energy of the foregoing chemical 
reaction, N is the number of molecules of lead 
chromate composing the lattice, No is the num- 
ber of molecules of oxygen in the gas surround- 
ing the crystal, and C is a constant.’ It follows 
that for the free energy to be a minimum n, 
must satisfy the following condition 


ny «= No exp (—,/3kT) « Po! exp (€/3RkT) 


since Np is proportional to the partial pressure of 
the oxygen and provided that 2.<Np. 

At a particular temperature, 1; of these m» lead 
atoms will be ionized: 


Pb°—Pbt+e. 


8 C is approximated by (2xmkT)!V/h® where m, k, T, 
and h have the usual meaning and V is the volume of the 





The number ny; of free electrons thus produced 
may be calculated with the assumption of 
thermodynamic equilibrium following the pro- 
cedure outlined above. This leads to the condition 


n;<n»' exp (—e,/2kT) 


provided n;<n,. Therefore, since the conduc- 
tivity o is proportional to m,, 


o x Pye #7 


where e=€;/2+¢,/6. 

We now see the significance of the activation 
energies of the previous section. They consist of 
two terms. One arises from the energy necessary 
to produce the chemical reaction and the other 
comes from the energy necessary to produce the 
ionization reaction. Both of these terms enter 
when the rate of temperature change is such as to 
maintain oxygen equilibrium, so that we must 
set e,/2+«,/6 equal to 1.48 ev. When the tem- 
perature is changed very rapidly, however, the 
slower chemical reaction is suppressed to a great 
extent, so that the activation energy depends 
mainly on the ionization reaction. Therefore 
le, lies between 0.96 and 0.65 ev. In the opinion 
of the writer the lower value is the most reliable. 
This is indicated by the fact that at low tem- 
peratures where the chemical reaction is almost 
completely suppressed the activation energy is 
about 0.60 ev (Fig. 3). Using 0.60 ev for 4e,; it 
follows that e, is 1.20 ev and ¢ is 5.28 ev. 

If the assumption that lead chromate dis- 
sociates according to the equations given above 
is valid, we should expect lead chromate to be 
an excess electronic conductor and should there- 
fore expect that the conductivity would decrease 
with increasing oxygen pressure. This in fact 
is found to be the case. Measurements of the 
oxygen dependence of the conductivity were 
made around 600°C. The results are shown in 
Fig. 5. It is found that the conductivity obeys 
the law: 

g x Py-t8-1 


which is very close to the foregoing theoretical 
law, namely: 
ca Pa. 


The reaction which leads to this theoretical law 
was chosen from a number of possible reactions 
as being the most likely. It should be noted that 





lead chromate becomes green on being heated 
in the absence of oxygen. This color must be 
assumed to come from trivalent chromium, 
which supports the preceding reaction. 

It is of interest to note that if most of the 
interstitial lead atoms are ionized, thermody- 
namic considerations would predict 


ox Py}, 


G. The Ionic or Electronic Nature of — 
the Conducting Carriers 


The foregoing discussion indicates that lead 
chromate is a normal electronic conductor. 
Transport measurements of the type employed 
by Tubandt?® have been made to test this directly. 
The results indicate that the ionic conduction 
must be less than 1 percent of the total conduc- 
tion. It must be admitted that there is a slight 
possibility that blocking layers between the 
pellets or between the individual grains pre- 
vented the ions from migrating to the electrodes. 
This seems very unlikely, however, since Tu- 
bandt’s methods have proved successful in so 
many cases.!° 

Hall effect measurements have also been 
attempted. However, the results to date have 
not been reliable because of the high resistance 
of the samples. 

Seebeck effect measurements, on the other 
hand, indicate clearly that lead chromate is an 
excess semiconductor. From qualitative measure- 
ments the Seebeck e.m.f. per degree is in the 
neighborhood of 3 millivolts at 50°C. The hotter 
electrode was positive, so that the conductors are 
electrons. As a check, the effect was measured in 
Cu,0O, which is known to be a hole conductor, 
and it was found that the hotter electrode was 
negative as it should be. It is interesting to note 
that the free electron model preducts that the 
Seebeck e.m.f. should be equal to $¢;/T. This 
gives for the Seebeck e.m.f. at 50°C the value 


°C. Tubandt, Handbuch der Experimentalphysik (1932), 
Vol. 12, p. 396. 

10 In support of the conclusions drawn here is the recently 
published article by Selma B. Brody in J. Chem. Phys. 10, 
650 (1942) in which it is concluded from x-ray investiga- 
tions that lead chromate has a close packed structure and 
therefore ionic migration is unlikely. 
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$(1.20/325) or about 2 millivolts which is in 
rough agreement with the measured value. 


H. Summary of Results 


In the measurement of the electrical conduc- 
tivity of lead chromate the best results have been 
obtained with the use of sintered pellets pos- 
sessing gold electrodes that were evaporated on 
in vacuum. (1) Such samples have a negligible 
contact resistance and obey Ohm’s law over a 
large range of ‘voltage and temperature. (2) 
The conductivity has been found to vary with 
the temperature in the range from 700°C to 
350°C in accordance with the Boltzmann-type 
law o <exp (—1.48ev/kT), whereas the conduc- 
tivity obeys the law o«exp (—0.60ev/kT) at 
lower temperatures. (3) At constant temperature 
the conductivity increases with decreasing oxy- 
gen pressure in accordance with the relation 
gx P,—/%1t, This relation indicates that lead 
chromate dissociates in accordance with the 
equations: 


2PbCrO.—Pb(CrO,z) ot Pb®°+ 202 
Pb°—Pbt +e. 


The activation energies for these reactions are 
probably 5.28 and 1.20 electron volts, respec- 
tively. These same reactions may explain the 
optical darkening process with the optical 
excitation replacing thermal excitation in produc- 
ing the darkening substance. (4) Hall effect 
measurements on sintered pellets have not been 
reliable ; however, transport measurements show 
that ionic conduction is probably less than 1 per- 
cent and Seebeck effect measurements indicate 
conduction by electrons with a Seebeck e.m.f. 
of about 3 millivolts per degree. (5) It is con- 
cluded that lead chromate is an excess electronic 
semiconductor. 

In conclusion I wish to thank Professor F. 
Seitz, Dr. P. H. Miller, and Dr. A. W. Lawson 
for suggesting this problem and for their helpful 
advice and criticism. I am also indebted to W. 
Feldman and members of Works Progress 
Administration Project 25792 for assistance in 
observation and the Krebs Pigment Department 
of E. I. du Pont de Nemours and Company for 
financial assistance. 
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The correlation effect on the molecular distance distribution in a Bose-Einstein gas is cal- 
culated as a function of volume and temperature, particularly for the case of degeneracy 
(condensation) of the gas. It is shown that, under certain conditions in the condensed state, 
the density as seen from any one molecule appears, over a rather extended region around it, 
to be considerably larger than the true density. This result seems to have direct relation to the 
fact that the low temperature modification of liquid helium has a negative coefficient of 
thermal expansion. 





1. INTRODUCTION 


HLENBECK AND GROPPER! have shown that the symmetry properties of the wave functions 
of an ideal Bose-Einstein and Fermi-Dirac gas involve peculiar distance correlations between 
the respective force-free gas molecules. In the limiting case of high temperatures they found that 
the probability of finding a particle within the distance interval (r, 7+dr) of a chosen particle is 
given by: 
(N/V) {1+exp (—22r?/n?) } 4ar?dr. 


Here the plus sign refers to the Bose-Einstein case, the minus sign to the Fermi-Dirac case and \ 
is the de Broglie wave-length corresponding to the temperature T, \=h[2xmkT }-}. 

The formulas above are valid only for high temperatures or small densities, namely, for V> NV)’. 
For the Fermi-Dirac gas, Wigner and Seitz? have given the distribution function for the opposite 
limiting case VN)* (strong degeneracy), where the distribution becomes independent of tem- 
perature: 





N 2rkor cos 2xkor —sin 2rkor}? 
—i1-—9 4rr*dr. 
V 2rkor 


Here ko is the wave number per cm at the top of the Fermi distribution i.e. (47/3)keV=N. 

A general scheme of the statistical thermodynamics of the density matrix has recently been given 
by Kofink.’ 

In this paper we will explicitly calculate this correlation effect for a Bose-Einstein gas at all tem- 
peratures. One may expect that the peculiar condensation phenomenon, by which the degeneracy of 
the B. E. gas is distinguished, will have some characteristic consequence on the distance correlations. 
We shall show that, under certain conditions in the condensed state, the density as seen from any one 
molecule appears, over a rather extended region around it, to be larger than the true density. This 
peculiar phenomenon seems to have direct relation to the singular fact that liquid helium has a 
negative coefficient of thermal expansion below the so-called \-point (2.186°K).4 While there are 
still other indications which support the view that the A transition of liquid helium corresponds in some 
way to the condensation phenomenon of a Bose-Einstein gas,*:** we will not, in this paper, attempt 
to substantiate this idea in greater detail and shall restrict ourselves to only qualitatively discussing 
the possible connection between the condensation phenomenon of the B.E. statistics and the 
anomalous sign of the coefficient of thermal expansion of liquid He. 


1G. E. Uhlenbeck and L. Gropper, Phys. Rev. 41, 79 (1932). 

2 E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933). 

3 W. Kofink, Ann. d. Physik 28, 264 (1937). 

4W. H. and A. P. Keesom, Physica 1, 128 (1933). 

5F. London, Nature 141, 643 (1938); Phys. Rev. 54, 947 (1938). 

5¢ |. Tisza, Nature 141, 913 (1938); Comptes Rendus 207, 1187 (1938). 
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2. DISTANCE DISTRIBUTION IN AN IDEAL BOSE-EINSTEIN GAS 


The symmetric eigenfunction ¥ of N identical B.E. particles, of which n, are assumed to be in the 
quantum state p with the single particle eigenfunction, y,, can be written in the form 


(ti, fe-++ 5 Z 


Wo, m- =[1/(N!no!n!n2!- -+)#] - Wo(Ti)Wo(t2) - + -Wo(tno)Wiltno+1)-+-Yiltnotm)---. (1) 


Here >°q@ means the sum over all N! permutations of the particle numbers. The factor outside is 
chosen so that the eigenfunction VW is normalized to the value 1 if the single particle eigenfunctions 
y, are normalized to the same value. 

Now let the y, be the eigenfunction of a free particle of the momentum k,h in a box of volume V 


with periodic boundary conditions 
y,= V—'e***k, -r (2) 


and let us calculate the probability P(ri, r2) of the configuration (ri, r2) of the two particles 1 and 2. 
Because of the symmetry of ¥ it is immaterial which ones are chosen: 








P(t, 12) = f soe f ogee ee > ne 
One obtains: 
P(r, 5, POO a . > n,n, cos 2x(k,—k,)-(r1—f2) }. (3) 
v2| N(N—1) p< 


Here it must oe observed that the >>,<. has to be taken so as to contain each pair of states p¥o just 
once. 

Let us now calculate the mean value of P(r., r2) at the temperature J. It is well known® that 
(NpNe)w = Niptie for pX%o. It would further be convenient if we could sum over p and o independently. 
We are allowed to do so in Eq. (3) if we subsequently subtract the terms p=o and divide by 2 for 
po. Thus we may write 


1 1 
(P(r, n= lol D 2 Ripe cos 2x(k,—k,)-(r1—f2) — ~ tip ‘|. 
Since the distribution in temperature equilibrium is isotropic we may first take the average over 
the orientation of k, and k, for fixed values of k,=|k,| and k,=|k,| separately. In this way the 
mean value over the angles: 


sin 2rk,|ri—fe2| sin 2rk,|r1—F 
(cos 24(k, —k,) - (t1—12))w= wk, | 01 | wk, |ti—T2| 





2xk,|t1—fe| 2rk,|ti1—f2| 


is easily calculated and one obtains: 


1 1 _ sin 2rk,r\? os 
Peon dm= Fl tea (E Gepy) —E*'] 


This expression depends only on the relative distance r= |r;—r2| of the two particles, not on their 
center of gravity. Integrating over the latter, we obtain a factor V. Finally, since everything is sym- 
metric in all particles, we may multiply by N—1 (or as well by N, since we assume N>>1) in order to 
obtain the mean number of particles per volume unit or the particle density mie at the distance r of 


any given particle: 
Di)= alt+s (= sin ——) -5 ‘|| (4) 
ie, a, A Gad = od 


6 FE. Schrédinger, Physik. Zeits. 27, 95 (1926). 
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We will now calculate the density D(r) by assuming for 7, the B.E. distribution 


n,=Lexp (Bep+a)—1]-'=Lexp (wA*k,? +a) —1 (5) 
where B= (kT) =2am)?/h? and €,=k,7h?/2m. The quantity a@ is defined by the condition: 
N=> ’,=> [exp (x\*k,?+a)—1]-'. (6) 
p p 


In this relation (6) we may replace the sum over p by the integral f - +» +49 Vk*dk. However, care 
0 


must be taken as to the term referring to the ground state ko=0, as this has to be considered 
separately. This state would obtain the incorrect statistical weight 0 instead of 1, if it were re- 
garded as included in the integral.’ For the other low states the discrepancy between the true 
value and the asymptotic formula, 4% Vk*dk, for the number of quantum states in the wave number 
interval (k, k+dk) is of no consequence. Thus we may write for (6) 


N=[e*—1}'+ J [exp (wd2k?+a) — 1-42 Vk2dk 
0 





: (6!) 
=[e*—1} + V/A & I4e-!«, 
l=1 
The term [e*—1]"' refers to the contribution iio of the particles in the ground state ko=0. 
Similarly, by separately taking account of the ground state, we may now calculate the sum 
sin 2rk,r . sin 2rkr 
ni,—_————_ = not f [exp (x\*k? +a) —1]}7! 4a Vk2dk 
, 2rk pr 0 2rkr 
a7 = ° : 
=fhy+— > aod | exp (—l\?k*) sin (2rkr)kdk 
Yr t=1 0 
=fio + V/d3 > lt exp (—la—-ar?/ld?) (7) 
l=1 
and the sum: 
i ata f [exp (xA*k?+a) —1 }-*44 Vke2dk 
p 0 
=ner?+V/r3 > (-—l-e—'« (8) 
l=1 


or, because of (6’), « 
= fig? + V/r*3 & [-4e-'*#— (N—fio). 


l=1 
We introduce into (4) the values given by (8) and (7) and obtain, with D,= N/V and x=ar’/d’, 
1 @ 2 oa 
D(x) =De| 1+ (net V/x3 > teteatt) — fig? — ion tN—- V/x3 - x rel}. (9) 
l=1 l=1 


When determining a (or ip =[e*—1]-") by Eq. (6’) we have to distinguish between two essenti- 
ally cifferent cases, depending on whether 


NdP/V <> I-§=2.612 


l=1 
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or NiX*/V>2.612. In terms of the ‘‘condensation temperature” To, 


h? N t 
T= ( ) (10) 
2armk \ 2.612 V 


the two cases correspond to the alternative 7 >7> and T <7», respectively. 

(a) If T>T>, then Eq. (6’) is satisfied by a of the order 1. In this case the term [e*—1 ]-!=%p in 
(6’) is also of the order of 1, i.e., very small compared with N and may be neglected and we may write 
the Eq. (6) simply: 





N=V/nd* > i-te-« 
l=1 
or, because of (10) 


1 
T)/T)i=——SLI-'e—'« for T >To. 6a 
TITY rT >To. (6a) 
Accordingly in (9) the quantity %» may, for T 27>, be completely ignored. 

(b) If, however, T<T>o, then (6) requires a to be extremely small, viz., of the order 1/N. In this 
case the sum of (6’) can be replaced by its value for a=0, viz., 


a) 


> -3=2.612, 


ind 
and fio = (e*—1)—! may be replaced by 1/a. One obtains for (6’) 
N=1/at2.612V/\ 
or, because of (10), 
1/a=no= N[1—(T/T>)!] for T< To, (6b) 


i.e., for T7<To, a finite fraction of all particles is in the single state ky=0. 
Considering (9) we notice that the series 


@ 
> l-te-'« 
l=1 


diverges for a=0. However, we need the series rather for very small a(~1/N), not for a=0. For such 
small values of a we may replace the >>; by an integral over (Ja) and calculate its asymptotic value: 


Z. Helena f (la)—4e-Cd (la) = (4/a)* (for aK1). 
0 . 


l=1 


Taking a from (6b) we obtain: 


> -he-'# = (x N) 31 —(T/To)? }} for T<To. (11) 
ll 
For temperatures T>T> the quantity a is of the order of one and consequently the value of the 
series }-]/-4e—'¢ itself is of the order of one. 
The term 
(V/N2n*) > “tet 
l=1 
in (9) will, therefore, be of the order N-? for T<T> and even of the order N— for T> To; in each case 
it is negligible compared to 1 and may be omitted in (8) and in (9). 
Let us finally define the function 
2.612V 
eV, T)= V =(T/To)'=V/Vo (12) 
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where Vo=(N/2.612)(h?/2rmkT)? is the volume at which, at the temperature 7, Bose-Einstein 
condensation appears; further put 
1 wo 
, x) =— ete! 13 
Q(a, x) rr - (13) 


Omitting in Eq. (9) the terms found to be negligible we now summarize the foregoing remarks and 
obtain for D(r): 
D.[1+ ¢°Q?(a, rr?/d?) ] for g 21 
D(r) = (14) 
D.[1+(1—¢+¢Q(0, rr*/d?))?—(1—¢)*] for g <1. 

In Fig. 1 the function D(r)/D,, is drawn for different temperatures and for a fixed value of V. 
The unit of length chosen is \(T) =h[l2armkT, |-= (2.612 V/N)'. Conceived as a function of g and 
 (i.e., of V and 7) for fixed r the function D(r, ¢, \) and its first derivatives are continuous every- 
where, but the second derivatives of D are discontinuous at the transition curve (g=1) for any value 
of r. 

The asymptotic behavior of D(r) for r>d is very different in the cases T>7T>) and T<T>. In either 
case we may obtain asymptotic formulas if we notice that, for x>1, the lowest terms of the series for 
Q(a, x) (Eq. (13)) with an index /<x are smaller than e-7/' and may be neglected. Further, the 
subsequent terms with an index /~x and />x vary so smoothly with / that we may replace 


Sty fina 
lal 0 
In this way we obtain, for x>1: 


1 
2.612 





1 ea) 
Ola, 2>1) = J [-le-!e—z1 1] = (/x)! exp [—2(ax)!}. (15) 


For a>0O or T>Tp>, therefore, D(r)—D, disappears exponentially for large r/X namely like 
1/r exp [—4r/A(ra)*]. For T=T o, however, D(r)—D., behaves asymptotically like 1/r*; finally, 
for T<To, D(r)—D,. behaves like 1/r. 

In the case 7 >To, the excess of particles around a chosen one is finite, i.e., independent of the 
total number NV. The number of these excess particles is easily calculated to be 


¢/2.612./8 for T>To, 
2re?/2.6127(¢—1) —1 for 0< T—To/To<1. 


In the case T= T> this excess is of the order N! and for T<Tp> it is of the order N?. 


3.- VAN DER WAALS EXCHANGE INTERACTION 


The distance correlations which we discussed in the preceding paragraph and the departure of 
the ideal B.E. gas from the behavior of an ordinary ideal gas are both only different sides of the same 
matter and thus far our present investigation can, of course, not add anything new to the thermo- 
dynamics of an ideal B.E. gas. These correlations may, however, make themselves felt if there is 
an intermolecular force field insofar as they may give a distinguished weight to those distances 
where the Van der Waals field is noticeable. This will affect the resultant interaction in a very 
different way in the cases of the degenerate and of the non-degenerate state. 

L. I. Schiff? has made a very general investigation of the degeneracy of a non-ideal B.E. gas in 
which he develops the grand partition function of the system in subsequent steps by considering, 
in close analogy to the method of Ursell, the consecutive multiple clusters of particles separately. 


m on Goat) Phys. Rev. 59, 751, 758, 927 (1941); ibid. 60, 362 (1941); see also W. E. Lamb jr. and A. Nordsiek, ibid. 
q 1941). 














208 F. LONDON 


This method is entirely rigorous and should, at least in principle, contain the answer toevery thing that 
could be asked about the influence of interaction forces on the B. E. condensation. However, it is 
scarcely possible actually to calculate a higher approximation than that which corresponds to the 
second virial coefficient (binary clusters). It is obvious, however, that cooperative phenomena, like 
a condensation, cannot satisfactorily be discussed on the basis of a consideration of binary clusters 
only. Their investigation certainly requires the consideration of the whole of all virial coefficients 
and their convergence properties similar to that which has been developed in the theory of con- 
densation by Mayer and others.*® 

Since we must probably forego any hope, for the present at least, of finding a higher approxi- 
mation of the Ursell-Schiff development we might try to choose an approach which, though less 
rigorous, allows one to account more directly for the simultaneous interaction of many neighboring 
particles. There is the method, which has been applied by Heisenberg and more specifically by Bloch® 
in the theory of ferromagnetism, which offers itself rather immediately to the present problem. 
While there is no question about the fact that this method is quite crude and, still worse, that it is 
rather difficult to say in which sense it can be considered as an approximation, it has, in the case of 
ferromagnetism at least, served to sketch a very simple picture of what certainly is the essential 
mechanism for the establishment of that cooperative phenomenon. This method consists in taking, 
as a basis, the first step of the quantum mechanical perturbation development in the following sense : 


any state of the ensemble is represented by an eigenfunction of ‘‘zero order,”’ built up, like our 
(ti, t2, ++ +ty) , , : , ‘ —— 
function Wno, m1, Ks in Eq. (1), by the eigenfunctions of isolated single particles. This eigen- 


function W itself does not account for the interaction forces except implicitly insofar as it has 
been given the ‘right symmetry” so that it might be thought of as continuously deforming into 
the true eigenfunction if the interaction were continuously introduced. Thus the state of the 
ensemble can still be described by the numbers u, giving the distribution of particles in the single 
particle state p and by the symmetry character of ¥. The presence of interaction forces is, however, 
explicitly considered by the energy attributed to this state, in that for this energy is taken the 
quantum mechanical mean value of the total energy H (including the interaction energy) in the 
configuration Wno, 21, n2---n,---. This mean value is given by 


E(no, N1, N2°**Ny*° )= f ¥en, S4°°**Ry** -(H, WVno, m-- ‘np |d71, dr2: . -dtn (16) 


which is, in fact, the so-called first-order approximation for the energy. One may now develop the 
statistics for the mo, m1: --m,- ++ in the usual way by utilizing the above approximation for the energy. 
Let us assume as Hamilton function 
1 ; 
H=—¥ pt+¥ Wri) (17) 
2m x i<k 

where W(r,) is the potential of the intermolecular forces, assumed as additive central forces, between 
the particles 7 and k. The approximate energy value (16) of a state characterized by the distribution 
No, M1°*-N,--+, 1.e., by the W-function (1), is given by 


E(no, M1, °**Np** -)=>on,€,+ U(no, M1, °° "Mp, ** -) 
where €,=k,*h?/2m and 
(ti, fe-+ +r 


) 
U(no, ma? *** -) = z. W(rix) | V no, ni- ily tee |°?dridt2° . -dty. 
i<k 
Because of the symmetry of H and ¥ in all particles all }N(N—1) terms of the sum > j<; have the 
8 J. E. Mayer, and S. F. Harrison, J. Chem. Phys. 6, 87, 101 (1938); M. Born and K. Fuchs, Proc. Roy. Soc. A166, 


391 (1938); B. Kahn and G. E. Uhlenbeck, Physica 5, 399 (1938). 
*F. Bloch, Zeits. F. Physik 57, 546 (1929). 
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same value. By using (2) and (3) we obtain: 


U(no, 11° me aN f Win {14 — : NpN_ COS 2r(k, —k,)- (rir) }dradra 


N(N—1) p<e 
Therefore we may write: 


N(N-1) 
E(no, N11, eee —————____— 


Wott X nneWoe (18) 


pHa 


Np ‘Y= n,e,+ 


where 


4nr 7°” 
w= [ W(r)r2dr (19). 
V Jo 
is the ‘‘classical’’ mean value of the van der Waals energy per pair, whereas the quantities 


1 
Wee "2 f wer) LW o(ri)Wo*(r1)Wp*(r2)Wo(r2) +o p*(ri)bolti)bo(r2)Wo*(r2) \dridre 
' (20) 
=> fro cos 2r(k,—k,)-rdr 


are the “‘exchange integrals” which account for the correlations of the B.E. statistics. In the Fermi- 
Dirac case they would appear in the energy expression with negative sign or, if there is a spin, with 
both signs depending on the relative spin-orientation, and would give rise, eventually, to ferro- 


magnetism. 


Since we assume the molecular forces to be central forces, we may integrate, in (20), over all 


orientations and obtain: 


4x (” 
W,.=— f Wir) 
V Jo 


4. PRELIMINARY DISCUSSION OF THE 
EQUILIBRIUM DISTRIBUTION 


Let us now consider the distribution for 
thermal equilibrium taking the energy (18) as a 
basis. The condition for the most probable dis- 


tribution, wo, %1, ---%,---, leads to the fol- 
lowing system of simultaneous equations: 
wa 
In —-—=6aE/an, +a 





p 


or with (18): 


A,+1 
kT In =e,+)> ioWe,takT. (21) 


Lp op 





It seems rather difficult to solve the system of 
Eqs. (21) in full generality. We may, however, 
gather some information from a discussion of 
their solution at 7=0°K, which solution can be 
found easily. We may presume that the W,, 


which, in this case, come in question belong to . 


very small momenta k, and k, so that |k,—k,|r 
will be very small compared with 1 for all values 
of r for which W(r) is noticeably different from 


sin 2r|k,—k.|r 





rdr. (20’) 


2r\k,—ke\r 





zero. Indeed we may verify this assumption once 
we have found the solution. Then we may replace 
the sin in (20’) by its argument and obtain 
Woe ~ Wo for all p and ¢ which come in question. 
Further it appears reasonable to presume that 
Wo is negative, Wo=—|Wo|, since it is an 
average taken over an, in general, attractive 
field. Thus we may write, for sufficiently low 
temperature, instead of (21) simply 
eT in =e =e,—|Wo| X a.takT 


i, exp 
=e,—|Wo|(N—A,)+akT 

or with the notation n= —akT+N! Wo 
ni,+1 


ty 





kT In 





=e,—utn,| Wo) (21’) 


from which equations, in the limit 7-0, the 
quantities %,—>n,° and w—y° are easily deter- 
mined: 
(u°—e,) 
Sle for €,< yu" 
| Wo} (22) 


0 for €, >p°. 
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Fic. 1. Density of an ideal Bose-Einstein gas around a 
given molecule for different temperatures and for fixed 
volume. 


The condition >>, ”,°=N determines the quan- 


tity p°: 
-(-) (2 N| Wo! 2/5 
Nam \Be iV ; 


The energy at absolute zero of the distribution 
(22), (23) can be evaluated easily : 


Exin= 1,°€,= (3/7) Np, 
Eexch= —3/ W,| He s,°n,° 


po 


= —3|Wo|X ,°(N—n,°) 


p 


=—}|Wo|(N*—¥ (n,")*) 





(23) 
(24a) 


(24b) 


= —2|Wo| N?+(2/7)Nu° 


1 
Eaiass = i Wo | : (24c) 


In total we obtain for the energy (18) of the 

distribution (22) (23): 
1 

B= — N*| Wel +(5/D)Nv?. (24) 

The distribution (22) appears very different 

from that of the ideal B.E. gas: Instead of being 

concentrated in a single quantum state the 

molecules are spread, somewhat similarly to a 

Fermi distribution, over many quantum states 

according to the condition e, <y®, which condition 


together with (23), yields the value k,, for the 
wave number of the highest occupied state, 


LONDON 





namely, 





2m|Wo| 15N\ 
) . (25) 


1 

=— 0)i— 
—e ( kh? 8nV 

However, this spreading over many states is 
of sub-microscopic extension so to speak, and 
has no noticeable effect, for example, on the 
kinetic energy. In order to obtain an estimate of 
the orders of magnitude in question let us notice 
that N(h?/2m)(15N/8rV)! is of the order of 
the thermal energy per mole at the transition 
temperature, say RO, and that N?| Wo| also can 
be presumed to be of the order of a small thermal 
energy per mole, say RO:, since it corresponds to 
the contribution of the van der Waals cohesion 
to the evaporation heat. With these notations 
the increase of kinetic energy (24a) due to the 
spreading over the states k,<k,, can be written 
in the form 


Exin =_ (3/7) N-?4 RO 3/56 ,2/5, 


This shows the kinetic energy is of the order 
N-'5 times a measurable thermal energy, i.e., 
negligible. 

The number of occupied states can be esti- 
mated easily in the same way. It is given by 


4r 5 (=?) 3/5 
—k,,? V=—N?2/5{ — rw N25, 
3 2 0; 


The number of molecules in one of the occupied 
states must accordingly be of the order N*/>, For 
example, from (21) and (22) we obtain for the 
number of molecules in the lowest state: 


no? = N2!5(Q,/O2)3!9, (26) 


While, with respect to the kinetic energy, the 
transition to the new distribution (22) involves 
no noticeable change, it has a considerable effect 
on the exchange energy. If all molecules were in 
a single state the exchange energy would be 
zero, whereas the apparently minute spread over 
the N?/> states represents a considerably more 
stable arrangement, since it has an exchange 
energy of about —(N?/2)| Wo]. 

For a temperature 7 >0 the exchange integrals 
of the excited states, W,,, will appear in (18), 
which quantities then can no longer be replaced 
by Wo. Since | Wye| <|Wo|, the exchange inter- 
action will have the tendency to favor the B.E. 
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condensation once the latter is initiated and this 
might lead to a more pronounced discontinuity 
at the condensation point than that of an ideal 
B.E. gas. 


5. DISTANCE CORRELATIONS IN THE CASE OF 
WEAK EXCHANGE INTERACTION 


We will now calculate the distance correlations 
for the modified distribution (22), which may be 
written in the form 


. awe for kp< Rm 
9 = 


(22’) 
0 for k,>km 


where k,, and n° are given by (25) and (26). If we 
introduce this distribution into (4), we may 
easily calculate the summations in question and 
obtain 


D(r) =D. [(3—(2akmr)?) sin (27kmr) 


{ts 
3(2tkm?) (2akmr) |? “| (27) 
TR?) COS (2tRm?r) ani" 


The term —(4/7)no°/N, which corresponds to 
the term (1/N?)>°, ,? in Eq. (4), is negligible 
compared with 1 since it is of the order N~?/5 
(see Eq. (26)). For 27kn»r<K1 we may develop 
the sin and cos in (27) and obtain: 


(2tkmr)? 4 


D(r)~D,42- —(2ek7)*—->- 
) 7 tm” 


for rK(2rkm)—. 


Seen from a given molecule it appears as if in 
the neighborhood within the radius r ~~ (27k) 
= N'5(V/N)"/3(61/02)?’9, i.e., of the order of N¥/5 
mean molecular distances, the density of the gas 
has just the double of its true value. The number 
of the excess molecules amounts to (4/7) mo° ~ N?/5. 
Outside of this region of anomalous density the 
excess must be, of course, compensated by a cor- 
responding decrease in density below its normal 
value. This decrease is accounted for by the 
small term —(4/7)n0°/N in (27), which, in this 
balance, therefore must not be ignored. 

It is a characteristic feature of the distribution 
(22) that its >, 2,” is of a smaller order than N? 
(namely N%/5 see Eq. (24b)) whereas for the 
ideal B.E. distribution in the case of degeneracy 
this sum is of the order N? (Eq. (8)). For higher 
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temperatures the distribution will certainly not 
be more concentrated; hence the >>, m,? can 
only be smaller than its value at 0°K. Thus, it 
is safe to conclude that >>, 1, is always to be 
considered as small compared with N? as soon 
as there is a however feeble interaction with a 
negative exchange integral. 

Accordingly, when we now consider the 
distance distribution for T>0°K, the term in 
Eq. (4) containing (1/N?)>°, n,? can certainly be 
regarded as negligibly small and may be omitted. 
On the other side the sum (7) will have to be 
recalculated. The integral over the states with 
an energy of the order of a thermal excitation 
will be about the same as before. The term fig in 
(7), however, will now have to be replaced by 
a sum (or integral) over the states of a ‘“‘sub- 
microscopic’ energy interval containing sin 
factors with various wave numbers of the 
order ~k, (Eq. (25)). For distances r&k,,~! 
these sin factors can again be replaced by their 
arguments and in this case, the sum (7) yields 
the same value as before, again with a finite 
fraction of N, corresponding to fig in (7), set 
apart from the integral. For r>k,,-!, however, 
these sin terms will cancel each other by inter- 
ference as in (27), i.e., for r>>k,»~! the term cor- 
responding to ip in (7) is to be omitted. Ob- 
serving the preceding remarks we obtain, instead 
of (14), for g=1: 


D..[1+(1—¢+¢Q(0, mr?/d?))?] 


Dir) for rKRy! (28a) 

Y)= 

' D.[i+(1—¢)2] for \Kr<Kkm-! (28b) 
D, for r>k»-!. (28c) 


For g=1 i.e. for the uncondensed state, the dis- 
tribution is the same as before: 


D(r) =D.[1+¢0%(a, wr?/d2)]. (28d) 


In Fig. 2 the distribution (28) has been drawn 
for different values of JT. We see that, in the 
condensed state, apart from a small region within 
a radius \, an observer fixed on a given molecule 
has the impression of being in a medium, not of 
the ‘true’ density D,, but rather of an ‘‘ap- 
parent” density D.[1+(1—¢)?]. Only in a 
distance of the order kn-'~N*(V/N)"/3 does 
this apparent density decrease to the true value 
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Fic. 2. Density of a Bose-Einstein gas around a given molecule for the limit of a weak van der Waals exchange interaction. 


D., (or rather to the slightly smaller value 
Dal1—(1/N*) din, )). 

The dualism of an ‘‘apparent” and a “‘true”’ 
density is not so difficult to picture as it might 
appear on first view. A simple model which shows 
distance correlations of the kind in question is 
represented by the molecular distribution in a 
piece of Swiss cheese full of holes. In such a cheese 
infinitely many more molecules are located in the 
interior of the massive substance, than are found 
on the surface of a hole. Accordingly the mean 
density around an average molecule is the density 
of a piece of hole-free cheese. This, howeve., 
holds only within a range which is given by the 
mean distance of a molecule from the next hole. 
Beyond this distance the mean density relative 
to a central molecule decreases to the “true 
value,” as it is given by averaging over the holes 
and the massive substance indiscriminately. 

It is significant that in our case the increase in 
density of the hole-free region, D.(1—¢)’, is 
independent of the strength of the intermolecular 
energy, but rather determined by the degree of 
degeneracy, as a pure statistical effect. Only the 
number and the size of the holes depend on the 
strength of the van der Waals field, which, there- 
fore, at least in the limit of very weak interaction, 
is reduced to the role of essentially a surface 
effect. 


CONCLUSION 


Comparing the two distribution functions, 
Fig. 1 and Fig. 2, we notice that there is, in the 





condensed state, a peculiar instability. Figure 2 
represents the limiting case of a very small 
negative exchange integral. It can further be 
shown that Fig. 1, which we calculated to be the 
distance distribution of an ideal B.E. gas, more- 
over represents also the limiting case of a very 
small positive exchange integral. At the transition 
curve (g=1), the discontinuities of D, conceived 
as a function of ¢, are of the same character for 
the two cases, W20 and W<0: D and dD/d¢ 
are continuous, but 0?D/d¢? is discontinuous for 
g=1 for any value of r (except for r=0 and 
W. <0 where D=2 independently of ¢). 

The increase of D(r) (Eq. (28)) to the value 2 
with r—0 is only a mathematical formality with- 
out any physical significance. Actually the inter- 
molecular repulsive forces prevent the molecules 
from coming nearer than a certain minimum 
distance, within which D(r) will in reality be 
practically zero. These repulsive forces are, of 
course, supposed to be included in the interaction 
potential W(r); however, they do not obtain 
explicit expression in the first-order approxima- 
tion (18) which has been considered exclusively 
so far. The higher approximations of the per- 
turbation method will, of course, account for the 
actual impenetrability of the molecules. 

On the other hand it is quite clear that our 
results, which are based on the assumption of 
wave functions of force free molecules, can cer- 
tainly not be trusted as valid in distances beyond 
the mean molecular free path. 

In view of these two limitations—with respect 
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to large as well as small distances—it might ap- 
pear that our results cannot have much relation 
to reality. There is, however, some chance that, 
notwithstanding these limitations, some residue 
of the distance correlations may be realized by 
helium—the only substance which is at least 
liquid for the temperature and density for which 
B.E. condensation should occur. Tisza®* has di- 
rected attention to the peculiarity of normal liquid 
helium (above the A-point) having a positive 
temperature coefficient of viscosity (dn/dT >0); 
this is distinctly different from the behavior of 
the viscosity of any other liquid (which all have 
a 0n/OT <0), but quite analogous to the viscosity 
of a gas. The above fact as well as the relatively 
very open structure of liquid helium as indicated 
by its x-ray analysis" have been taken to support 
the view® that it might be justified to consider 
liquid helium (below and above the )-point) 
somewhat more similar to a gas than one would 
be inclined to believe reasonable for liquids in 
general; thus one might try to interpret the vis- 
cosity above the A-point as due to a typical gas 
viscosity mechanism with momentum transfer 
over a mean free path rather than of the liquid- 
viscosity type (with momentum transfer over po- 
tential barriers) and, in particular, interpret the 
peculiar state of He below the )-point as corre- 
sponding to the condensed phase of an ideal 
B.E. gas. 

If then it appears reasonable to speak, with 
all due reserve, of a mean free path of the 
molecules in liquid helium it might be worth 
while to discuss, at least qualitatively, what 
implications would follow from assuming that 
there actually exists, within the range of the 
mean free path around every single molecule, a 
tendency toward such an increase of density as 
it is described by our formula (28b) for a de- 
generating B.E. gas in medium distances from a 
given molecule. 

It is true, the repulsion between the neighbors 
of the central molecule will actually exert a re- 
sistance to such an enormous increase of density 
and, in addition, the repulsion due to the zero- 


1 W. H. Keesom and K. W. Taconis, Physica 5, 270 
(1938); F. London, J. Phys. Chem. 43, 49 (1939). 





point energy coming from the quantization of the 
mean free path" will also act in the same direction. 

However, this opposition will not necessarily 
have a prohibitive character. It can rather be 
expected that a compromise between the two 
tendencies will be reached: For a fixed volume, 
the crowding around each central molecule will 
not be of quite as large an amount and of as wide 
a range as given by (28b); but there will be at 
least some crowding. On the other hand, in this 
somewhat crowded condition the repulsion 
between the neighbors and the action of the zero- 
point motion will have an increased effectiveness 
and this will be felt as an additional pressure 
necessary for maintaining the volume. This 
pressure will increase with decreasing tempera- 
ture, namely with progressing degeneracy. Thus 
(dp/dT)y will be negative for T<To. 

Or to say it differently inasmuch as it is the 
density relative to any one molecule which is essen- 
tially significant for the operation of the inter- 
molecular forces, the system will tend to keep the 
“apparent” density about constant at a value 
near that realized above the condensation point; 
to accomplish this the true density correspond- 
ingly has to decrease. For constant external pres- 
sure, therefore, the system will, with progressing 
degeneracy, expand to a greater volume. Thus 
also (0V/dT), will be negative for T<T>. 

It is certainly very unusual if a substance 
shows a negative coefficient of thermal expansion. 
Since (8V/dT),= —(0S/dp)r, this property im- 
plies that by isothermally compressing the sub- 
stance in question, entropy can be increased or 
disorder generated. On first view one should 
think that by compressing a system of spheri- 
cally symmetrical molecules it could only be 
brought into an arrangement of higher order. 
In fact liquid He below the \-point is one of the 
few examples where the anomalous (negative) 
sign of (9V/dT), and (0p/dT)y actually has been 
found. It seems hardly possible to think that any 
mechanism, other than the one we have dis- 
cussed above, could be made responsible for this 
singular behavior of liquid helium. 


1 F, London, Proc. Roy. Soc. A155, 576 (1936). 
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On the Structure of Mo(CN),;-* 


GrtuLio RACAH 
The Hebrew University, Jerusalem, Palestine 

January 27, 1943 
T was shown by Hultgreen! that no more than six equiv- 
alent cylindrical bond functions may be constructed by 
s—p—d hybridization, and it was generally assumed*? that 
f eigenfunctions are needed to construct eight equivalent 
bonds. Since the restriction to cylindrical bond functions 
was only based on the feeling that other bond functions 
would have a considerably smaller strength, it seemed 
worth while to us to show that eight equivalent non- 
cylindrical bond functions may be constructed by s—p—d 
hybridization, and to calculate the strength that may be 

reached by them. 

It is easy to see with the help of the group-theory meth- 
ods used by Mulliken‘ and Van Vleck? that this construc- 
tion is possible only if the symmetry group of the bonds is 
the point group D,, and that the bond functions are 


(p2+dz4,+d2,+scosp+d,sing) /8} 

+[(p2+d,2)cosa+ (py+d,2)sina ]/2 
and the seven others which are obtained from it by the 
symmetry operations of D,. The parameters a and £6 are 
arbitrary; if we fix them by the condition of giving the 
maximal strength to the bonds, we obtain a=22°30' and 
8=18°; the bonds are then directed to the corners of a 
square antiprism, make an angle of 60°54’ with the tetra- 
gonal axis and have the strength 2.9886, while the best 
s—p-—d bond has the strength 3. The restriction to cylin- 
drical bond functions on ground of strength considerations 
is therefore not justified. 

The obtained structure, however, does not agree with 
the experimental results of Hoard and Nordsieck,® which 
found that the eight bonds of Mo(CN)s~ are not equiva- 
lent, but split into two sets, each consisting of four approxi- 
mately equivalent bonds and having very nearly the 
symmetry Dea; the bonds of the first set make angles of 
about 34° with the z axis, those of the second set angles of 
about 73°. 

This structure also may be obtained by s—p—d hy- 
bridization, assuming for the first set the bond function 


(ssina+d,cosa+d,4,sinB+ p,cosp) /2 
+ (pzsiny+d,.cosy) /24 


and the three others which are obtained from it by the 
symmetry operations of D2z, and assuming for the second 
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set the same expressions with the parameters a, 8, and y 
increased by 90°. The sum of the strengths reaches its 
maximum when a=40°13’, B= 19°35’, y =30°57’; for these 
values of the parameters the bonds of the first set make an 
angle of 34°33’ with the z axis and have the strength 2.9954, 
those of the second set make an angle of 72°47’ and have 
the strength 2.9676. 

The strengths are also in this case very satisfactory and 
the angles agree with the experimental results. 

1R. Hultgreen, Phys. Rev. 40, 891 (1932). 

2 J. H. Van Vleck, J. Chem. Phys. 3, 803 (1935). 

3W. G. Penney and J. S. Anderson, Trans. Faraday Soc. 33, 1363 
OTR S. Mulliken, Phys. Rev. 43, 765 (1933). 


5 J. L. Hoard and H. H. Nordsieck, J. Am. Chem. Soc. 61, 2853 
(1939). 





The Inhibitory Effect of Iodine on the Photolysis 
of Gaseous Hydrogen Iodide 


RIcHARD A. OGG, JR. AND RUSSELL R. WILLIAMS, JR. 


Department of Chemistry, Stanford University, 
Stanford University, California 


March 1, 1943 


HILE seeking a precision actinometric method for 
study of another photochemical reaction, we have 
reinvestigated the supposedly classical photolysis of gase- 
ous hydrogen iodide. Monochromatic light of \2537A was 
used, and the reaction was followed by measurement of the 
quantity of hydrogen produced. The partial pressure of 
hydrogen iodide was varied from some 50-150 mm Hg, and 
that of iodine from 0-20 mm Hg. Twelve separate experi- 
ments at 108°C are in excellent agreement with the widely 
accepted mechanism (analogous to that for photolysis of 
hydrogen bromide): 


HI+/v—>H +I, 
H+HI—H2+I, ke 
H+I.—HI-+I, ks 
21+M—I.+M. 


The ratio o. rate constants k3/ke was found to have the 
value (at 108°C) of 4.9+0.3. Over a temperature range of 
some 50°C there appeared-to be no appreciable change in 
this ratio. ; 

The only previous quantitative investigation of this 
iodine inhibition appears to be that of Bonhoeffer and 
Farkas,! who estimated the ratio k3/k2 to be about 100— 
these authors claiming that 15 percent of iodine “‘com- 
pletely” inhibits the photolysis of hydrogen iodide. This 
figure of 100 has been rather widely used in interpretation 
of related reaction rate experiments.? However, the experi- 
ments of Bonhoeffer and Farkas were conducted at ex- 
tremely small gas pressures, and there seems little doubt 
that the secondary reactions occurred principally on the 
walls of the vessel. The value of k3/k2 reported here cer- 
tainly corresponds to the true homogeneous reactions. 

The ratio k3/k2 for the analogous reactions of Cle and 
HCI is known qualitatively to be rather large—at least 
100.2 For the analogous reactions of Bre and HBr the ratio 
is accurately known to be 8.4.3 The above value of 4.9 for 
the reactions of Iz and HI thus falls into a logical sequence, 
whereas the previously estimated value seemed highly 
anomalous. 
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The practical independence of the ratio k3/k2 on tem- 
perature indicates that the above reactions 2 and 3 have 
practically identical activation energies. However, the pre- 
vious estimate? of zero for the activation energy of reaction 
3 was based upon the erroneous estimate of k3/k2 as 100. 
It seems probable that reaction 3 must have a finite, 
although undoubtedly small, activation energy. 
qoash Bonhoeffer and L. Farkas, Zeits. f. physik. Chemie 132, 235 


2 J. C. Morris and R. N. Pease, J. Chem. Phys. 3, 796 (1935). 
3M. Bodenstein and G. Jung, Zeits. f. physik. Chemie 121, 127 (1926). 





Energy Relations in Film Penetration 


Jack SCHUBERT AND G. E. Boyp 
George Herbert Jones Laboratory, University of Chicago, Chicago, Illinois 
March 10, 1943 


INCE the energy relations in film penetration’? are 
largely unknown it has appeared worth while to in- 
vestigate these relations for an insoluble film of cetyl 
alcohol penetrated by sodium cetyl sulfate. 
Information on the composition of a mixed film in equi- 
librium with a cetyl alcohol crystal may be obtained by 
use of the Gibbs equation: 


dye= —T dui —T2du2—T sdys, (1) 


where the subscripts designate the species water, sodium 
cetyl sulfate, and cetyl alcohol, respectively; T and y» the 
surface concentrations and chemical potentials; and y. the 
surface tension of the equilibrium film. 

The position of the Gibbs plane may be chosen in a 
manner such that the surface concentration of water, T,, 
is zero. Equation (1) thus becomes: 


dy.= —T2%dpu2— T3dys3. (2) 


Additionally, if the underlying solution is saturated with 
cetyl alcohol, as is likely to be the case when spreading 
results from excess crystals, the quantity us may be as- 
sumed constant, resulting in the further specialization 
of (1): 

(Oye/Om2) y= —T2 (Ou2/dme2)y5, (3) 


where mz is the concentration of the sodium alkyl sulfate 
in the bulk phase. 
Since 
du2= RTd |n mof2, 


where fz is the activity coefficient, Eq. (3) becomes: 


Ove) [4 (Gag |: 
(3) 47 ial i dm: * (4) 
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If it be assumed for highly dilute solutions that 
(0 In fi/Om2),,=9, 


(4) changes to a convenient equation for the determination 
of the adsorption of sodium cetyl sulfate in the film: 


ro = —(m2/RT)(dve/Om2) y;- (5) 


Thus, a method of determining the surface concentration 
of alkyl sulfate is made possible. The information necessary 
is the variation of the surface tension of the mixed film 
with me, the concentration of the alkyl sulfate. 

The free energy change in the penetration process ap- 
pears to be comprised of two parts in the main: (1) the 
reversible work of adsorption of the sodium cetyl sulfate 
from the body of the solution, and (2) the mechanical work 
necessary to compress the cetyl alcohol monolayer suffi- 
ciently to make room for the penetrating molecules. Thus, 


AF=RT In (c°/c*) + | mde, (6) 


where c’ and c* are the surface and bulk concentrations of 
sodium cetyl sulfate, x, the equilibrium film pressure,’ and 
o the molar area of alkyl sulfate in the mixed film resulting 
from penetration. 

The free energy may be obtained also from the equation 


(A Fin) = (AHm)s—(Qm)s, (7) 


where the molar latent heat (Q,,). is given by an application 
of the Clapeyron equation 


(Qm) s= T(m-/dT )o,(om—os) (8) 
and the enthalpy by* 
eee oo O(1re/ T) 
(AHn).=— J. [sary ae |, to. (9) 


Experiments have been completed on the variation of 
the surface tension with temperature and with composition 
of a system consisting of excess cetyl alcohol crystals in 
contact with the surface of a solution of sodium cetyl sul- 
fate. It is planned to publish details and calculations when 
time permits. 

1 (a) J. H. Schulman and A. H. Hughes, Biochem. J. 29, 1236, 1243 
(1935); (b) J. H. Schulman and E. Stenhagen, Proc. Roy. Soc. B126, 
356 (1938); (c) J. H. Schulman and E, K. Rideal, Proc. Roy. Soc. B122, 
29, 46 (1937); (d) F. Sebba and E. K. Rideal, Trans. Faraday Soc. 27, 
273 (1941). 

2 (a) N. K. Adam, F. A. Askew, and K. G. A. Pankhurst, Proc. Roy. 


Soc. A170, 485 (1939); (b) K. G. A. Pankhurst, ibid. A179, 393 (1942). 
3 The equilibrium spreading pressure we is defined by Teo@¥,.—Ye, 


where y, may be either the surface tension of the solution or of the clean 


water at the temperature in question and +e is the surface tension of the 
aqueous or solution surface covered by the insoluble monolayer in 
equilibrium with the crystal. 

4W. D. Harkins, T. F. Young, and G. E. Boyd, J. Chem. Phys. 8, 
954 (1940). 











